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1. Neutrino Astronomy: A Brief History 


Soon after the 1956 observation of the neutrino,! the idea emerged that it repre- 
sented the ideal astronomical messenger. Neutrinos travel from the edge of the 
Universe without absorption and with no deflection by magnetic fields. Having 
essentially no mass and no electric charge, the neutrino is similar to the photon, ex- 
cept for one important attribute: its interactions with matter are extremely feeble. 
So, high-energy neutrinos may reach us unscathed from cosmic distances: from the 
inner neighborhood of black holes and from the nuclear furnaces where cosmic rays 
are born. But, their weak interactions also make cosmic neutrinos very difficult 
to detect. Immense particle detectors are required to collect cosmic neutrinos in 
statistically significant numbers.” By the 1970s, it was clear that a cubic-kilometer 
detector was needed to observe cosmic neutrinos produced in the interactions of cos- 
mic rays with background microwave photons.” Subsequent estimates for observing 
potential cosmic accelerators such as Galactic supernova remnants and gamma-ray 
bursts unfortunately pointed to the same exigent requirement.* © Building a neu- 
trino telescope has been a daunting technical challenge. 

Given the detector’s required size, early efforts concentrated on transforming 
large volumes of natural water into Cherenkov detectors that collect the light pro- 
duced when neutrinos interact with nuclei in or near the detector.’ After a two- 
decade-long effort, building the Deep Underwater Muon and Neutrino Detector 
(DUMAND) in the sea off the main island of Hawaii unfortunately failed.” How- 
ever, DUMAND paved the way for later efforts by pioneering many of the detector 
technologies in use today, and by inspiring the deployment of a smaller instrument 
in Lake Baikal® as well as efforts to commission neutrino telescopes in the Mediter- 
ranean.!°!? These efforts in turn have led towards the construction of KM3NeT. 
But the first telescope on the scale envisaged by the DUMAND collaboration was 
realized instead by transforming a large volume of transparent natural Antarc- 
tic ice into a particle detector, the Antarctic Muon and Neutrino Detector Array 
(AMANDA). In operation beginning in 2000, it represented a proof of concept for 
the kilometer-scale neutrino observatory, IceCube.!*:!4 

Neutrino astronomy has achieved spectacular successes in the past: neutrino de- 
tectors have “seen” the Sun and detected a supernova in the Large Magellanic Cloud 
in 1987. Both observations were of tremendous importance; the former showed that 
neutrinos have mass, opening the first crack in the Standard Model of particle 
physics, and the latter confirmed the basic nuclear physics of the death of stars. 
Fig. 1 illustrates the neutrino energy spectrum covering an enormous range, from 
microwave energies (1071? eV) to 107° eV.'° The figure is a mixture of observations 
and theoretical predictions. At low energy, the neutrino sky is dominated by neutri- 
nos produced in the Big Bang. Nuclear fusion in the sun generates neutrinos with 
keV energy and a spectrum that extends to MeV. At MeV energy, neutrinos are 
produced by supernova explosions; the flux from the 1987 event is shown. At yet 
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Fig. 1. The cosmic-neutrino spectrum. Sources are the Big Bang (CvB), the Sun, supernovae 
(SN), atmospheric neutrinos, active galactic nuclei (AGN) galaxies, and GZK neutrinos. 


higher energies, the figure displays the atmospheric-neutrino flux, up to energies of 
100 TeV which were measured at the Frejus underground laboratory! AMANDA 
experiment.!’ Atmospheric neutrinos are a main player in our story, because they 
are a dominant background for extraterrestrial searches. The flux of atmospheric 
neutrinos falls with increasing energy; events above 100 TeV are rare, leaving even- 
tually a clear field of view for extraterrestrial sources at the highest energies. 

The highest energy neutrinos in Fig. 1 are the decay products of pions produced 
by the interactions of cosmic rays with microwave photons.! Above a threshold of 
~ 4x 101° eV, cosmic rays interact with the microwave background introducing an 
absorption feature in the cosmic-ray flux, the Greisen-Zatsepin-Kuzmin (GZK) cut- 
off. As a consequence, the mean free path of extragalactic cosmic rays propagating 
in the microwave background is limited to roughly 75 megaparsecs, and, therefore, 
the secondary neutrinos are the only probe of the still enigmatic sources at longer 
distances. What they will reveal is a matter of speculation. The calculation of 
the neutrino flux associated with the observed flux of extragalactic cosmic rays is 
straightforward and yields one event per year in a kilometer-scale detector. The 
flux, labeled GZK in Fig. 1, shares the high-energy neutrino sky with neutrinos 
anticipated from gamma-ray bursts and active galactic nuclei (AGN).* © !° 

A population of extragalactic cosmic neutrinos with energies of 60'TeV—1 PeV 
was revealed by the first two years of IceCube data. We will review the present 
status of the observations, the identification of the sources by multimessenger as- 
tronomy, as well as attempts to decipher the phenomenology of the heavenly beam 
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dumps producing neutrinos. Subsequently, we will describe the status of the search 
for Galactic neutrino sources and conclude with brief discussions of other uses of 
neutrino telescopes. 


2. Rationale for the Construction of Kilometer-Scale Neutrino De- 
tectors 


The construction of kilometer-scale neutrino detectors was primarily motivated by 
the prospect of detecting neutrinos associated with the sources of high-energy cosmic 
rays. Cosmic accelerators produce particles with energies in excess of 100 EeV; we 
still do not know where or how;?° see Fig. 2*. The bulk of the cosmic rays are 
Galactic in origin. Any association with our Galaxy presumably disappears at EeV 
energy when the gyroradius of a proton in the Galactic magnetic field exceeds its 
size. The cosmic-ray spectrum exhibits a rich structure above an energy of a few 
PeV, the so-called “knee” in the spectrum, but where exactly the transition to 
extragalactic cosmic rays occurs is a matter of debate. 


2.1. Cosmic-Ray Accelerators 


The detailed blueprint for a cosmic-ray accelerator must meet two challenges: the 
highest-energy particles in the beam must reach energies beyond 10? TeV (108 TeV) 
for Galactic (extragalactic) sources and their luminosities must accommodate the 
observed flux. Both requirements represent severe constraints that have guided the- 
oretical speculations. Acceleration of protons (or nuclei) to TeV energy and above 
requires massive bulk flows of relativistic charged particles. Instead of accelerating 
protons, cosmic accelerators can boost really large masses to relativistic velocities. 
The radio emission reveals that the plasma in the jets of active galaxies flows with 
velocities of 0.99c. A fraction of a solar mass per year can be accelerated to rela- 
tivistic Lorentz factors of order 10 leading to luminosities of 10° erg/s close to the 
Eddington limit. In the collapse of very massive stars 10°! ~ 10° erg/s is released 
in a fireball that expands with velocities of 0.999 c. 

The blueprint of the accelerator can be copied from solar flares where particles 
are accelerated to GeV energy by shocks and, possibly, magnetic reconnection; see 
Fig. 3. Requiring that the gyroradius of the accelerated particle be contained within 
the accelerating B-field region, E/ZevB < R, leads to an upper limit on the energy 
of the particle, the Hillas?” formula 


E<ZevBR. (1) 


Reaching energies much above 10 GeV in solar flares is dimensionally impossible. 
In a solar flare, the extent R of the accelerating region and the magnitude of the 


*We will use energy units TeV, PeV and EeV, increasing by factors of 1000 from GeV energy. 
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Fig. 2. At the energies of interest here, the cosmic-ray spectrum follows a sequence of three power 
laws. The first two are separated by the “knee,” the second and third by the “ankle.” Cosmic rays 
beyond the ankle are a new population of particles produced in extragalactic sources. Note that 
the spectrum F(E)(= dN/dE) has been multiplied by a power E®-7 in order to visually enhance 
the structure in the spectrum. Figure from Particle Data Group.?! 


magnetic fields B are not large enough to accelerate particles of charge Ze to energies 
beyond GeV even if their velocity is taken to be the speed of light, c. Another way 
to view the dimensional argument is by estimating the energy of a particle from the 
Lorentz force 


R AB .,10°s 


d® dB 2 
10km) (Gral At )» 


E = —e— = erR? — ~ 1GeV( 


dt dt (2) 


where ® is the magnetic flux set up in the loops of gyrating particles in Fig. 3 that 
can reach several thousand Gauss in a time of order one day. The result follows 
from the loops have radii of more than 104km. In the spirit of dimensional analysis 
the two estimates above are the same by identifying the velocity v with dR/dt. 
While it is not a challenge to find astronomical sources with larger B and R, 
the other challenge is that the luminosity of the cosmic ray sources is large as well, 
and here a central idea for accommodating the high luminosities of the Galactic 
and extragalactic cosmic rays observed is that a fraction of the gravitational energy 
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released in a stellar collapse is converted into particle acceleration, presumably by 
shocks. 


Fig. 3. Opportunities exist near intense charged particle flows, seen as filaments in this X-ray 
picture of a solar flare, for solar particles to accelerate to GeV energy. 


Baade and Zwicky?’ suggested as early as 1934 that supernova remnants could 
be sources of the Galactic cosmic rays. It is assumed that, after the collapse, 
~ 10° erg of energy is transformed into particle acceleration by diffusive shocks 
associated with young (~1000 year old) supernova remnants expanding into the 
interstellar medium. Like a snowplow, the shock sweeps up the ~ 1 proton/cm? 
density of hydrogen in the Galactic plane. The accumulation of dense filaments 
of particles in the outer reaches of the shock, clearly visible as sources of intense 
X-ray emission, are the sites of high magnetic fields; see Fig. 4. It is theorized that 
particles crossing these structures multiple times can be accelerated to high energies 
following an approximate power-law spectrum dN/dE ~ E~?. The mechanism 
copies solar flares where filaments of high magnetic fields, visible in Fig. 3, are the 
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Fig. 4. This X-ray picture of the supernova remnant CasA reveals strong particle flows near its 
periphery. We believe they are the site for accelerating Galactic cosmic rays to energies reaching 
the “knee” in the spectrum. 


sites for accelerating nuclear particles to tens of GeV. The higher energies reached 
in supernova remnants are the consequence of particle flows of much larger intensity 
powered by the gravitational energy released in the stellar collapse. 

This idea has been widely accepted despite the fact that to date no source has 
been conclusively identified, neither by cosmic rays nor by accompanying gamma 
rays and neutrinos produced when the cosmic rays interact with Galactic hydro- 
gen. Galactic cosmic rays reach energies of at least several PeV, the “knee” in the 
spectrum; therefore, their interactions should generate gamma rays and neutrinos 
from the decay of secondary pions reaching hundreds of TeV. Such sources, referred 
to as PeVatrons, have not been found; see, however, reference 24. Nevertheless, 
Zwicky’s suggestion has become the stuff of textbooks, and the reason is energetics: 
three Galactic supernova explosions per century converting a reasonable fraction of 
a solar mass into particle acceleration can accommodate the steady flux of cosmic 
rays in the Galaxy. It is interesting to note that Zwicky originally assumed that 
the sources were extragalactic since the most recent supernova in the Milky Way 
was in 1572. After diffusion in the interstellar medium was understood, supernova 
explosions in the Milky Way became the source of choice for the origin of Galactic 
cosmic rays,”° although after more than 50 years the issue is still debated.7° 

Energetics also guides speculations on the origin of extragalactic cosmic rays. 
By integrating the cosmic-ray spectrum above the ankle at ~ 4 EeV, it is possible to 
estimate?’ the energy density in extragalactic cosmic rays as ~ 3 x 10719 erg cm™3. 
This value is rather uncertain because of our ignorance of the energy where the 
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Fig. 5. Colliding shocks in the simulation of a gamma-ray burst (GRB) fireball may 
accelerate cosmic rays to the highest energies observed. The filaments in the parti- 
cle flow are directed along the rotation axis of the black hole. Animated view at 
http: //www.nasa.gov/centers/goddard/news/topstory /2003/0618rosettaburst.html. 


transition from Galactic to extragalactic sources occurs. The power required for 
a population of sources to generate this energy density over the Hubble time of 
101° years is 2 x 10°” erg s~! per Mpc?. Long-duration gamma-ray bursts have 
been associated with the collapse of massive stars to black holes, and not to neutron 
stars, as is the case in a collapse powering a supernova remnant. A gamma-ray- 
burst fireball converts a fraction of a solar mass into the acceleration of electrons, 
seen as synchrotron photons. The observed energy in extragalactic cosmic rays can 
be accommodated with the reasonable assumption that shocks in the expanding 
gamma-ray burst (GRB) fireball convert roughly equal energy into the acceleration 
of electrons and cosmic rays;*° see Fig. 5. It so happens that 2 x 10°! erg per GRB 
will yield the observed energy density in cosmic rays after 10!° years, given that 
their rate is on the order of 300 per Gpc? per year. Hundreds of bursts per year over 
a Hubble time produce the observed cosmic-ray density, just as three supernovae 
per century accommodate the steady flux in the Galaxy. 

Problem solved? Not really: it turns out that the same result can be achieved 
assuming that active galactic nuclei convert, on average, 2 x 1044 erg s~! each into 
particle acceleration.'° This is an amount that matches their output in electromag- 
netic radiation. 

In contrast with our own Galaxy where the black hole is mostly dormant, in an 
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Fig. 6. Anticipated astrophysical neutrino fluxes compared with measured and calculated fluxes 
of atmospheric neutrinos. Measurements of v, from Super-K,?° Frejus,?? AMANDA,°!°? and 
IceCube??)* are shown along with the electron-neutrino spectrum at high energy from reference?" 
(green open triangles). Calculations of conventional ve (red line) and v, (blue line) from Honda 
et al.,°° ve (red dotted line) from Bartol,?” and charm-induced neutrinos (magenta band)*® are 


also shown. 


active galaxy the supermassive black hole is absorbing the matter in its host galaxy 
at a very high rate. An active galactic nucleus (AGN) hosts a rotating supermassive 
black hole. Fast spinning matter falling onto it swirls around the black hole in an 
accretion disk, like the water approaching the drain of your bath tub. When the 
accretion disk comes in contact with the rotating black hole space-time drags on 
the magnetic field winding it into a tight cone around the rotation axis into a jet 
of particles; see Fig. 7. Not just particles but huge “blobs” of plasma from the 
accretion disk are flung out along these field lines. It is not clear whether it is the 
rotation energy of the black hole or the magnetic energy in the rotating plasma that 
powers the accelerator. When this jet runs into a target material, for instance the 
ubiquitous 10 eV ultraviolet photons in some galaxies, neutrinos can be produced. 

Active galaxies are actually complex systems with many possible sites for accel- 
erating cosmic rays and for targets to produce neutrinos. Acceleration of particles 
may occur at the spectacular termination shocks of the jets in intergalactic space’? 
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at distances of hundreds of Mpc from the center where there would be little target 
material available and few neutrinos produced. In contrast, production of neutri- 
nos near the black hole,*° or in collisions with interstellar matter of the accelerated 
particles diffusing in the magnetic field of the galaxy hosting the black hole,*! could 
yield fluxes at the level observed. We will work through these examples further on. 


— Jet of high- 
particles 


Accretion 
disk 


Fig. 7. The accretion disk meets the spinning black holes that winds up the disk’s magnetic field 
lines. Credit: http://www.physics.ucc.ie/gabuzda/Nadia.html. 
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2.2. Neutrinos and Gamma Rays Associated with Cosmic Rays 


Neutrinos must be produced at some level in association with the cosmic-ray beam. 
Cosmic rays accelerated in regions of high magnetic fields near black holes or neutron 
stars inevitably interact with the matter or radiation surrounding them. Thus, 
cosmic-ray accelerators are also part of a “beam dump”, like the ones producing 
neutrino beams in accelerator laboratories: the beam is dumped is a target where 
it produces pions and kaons that decay into neutrinos. Dense targets absorb all 
secondary particles except for the neutrinos. This is typically not the case for a 
cosmic beam dump where neutrinos are expected to be accompanied by other stable 
particles: protons, neutrons and photons. For example, cosmic rays accelerated in 
supernova shocks interact with gas in the Galactic disk, producing equal numbers 
of pions of all three charges that decay into pionic photons and neutrinos. A larger 
source of secondaries is likely to be produced by the interaction of accelerated 
particles with the gas near the sources, for example cosmic rays interacting with 
high-density molecular clouds that are ubiquitous in the star-forming regions where 
supernovae are more likely to explode. For extragalactic sources, the neutrino- 
producing target may be electromagnetic, for instance photons radiated by the 
accretion disk of an AGN, or synchrotron photons that coexist with protons in the 
expanding fireball producing a GRB. 

How many neutrinos and, inevitably, gamma rays are produced in association 
with the cosmic-ray beam? A Galactic supernova shock is an example of a hadronic 
beam dump. Cosmic rays interact with the hydrogen in the Galactic disk, producing 
equal numbers of pions of all three charges in hadronic collisions p + p > ny [7° + 
mt +r] +X; n, is the pion multiplicity. In the case of a photon target, neutral 
and charged pion secondaries are produced by the photoproduction processes 

pty or Ator4+p and ptyo At >rt tn. (3) 
Only the neutrinos and neutrons will escape the source. While secondary protons 
may remain trapped and loose energy in the high magnetic fields of the accelerator, 
neutrons will decay and the decay products escape the dump with high energy. 
The energy escaping the source is therefore distributed among cosmic rays, gamma 
rays and neutrinos, particles produced by the decay of neutrons, neutral pions and 
charged pions, respectively. Photoproduction produces charged and neutral pions 
according to Eq. 3, with probabilities of 2/3 and 1/3, respectively. Subsequently, 
the pions decay into gamma rays and neutrinos that carry, on average, 1/2 and 1/4 
of the energy of the parent pion. It is a good approximation to assume that, on 
average, the four leptons in the decay tt > vy + ut > vp + (et + ve + Du) equally 
share the charged pion’s energy. The energy of the pionic leptons relative to the 


proton is: 
E, 1 1 
v = t © 4 
Se T (4) 
and m j i 
he = 
t= ee T To (5) 
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Er 
£r = (—) & 0.2 6 
n= (FR) (6) 
is the average energy transferred from the proton to the pion that produces the 
neutrino.*? 


Interestingly, these relations are approximately valid whether the pions are pro- 
duced in pp or py interactions. Pion production is usually described in terms of 
their average multiplicity ny and a pion’s average energy E, in the final state. Af- 
ter interaction, the initial state energy is distributed between the proton and the 
production of pions. The fraction of energy going into the production of pions is re- 
ferred to as the inelasticity « which is 0.5(0.2) for pp(py) interactions. The fraction 
of energy going into a single pion x, is 

K 
Te = A ~ 0.2. (7) 
It turns out that 2, ~ 0.2 yields the pion energy for both pp and py interactions, 
despite the very different particle processes. For photoproduction, the charged pion 
takes 0.2 of the initial energy and the pion multiplicity is 1. 

While both gamma-ray and neutrino fluxes can be calculated knowing the lu- 
minosity of the accelerated protons and the density of the target material, their 
relative flux is independent of the details of the production mechanism. Their pro- 
duction rates dN/dEdt of neutrinos and gamma rays are related by known particle 
physics. The above discussion can be summarized as: 

1 dN, Kr dN, 

pL Ba RE) ~ Soap (8) 
Here, N and E denote the number and energy of neutrinos and gamma rays and va 
stands for the neutrino flavor. Note that this relation is solid and depends only on 
the charged-to-neutral secondary pion ratio, with K, = 1(2) for py(pp) neutrino- 
producing interactions. In deriving the relative number of neutrinos and gamma 
rays, one must be aware of the fact that the neutrino flux represents the sum of the 
neutrinos and antineutrinos, which cannot be separated by current experiments: in 
short, a 7° produces two y rays for every charged pion producing a v,, +, pair. A 
more detailed discussion of this relation will follow in the context of multimessenger 
astronomy; see Section 7. 

The production rate of gamma rays at their origin described by Eq. 8 is not 
necessarily the emission rate observed. For instance, in cosmic accelerators that 
efficiently produce neutrinos via py interactions, the target photon field can also 
efficiently reduce the energy of the pionic gamma rays produced via pair production. 
Gamma rays with energies above the threshold for pair production will lose energy 
in the source. Their maximum energy, in the comoving frame, is determined by the 
energy of the target photons: 


By Eln < me. (9) 
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This is a calorimetric process that will, however, conserve the total energy of 
hadronic gamma rays. The production of photons in association with cosmic neutri- 
nos is inevitable but unlike neutrinos, photons may reach Earth with reduced energy 
after losing energy in the target and after propagation in the universal microwave 
and infrared photon backgrounds. However, their cascaded energy must appear in 
some electromagnetic wave band because of energy conservation. Furthermore, one 
must be aware of the fact that inverse-Compton scattering and synchrotron emis- 
sion by accelerated electrons in magnetic fields in the source have the potential to 
produce gamma rays; not every high-energy gamma ray is pionic. 

A couple of decades of modeling potential neutrino sources yielded generic pre- 
dictions shown in Fig. 6, estimates of astrophysical neutrino fluxes are compared 
with measurements of atmospheric neutrinos. While the models varied, the result 
were typically dictated by the relation between the gamma ray and neutrino fluxes 
discussed above. The shaded band indicates the level of model-dependent expecta- 
tions for high-energy neutrinos of astrophysical origin. The estimates that we will 
discuss in more detail further on predicted a neutrino flux at a level of 


E? dN, /dE, ~ 1078 GeV cm~?s~'sr7* (10) 


per flavor. The figure illustrates the rationale for building a kilometer-scale detector 
because it yields about 100 neutrino events per year in a cubic kilometer detector. 
We now know that this is indeed the magnitude of the cosmic component of the 
neutrino spectrum above 100TeV revealed by IceCube’s data although we have 
found no evidence for the specific sources shown in the picture! The ongoing search 
by IceCube for neutrinos in coincidence with and in the direction of GRB alerts 
issued by astronomical telescopes has limited the GRB neutrino flux to less that 
1% of the diffuse cosmic neutrino flux actually observed by the experiment.*?“4 
However, this may not conclusively rule out GRBs as a source of cosmic rays; the 
events that produce the spectacular photon displays catalogued by astronomers as 
GRBs may not be the stellar collapses that are sources of high-energy neutrinos. 
We will return to this point further on when we discuss acceleration of cosmic rays 
in GRB fireballs. 

The failure of IceCube to observe neutrinos from GRBs has lately promoted 
AGNs as the best-bet source of the cosmic neutrinos observed. Here again 
Gaisser*?*° has emphasized the relation between IceCube limits and the electro- 
magnetic energy of the sources; see Fig. 8.4 In this context, we introduce Fig. 847 
showing IceCube upper limits*® on the neutrino flux from nearby AGNs as a func- 
tion of their distance. The sources at red shifts between 0.03 and 0.2 are Northern 
Hemisphere blazars for which distances and intensities are listed in TeVCat*® and 
for which IceCube also has upper limits. In several cases, the muon-neutrino limits 
have reached the level of the TeV photon flux. One can sum the sources shown in 
the figure into a diffuse flux. The result, after accounting for the distances and lumi- 
nosities, is 3 x 1079 GeV cm~? s7! sr™!, or approximately 1078 GeV cm~? s7! sr7! 
for all neutrino flavors. This is at the level of the generic astrophysical neutrino flux 
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of Eq. 10. At this intensity, neutrinos from theorized cosmic-ray accelerators will 
cross the steeply falling atmospheric neutrino flux above an energy of ~ 300 TeV; 
see Fig. 6. The level of events observed in a cubic-kilometer neutrino detector is 
10 ~ 100 v,,-induced events per year. Such estimates reinforce the logic for building 
a cubic kilometer neutrino detector.°? 
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Fig. 8. Limits on the neutrino flux from selected active galaxies derived from IceCube data taken 
during construction, when the instrument was operating with 40 and 59 strings of the total 86 
instrumented strings of DOMs.*® These are compared with the TeV photon flux for nearby AGNs. 
Note that energy units are in ergs, not TeV. 


3. IceCube 


3.1. Detecting Very High Energy Neutrinos 


Cosmic rays have been studied for more than a century. They reach energies in 
excess of 10° TeV, populating the extreme universe that is opaque to photons be- 
cause they interact with the background radiation fields, mostly microwave photons, 
before reaching Earth. We don’t yet know where or how cosmic rays are acceler- 
ated to these extreme energies, and with the recent observation of a the rotating 
supermassive black hole TXS 0506+056 in coincidence with the direction and time 
of a very high energy muon neutrino, neutrino astronomy might have taken a first 
step in solving this puzzle.°!°? The rationale is however simple: near neutron stars 


page 15 


February 3, 2022 1:35 ws-rv961x669 Neutrino Physics and Astrophysics main page 16 


16 IceCube Neutrinos 


cosmic 


ray \ 
atmospheric 
muon 


cosmic 
neutrino IceCube 


TYY 
on 


$ upone A l l l 
Cherenkov light detection 


in optical modules 
cosmic a 


ray 


t down-go atmospheric 


neutrino 


Fig. 9. The concept of a neutrino telescope from the point of view of IceCube located at the 
South Pole. Muon neutrinos interact with the nuclei of atoms in a transparent medium like water 
or ice and produce a secondary muon track that is detected by the wake of Cherenkov photons it 
leaves inside the detector. The background of high-energy muons (solid blue arrows) produced in 
the atmosphere can be reduced by placing the detector underground. The surviving fraction of 
muons is further reduced by looking for upgoing muon tracks that originate from muon neutrinos 
that have traversed the Earth (dashed blue arrows) interacting close to the detector. This still 
leaves the contribution of muons generated by atmospheric muon neutrino interactions. This 
contribution can be separated from the diffuse cosmic neutrino emission by an analysis of the 
combined neutrino spectrum. 


and black holes, gravitational energy released in the accretion of matter or binary 
mergers can power the acceleration of protons or heavier nuclei that subsequently 
interact with gas (“pp”) or ambient radiation (“py”). Neutrinos are produced by 
cosmic-ray interactions at various epochs: in their sources during their accelera- 
tion, in the source environment after their release, and while propagating through 
universal radiation backgrounds from the source to Earth. 

Because of their weak interactions, high-energy neutrinos will reach our detectors 
without deflection or absorption. They essentially act like photons; their small mass 
is negligible relative to the TeV to EeV energies targeted by neutrino telescopes. 
They do however oscillate over cosmic distances. For instance, for an initial neutrino 
flavor ratio of Ve : Va : Vr ~ 1 : 2 : 0 from the decay of pions and muons, the 
oscillation-averaged composition arriving at the detector is approximately an equal 
mix of electron, muon, and tau neutrino flavors, Ve : Vp : Vr S 1:1:1. 

High-energy neutrinos interact predominantly with matter via deep inelastic 
scattering off nucleons: the neutrino scatters off quarks in the target nucleus by 
the exchange of a Z or W weak boson, referred to as neutral current (NC) and 
charged current (CC) interactions, respectively. Whereas the NC interaction leaves 
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the neutrino state intact, in a CC interaction a charged lepton is produced that 
shares the initial neutrino flavor. The average relative energy fraction transferred 
from the neutrino to the lepton is at the level of 80% at high energies. The inelastic 
CC cross section on protons is at the level of 10~°3 cm? at a neutrino energy of 
10? TeV and grows with neutrino energy as Otot X H930 945 The struck nucleus 
does not remain intact and its high-energy fragments typically initiate hadronic 
showers in the target medium. 

Immense particle detectors are required to collect cosmic neutrinos in statisti- 
cally significant numbers. Already by the 1970s, it had been understood? that a 
kilometer-scale detector was needed to observe the cosmogenic neutrinos produced 
in the interactions of cosmic rays with background microwave photons.°°°’ A va- 
riety of methods are used to detect the high-energy secondary particles created in 
CC and NC neutrino interactions. One particularly effective method observes the 
radiation of optical Cherenkov light radiated by secondary charged particles pro- 
duced in CC and NC interactions that travel faster than the speed of light in the 
medium. 

The detection concept is that of a conventional Cherenkov detector, a transpar- 
ent medium is instrumented with photomultipliers that transform the Cherenkov 
light into electrical signals by the photoelectric effect; see Figs. 9 and 10. IceCube 
consists of 80 strings, each instrumented with 60 10-inch photomultipliers spaced 
by 17m over a total length of 1 kilometer. The deepest module is located at a 
depth of 2.450km so that the instrument is shielded from the large background of 
cosmic rays at the surface by approximately 1.5km of ice. Strings are arranged at 
apexes of equilateral triangles that are 125m on a side. The instrumented detector 
volume is a cubic kilometer of dark, highly transparent and sterile Antarctic ice. 
The radioactive background in the detector is dominated by the instrumentation 
deployed into this natural ice. 

Each optical sensor consists of a glass sphere containing the photomultiplier and 
the electronics board that captures and digitizes the signals locally using an on- 
board computer. The digitized signals are given a global time stamp with residuals 
accurate to less than 3ns and are subsequently transmitted to the surface. Proces- 
sors at the surface continuously collect the time-stamped signals from the optical 
modules, each of which functions independently. The digital messages are sent to 
a string processor and a global event trigger. They are subsequently sorted into 
the Cherenkov patterns emitted by secondary muon tracks, or particle showers for 
electron and tau neutrinos, that reveal the flavor, energy and direction of the in- 
cident neutrino.°® There are two principle classes of Cherenkov events that must 
be separated by the detector, “tracks” and “cascades”. The two basic topologies 
are illustrated in Fig. 11: tracks initiated by v, neutrinos and cascades from ve, Vr 
neutrinos as well as the neutral current interactions from all flavors. On the scale 
of IceCube, PeV cascades, with a length of less than 10 m, are therefore essentially 
point sources of Cherenkov light in a detector of kilometer size. The term “tracks” 
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Fig. 10. Sketch of the IceCube observatory (left) and the digital optical module (right). 
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Fig. 11. Contrasting Cherenkov light patterns produced by muons (left) and by showers initiated 
by electron and tau neutrinos (right) and by neutral current interactions. The patterns are re- 
ferred to as tracks and cascades (or showers). Cascades are produced by the radiation of particle 
showers, whose dimensions are less than 10m, i.e., essentially a point source of light relative to 
the dimensions of the detector. 


refers to the Cherenkov emission of long-lived muons produced in CC interactions 
of muon neutrinos interacting inside or in the vicinity of the detector. Energetic 
electrons and taus produced in CC interactions of electron and tau neutrino in- 
teractions, respectively, will not produce elongated tracks: electrons will initiate 
electromagnetic showers in the ice and, because of the relatively short lifetime, taus 
will decay and also produce a shower in the ice. Because of the large background of 
muons produced by cosmic ray interactions in the atmosphere, the observation of 
muon neutrinos is typically limited to upgoing muon tracks that are produced in in- 
teractions inside or close to the detector by neutrinos that have passed through the 
Earth as illustrated in Fig. 9. The remaining background consists of atmospheric 
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neutrinos, which are indistinguishable from cosmic neutrinos on an event-by-event 
basis. However, the steeply falling spectrum (x Æ") of atmospheric neutrinos 
allows identifying diffuse astrophysical neutrino by a spectral analysis as we will 
highlight in the following sections. Above ~ 300TeV atmospheric neutrinos are 
relatively rare, even in a cubic kilometer detector, and every neutrino exceeding 
this energy is likely to be of cosmic origin. The atmospheric background is fur- 
ther reduced when looking in the specific direction of a point source, in particular 
transient neutrino sources like GRB. 

The hadronic particle shower generated by the target struck by a neutrino also 
radiates Cherenkov photons. Because of the large multiplicity of secondary particles 
at these energies and the repeated scattering of the Cherenkov photons in the ice, the 
light pattern is essential spherical tens of meters from its point of origin. The light 
patterns produced by the particle showers initiated by the electron or tau produced 
in CC interactions of electron or tau neutrinos, respectively, will be superimposed 
on the hadronic cascade. The two are not typically not separated. The direction 
of the initial neutrino can only be reconstructed from the Cherenkov emission of 
secondary particles produced close to the neutrino interaction point, and the angular 
resolution is inferior to that for track events. 

In contrast, the energy resolution of the neutrino is superior for cascades than 
for tracks. For both, the observable energy of the secondaries can be estimated 
from the total number of Cherenkov photons after accounting for kinematic effects 
and detection efficiencies. The Cherenkov light observed in cascades is proportional 
to the energy transferred to the cascade and is often fully contained in the instru- 
mented volume. It is actually sufficient that the cascade is partially contained, as 
long as one can reconstruct its actuals size. In contrast, muons produced by CC 
muon neutrino interactions lose energy gradually by ionization, bremsstrahlung, 
pair production, and photo-nuclear interactions while passing through the detector. 
The secondary charged particles produced in each energy-loss interaction radiate 
Cherenkov photons that allow for a measurement of the total energy lost by the 
muon in the detector. The energy deposited by the muon in the detector represents 
a lower limit on the initial neutrino energy. 

Catastrophically losing energy by the process mentioned above, muon tracks 
range out, over kilometers at TeV energy to tens of kilometers at EeV energy. 
Because the energy of the muon thus degrades along its track, the energy of the 
secondary showers decreases, which reduces the distance from the track over which 
the associated Cherenkov light can trigger a PMT. The geometry of the light pool 
surrounding the muon track is therefore a kilometer-long cone with a gradually 
decreasing radius. On average, in its first kilometer, a high-energy muon loses 
energy in a couple of showers with one-tenth of the muon’s initial energy. So the 
initial radius of the cone is the radius of a shower with 10% of the muon energy. At 
lower energies of hundreds of GeV and less, the muon becomes minimum-ionizing. 

Because of the stochastic nature of the muon’s energy loss, the relationship 
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between the observed energy loss inside the detector and the muon energy varies 
from muon to muon. Additionally, only the muon energy lost in the detector can 
be determined; we do not know the energy lost before entering the instrumented 
volume, nor how much energy it carries out upon exiting. An unfolding process 
is required to derive the neutrino energy from the energy lost inside the detector; 
fortunately, it is described by well-understood Standard Model physics. One derives 
a probability distribution for the energy of the initial neutrino that determines its 
most probable value. The neutrino energy may only be determined within a factor 
of 2 or thereabout, depending on the energy, but the uncertainties drop out when 
measuring a neutrino spectrum involving multiple events. In contrast, for ve and 
vz, the detector is a total energy calorimeter capturing all or most of the Cherenkov 
light produced, and the determination of their energy is superior. 

The different topologies each have advantages and disadvantages. For v, 
charged-current interactions, the long lever arm of muon tracks allows for a mea- 
surement of the muon direction with an angular resolution of better than 0.4°. 
Superior angular resolution can be reached for selected high-energy events. At the 
highest energies the neutrino is aligned with the muon within the angular resolution 
and the sensitivity to point sources thus maximized. The disadvantages are a large 
background of atmospheric neutrinos below 100 TeV, and of cosmic-ray muons at 
all energies, and the indirect determination of the neutrino energy that must be 
inferred from sampling the energy loss of the muon when it transits the detector. 

Observation of ve and v, flavors represents significant advantages. They are 
detected from both Northern and Southern Hemispheres. (This is also true for v, 
with energy in excess of several hundred TeV, where the background from the steeply 
falling atmospheric spectrum becomes negligible.) At TeV energies and above, the 
background of atmospheric Ve is lower by over an order of magnitude because long- 
lived pions, the source of atmospheric ve, no longer decay, and relatively rare K- 
decays become the dominant source of background ve. Atmospheric v+, produced 
by oscillations, are rare above an energy of 100 GeV. High-energy v, are therefore 
of cosmic origin; one such event with an energy of 100 TeV has been identified and 
represents an independent discovery of cosmic neutrinos.°®? Furthermore, one can 
establish the cosmic origin of a single single cascade event by demonstrating that 
the energy cannot be reached by muons and neutrinos of atmospheric origin. 

Finally, v, are not absorbed by the Earth:°° v, interacting in the Earth produce 
a secondary v, of lower energy, either directly in a neutral current interaction or 
via the decay of a secondary tau lepton produced in a charged-current interaction. 
High-energy v, will thus cascade down to energies of hundred of TeV where the 
Earth becomes transparent. In other words, they are detected with a reduced 
energy, but not absorbed. By this mechanism GZK v, with EeV energies produce 
a signal at PeV energy.°! 

Although cascades are nearly point-like sources of Cherenkov light and, in prac- 
tice, spatially isotropic, the pattern of arrival times of the photons at individual 
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optical modules reveals the direction of the secondary lepton. While a fraction of 
cascade events may be reconstructed to within a degree,°? the precision is inferior to 
that reached for v,, events, typically 8° using the present techniques. A campaign is 
underway to better characterize the optical properties of the ice and the calibration 
of the detector in order to improve this resolution. 

At energies above about 100 PeV, electromagnetic showers begin to elongate 
because of the Landau-Pomeranchuk-Migdal effect. An extended length scale, 
associated with the abundant radiation of soft photons, results in the interaction of 
the secondary shower particles with two atoms. Negative interference of this process 
relative to interactions with a single atom results in a reduction of the energy loss. 


3.2. The Neutrino Detector’s Telescope Area 


Cosmic neutrinos must be separated from the large backgrounds of atmospheric 
neutrinos and atmospheric cosmic-ray muons. Two principal methods have been 
developed: isolating neutrinos that interact inside the instrumented volume (“start- 
ing events”), and specializing to events where a muon enters the detector from be- 
low, created by a neutrino that has traversed the Earth (“throughgoing events”), 
thus pointing back to its origin. In the latter case, the Earth is used as a filter for 
cosmic-ray muons. 

For starting events, neutrinos are detected provided they interact within the 
detector volume, i.e., within the instrumented volume of one cubic kilometer. That 
probability is 


P(E.) =1- exp[—1/A,(E_)] ~ l/r. (Ev) , (11) 


where (0) is the path length traversed by a neutrino with zenith angle 0 within 
the detector volume and A, (Ep) = [pice Na ovn(E_)|~! is the mean free path in 
ice for a neutrino of energy E,. Here, pice = 0.9 gem~? is the density of the ice, 
N4 = 6.022 x 10? is Avogadro’s number, and o,,y(E_) is the neutrino-nucleon cross 
section. The path length is determined by the detector’s geometry and is typically 
much shorter than the neutrino mean free path Ap. 

A neutrino flux dN/dE, (with typical units neutrinos per GeV or erg per cm? 
per second) crossing a detector with energy threshold Æt! and cross sectional area 
A(E,) facing the incident neutrino beam will produce 


New =T | ACB.) P(B,) aF, (12) 
Eth v 

events after a time T. The “effective” detector area A(E,) is a function of neutrino 
energy and of the zenith angle 0. It isn’t strictly equal to the geometric cross section 
of the instrumented volume facing the incoming neutrino, because even neutrinos 
interacting outside the instrumented volume may produce enough light inside the 
detector to be detected. In practice, A(£,) is determined as a function of the 
incident neutrino direction and zenith angle by a full-detector simulation, including 
the trigger and the cuts that are used to isolate “starting” events. 
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In contrast, throughgoing muon neutrinos will be detected provided the sec- 
ondary muon reaches the detector with sufficient energy to trigger it. Because 
the muon travels kilometers at TeV energy and tens of kilometers at PeV energy, 
neutrinos are detected outside the instrumented volume with a probability 


P = Au /àv (13) 
obtained by the substitution 
l> An, (14) 


in Eq. 11. Here, A, is the range of the muon determined by its energy losses. Values 
for the neutrino nucleon cross section and the range of the muon can be found in 
reference 64. 

At energies above tens of TeV, one has to account for the fact that the neutrinos 
may be absorbed in the Earth before reaching the detector. The reduced flux of 
v,,-induced muons reaching the detector is given by:*° 


oA 8) = | _ P(E,, Ep”) exp[—otot(E.) Na X (0) (By, 0)dE,. (15) 
Amin 
The additional exponential factor accounts for the absorption of neutrinos along a 
chord through the Earth of length X(6) at zenith angle 6. 
For back-of-the-envelope calculations, the P-function describing the probability 
that a neutrino is detected can be approximated by 


P~1.3 x 10-°E?? for E = 1073-1 TeV, (16) 
x 1.3 x 10-8E°? for E = 1-108 TeV. (17) 


At EeV energy, the increase is reduced to only E°4. The parametrization describes 
how neutrinos of higher energy are more likely to be detected because of the increase 
with energy of both the cross section and muon range. But at neutrino energies of 
tens of TeV and above, this gain is partially mitigated by absorption in the Earth. 
The highest energy neutrinos reach the detector from zenith angles near the horizon. 

Tau neutrinos interacting outside the detector can be observed provided the 
secondary tau lepton reaches the instrumented volume within its lifetime. In Eq. 11, 
l is replaced by 


E 
l>Trer= —cr, (18) 
m 


where m, T, and E are the mass, lifetime, and energy of the tau, respectively. The 
tau’s decay length A, = [cr ~ 50 m x (E,/PeV) grows linearly with energy and 
exceeds the range of the muon near 1 EeV. At yet higher energies, the tau eventually 
ranges out by catastrophic interactions, just like the muon, despite the reduction of 
the leading energy-loss cross sections by a factor of (m,/m-,)?. 

At sub-PeV energies, tracks and showers produced by tau neutrinos are difficult 
to distinguish from those initiated by muon and electron neutrinos, respectively. 
Only at PeV energies it is possible to detect both the initial neutrino interaction 
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and the subsequent tau decay that are separated by tens of meters. Additionally, 
both must be contained within the detector volume; for a cubic-kilometer detector, 
this can be realized for neutrinos with energies from a few hundreds of TeV to a 
few tens of PeV.®° 

For an in-depth discussion of neutrino detection, energy measurement, and fla- 
vor separation, and for detailed references, see the IceCube Preliminary Design 


Document!® and reference 2. 


3.3. Atmospheric Neutrinos: Calibration and Background 


The 3-kHz trigger rate of the IceCube detector is dominated by atmospheric muons 
from the decay of pions and kaons produced in the atmosphere above the detector. 
Their distribution peaks near the zenith and decreases with increasing angle because 
the muon energy required to reach the deep detector increases. Most atmospheric 
muons are identified as tracks entering the detector from above and are rejected 
because the Earth shields the detector from atmospheric muons in the Northern 
Hemisphere. Even after the removal of cosmic ray muons, the neutrinos from decay 
of mesons produced by cosmic-ray interactions in the atmosphere are a residual 
background in the search for neutrinos of extraterrestrial origin. Because of the 
large ratio of atmospheric muons to neutrinos misreconstructed atmospheric muons 
remain an important source of background for most searches. 

Measurement of the relatively well-established spectrum of atmospheric neutri- 
nos is an important benchmark and a useful calibration tool for a neutrino telescope. 
IceCube detects an atmospheric neutrino every few minutes, more than one hundred 
thousand per year. The spectrum of atmospheric v,, has been measured by unfold- 
ing the measured rate and energy deposition of neutrino-induced muons entering 
the detector from below the horizon,** as shown in Fig. 6. More challenging is the 
measurement of the flux of atmospheric electron neutrinos. This has been achieved 
by making use of DeepCore, the more densely instrumented subarray in the deep 
center of IceCube, to identify contained shower events. The measured spectrum 
of vy, is used to calculate the contribution of neutral current interactions to the 
observed rate of showers. Subtracting the neutral current contribution leads to the 
measurement of the spectrum of atmospheric electron neutrinos from 100 GeV to 
10 TeV;°° see Fig. 6. 

In general, atmospheric neutrinos are indistinguishable from astrophysical neu- 
trinos. An important exception is when muon neutrinos reaching the detector from 
above can be tagged as atmospheric by detecting the muon produced in the same 
decay as the neutrino. The neutrino energy must be sufficiently high and the zenith 
angle sufficiently small that this muon to reach the detector.°° Monte Carlo simu- 
lation are performed to evaluate the rejection rate or, equivalently, the atmospheric 
neutrino passing rate. Also high-energy muons other than the muon associated 
directly with the neutrino and produced in the same cosmic-ray shower as the neu- 
trino are included in the veto. In this way, the method can be extended to electron 
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neutrinos. In practice, the passing rate is significantly reduced for EF, > 100 TeV 
and zenith angles 0 < 70° where the ice overburden is not too large. 

The spectrum of atmospheric neutrinos becomes one power steeper than the 
spectrum of primary nucleons at high energy because the competition between 
interaction in the atmosphere and the decay of pions and kaons increasingly sup- 
presses their decay. For the kaon channel, dominant at high energies, the charac- 
teristic energy for the steepening is E, ~ 1 TeV/cos@. A further steepening occurs 
above 100 TeV as a consequence of the knee in the primary cosmic ray spectrum. 
In contrast, astrophysical neutrinos should reflect the cosmic-ray spectrum of the 
cosmic accelerator expected to be significantly harder spectrum than atmospheric 
neutrinos. Establishing an astrophysical signal above the steep atmospheric back- 
ground requires an understanding of the atmospheric neutrino spectrum at energies 
of 100 TeV and above. 

Although there is some uncertainty associated with the composition through 
the knee region,°’ the major uncertainty in the spectrum of atmospheric neutrinos 
at high energy is the level of charm production. The short-lived charmed hadrons 
preferentially decay, up to a characteristic energy of 10° GeV, producing “prompt” 
muons and neutrinos with the same spectrum as their parent cosmic rays. This 


68,69 allow a factor 
38 


prompt flux of leptons has yet to be measured. Existing limits 
of two or three around the level predicted by model calculations. 
able assumptions, the charm contribution is expected to dominate the conventional 


spectrum above ~ 10 TeV for ve, above ~ 100 TeV for v,,, and above ~ 1 PeV for 
70 


For reason- 


muons. 
The expected hardening in the spectrum of atmospheric neutrinos due to prompt 
neutrinos, is partially degenerate with a hard astrophysical component. However, 
the spectrum of astrophysical neutrinos should reflect the spectrum of cosmic rays 
at their sources, which is expected to be harder than the spectrum of cosmic rays 
at Earth. It should eventually be possible with IceCube to measure the charm 
contribution by requiring a consistent interpretation of neutrino flavors and cosmic- 
ray muons for which there is no astrophysical component. An additional signature 
of atmospheric charm is the absence of seasonal variations for this component.’! 
As we will discuss further on, with a good understanding of the energy and zenith 
angle dependence of the atmospheric neutrino spectrum, supplemented by the veto 
technique and the shielding of the muons by the Earth, IceCube has been successful 
in separating a diffuse flux of cosmic neutrinos from the atmospheric backgrounds. 
The next step is to find the origin of this flux by identifying individual sources. 
The strategy of searching for neutrino sources is to look for spatial clustering 
in the arrival direction of neutrinos to find any excess over the expected isotropic 
distribution of background. The technique used by IceCube to search for point 
sources is described in reference 72. In this method, an unbinned maximum likeli- 
hood is constructed to search for spatial clustering of the events. Significances are 
estimated by repeating each hypothesis test on data sets that are randomized in 
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right ascension and dominated by background. This provides robust p-values that 
are largely independent of detector systematic uncertainties. 

The unbinned maximum likelihood ratio method used to look for a localized, 
statistically significant excess of events above the background allows full use of 
spatial and spectral information from the data. The data are hypothesized to be a 
mixture of events from signal and background. 

For an event with reconstructed direction 7; = (a;,6;), the probability of orig- 
inating from the source at x, is modeled as a circular two-dimensional Gaussian. 
The signal probability distribution function (PDF) S; incorporates directional in- 
formation for each individual event and its angular uncertainty, 0;, and the angular 
difference between the reconstructed direction of the event and the source: 


Si = Si(|Xi — Psl, 0%) E:(0i, Fi, 7), (19) 
where the spatial distribution is modeled as a two-dimensional Gaussian 
rae 1 |Z; — Zl? 
Sill = Bla) = zg OP ( = ) (20) 


The background PDF, 6;, contains similar terms that describe the angular and 
energy distributions of background events. The likelihood function for a point source 
is defined as 


events 
Ts N= Ns 
L(ns, £s,7) = Il (sila — zs|, oi, Ei, y) + TB, E) (21) 


i 
The likelihood ratio test statistic (TS) is used to perform statistical tests. The 
results of these tests can be clearly defined in the context of testing between two 
hypotheses: the null hypothesis Ho and the alternative hypothesis Hı that signal 
events exceed the background. Hp represents the case of background only (ns = 0). 
After maximizing and determining the best fit number of signal events ñs and 
their spectral index ĝĵ, the test statistic (TS) is defined as the log likelihood ratio 
between the null and signal hypothesis. In this case, the null hypothesis is that all 
events are generated from the isotropic background distribution, i.e., ns = 0. The 
alternative hypothesis is that n, neutrinos originate from the source. As is the case 
for the diffuse flux, a harder spectrum can indicate a signal. The TS is calculated 
as: 
L(is,4) 
TS = 2log len, = 0: (22) 
The significance of an observation is determined by comparing the observed TS 
to the TS distribution from data sets randomized in right ascension. The TS distri- 
bution for randomized data sets represents the probability that a given observation 
could occur by random chance within the data set. For large sample sizes, this 
distribution approximately follows a chi-squared distribution, where the number of 
degrees of freedom corresponds to the difference in the number of free parameters 
between the null hypothesis and the alternate hypothesis. 
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In addition to triggered and untriggered searches for neutrino sources, stacking 
a collection of candidate sources could be an effective way to enhance the discovery 
potential. in stacking searches, the correlation to a catalog of sources is tested 
instead of searching in the direction of a single source. The stacking likelihood is 
defined as 


events M 


Ns Wij N- Ns 
L(ns, Y = II LN Po (lai — zil, 0%, Eiry) + — p 8il Bi), (23) 


where B; puede the isotropic background PDF, and S; the signal PDF, for 
each event. M is the number of sources in the catalog and wj the normalized 
theoretical weight for each source. This weight could be associated with properties 
of the individual sources in the catalog such as distance or the magnitude of their 
flux at some wavelength. 

Either the point source search or the stacking search can be modified in order 
to search for the extended sources of neutrinos. For this purpose, the angular 
uncertainty in the spatial PDF has to be modified using the extension of the source. 
In this case the effective angular uncertainty for an extended source is given by 


Oot = 4/0? +024 (24) 


The technique described above can be extended to a time-dependent search for 
transient neutrino sources.” A time component is added to the signal function in 
Eq. (19). This component will take into account the temporal structure of the neu- 
trino emission. This temporal dependency can be modeled as a Heaviside function 
or a Gaussian, for example. 

Time-dependent searches are generally more sensitive than time-integrated 
searches because they accumulate lower background rates. These searches, there- 
fore, offer an alternative opportunity to pinpoint the origin of cosmic neutrinos. 


4. The Discovery of Cosmic Neutrinos 


For neutrino astronomy, the first challenge is to select a pure sample of neutrinos, 
more than 100,000 per year above a threshold of 0.1 TeV for IceCube, in a back- 
ground of ten billion cosmic-ray muons, while the second is to identify the small 
fraction of these neutrinos that is astrophysical in origin, expected to be at the 
level of ten to hundred events per year according to the expectations of Fig. 6. 
Atmospheric neutrinos are a background for cosmic neutrinos, at least at neutrino 
energies below ~ 300 TeV. Above this energy, the atmospheric neutrino flux reduces 
to less than one event per year, even in a kilometer-scale detector, and thus events 
in that energy range are predominantly cosmic in origin. 

Searching for high-energy neutrinos of cosmic origin, IceCube continuously mon- 
itors the whole sky collecting very high statistics data sets of atmospheric neutrinos. 
Neutrino energies cover more than six orders of magnitude, from 5 GeV in the highly 
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Fig. 12. Left Panel: Light pool produced in IceCube by a shower initiated by an electron or tau 
neutrino of 1.14 PeV, which represents a lower limit on the energy of the neutrino that initiated 
the shower. White dots represent sensors with no signal. For the colored dots, the color indicates 
arrival time, from red (early) to purple (late) following the rainbow, and size reflects the number 
of photons detected. Right Panel: A muon track coming up through the Earth, traverses the 
detector at an angle of 11° below the horizon. The deposited energy, i.e., the energy equivalent of 
the total Cherenkov light of all charged secondary particles inside the detector, is 2.6 PeV. 


instrumented inner core (DeepCore) to beyond 10 PeV. Soon after the completion 
of the detector, with two years of data, IceCube discovered an extragalactic flux 
of cosmic neutrinos with an energy flux in the local universe that is, surprisingly, 
similar to that in gamma rays; see references 74 and 75. 

Two principal methods are used to identify cosmic neutrinos. The first method 
reconstructs upgoing muon tracks initiated by muon neutrinos and the second iden- 
tifies neutrinos of all flavors interacting inside the instrumented volume of the de- 
tector. For illustration, the Cherenkov patterns initiated by an electron (or tau) 
neutrino of 1 PeV energy and a neutrino-induced muon losing 2.6 PeV energy while 
traversing the detector are contrasted in Fig. 12. 

The kilometer-long muon range makes it possible to identify neutrinos that 
interact outside the detector and to separate them from the the background of at- 
mospheric muons using Earth as a filter. Using this method, IceCube has measured 
the background atmospheric neutrino flux over more than five orders of magnitude 
in energy with a result that is consistent with theoretical calculations. However, 
with eight years of data, IceCube has observed an excess of neutrino events at en- 
ergies beyond 100 TeV *® that cannot be accounted for by the atmospheric flux; 
see Fig. 13. Although the detector only records the energy of the secondary muon 
inside the detector, from Standard Model physics we can infer the energy spectrum 
of the parent neutrino. The high-energy cosmic muon neutrino flux is well described 
by a power law with a spectral index of 2.28+0.10 and a normalization at 100 TeV 
neutrino energy of (1.447933) x 10718 GeV—!cm~?sr-!.° 

The arrival directions of the muon tracks have been analyzed by a range of 
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Fig. 13. Distribution of muon tracks arriving from the Northern Hemisphere, i.e., with declination 
greater than -5°, as a function of median neutrino energy. The black crosses correspond to observed 
data. The blue band represents the expected conventional atmospheric neutrino. The red lines 
shows the high-energy cosmic neutrino flux with the best fit spectral index of 2.13. The black 
dashed line demonstrates the probability density function for the highest energy neutrino event 
when the best fit spectrum is assumed. 


statistical methods, yielding a first surprise: there is no evidence for any correlation 
to sources in the Galactic plane. IceCube is recording a diffuse flux of extragalactic 
sources. Only after analyzing 10 years of data®° has evidence emerged at the 3c 
level that the neutrino sky map might not be isotropic. The anisotropy results from 
four sources-TXS 0506+056 among them (you will hear more about that source 
later on), that emerge as point sources at the 4ø level (pretrial); see Fig. 14. The 
strongest of these sources is the nearby active galaxy NGC 1068, also known as 
Messier 77. 

NGC 1068 is one of the best-studied Seyfert 2 galaxies that motivated the AGN 
unification model. Its measured infrared luminosity implies a high level of star- 
burst activity. In addition, NGC 1068 has a heavily obscured nucleus and the 
column density of the gas surrounding it reaches to ~ 10?°cm~?.°! Therefore, the 
high-energy electromagnetic emission is absorbed in the Compton thick molecular 
gas, which even makes measuring its intrinsic X-ray emission challenging.*!*? As 
such, the optically thick environment at the core of NGC 1068 provides a favor- 
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Fig. 14. Upper limits on the flux from candidate point sources of neutrinos in 10 years of IceCube 
data assuming two spectral indices of the flux. Four sources exceed the 4c level (pretrial) and 
collectively result into a 30 anisotropy of the sky map. 


able environment for the efficient production of high-energy neutrinos. The ~51 
signal neutrinos identified in the direction of NGC 1068°° implies a neutrino flux 
that considerably exceeds the gamma-ray emission measured by the Fermi satellite, 
implying a very efficient neutrino emission concurrent with suppression of the very 
high energy gamma rays. 

Modeling of the high-energy neutrino emission from the corona of AG 
can accommodate these conditions. In Seyfert galaxies and quasars, a magnetized 
corona is formed above the accretion disk on the central black hole; see Miller & 
Stone (2000)°° for details. The hot, turbulent, and highly magnetized corona of 
AGN facilitates particle acceleration. Combined with the high density and abun- 
dance of target gas and radiation in the vicinity of the AGN core, a sizable neutrino 
flux can be expected.°° The accompanying pionic photons will cascade down due to 
the high opacity of the target. The efficient production of neutrinos in the vicinity 
of supermassive black holes by this mechanism can accommodate the measured flux 
of high-energy neutrinos. We will discuss this further in Sec. 6 where we highlight 
implications of astrophysical beam dumps for neutrino emission, especially in the 
presence of major accretion activity onto a supermassive black hole. 

The second method for separating cosmic from atmospheric neutrinos exclusively 
identifies high-energy neutrinos interacting inside the detector, so-called high-energy 
starting events. It divides the instrumented volume of ice into an outer veto shield 
and a ~ 500-megaton inner fiducial volume. The advantage of focusing on neutrinos 
interacting inside the instrumented volume of ice is that the detector functions as a 
total absorption calorimeter®’ allowing for a good energy measurement that sepa- 
rates cosmic from lower-energy atmospheric neutrinos. With this method, neutrinos 
from all directions in the sky and of all flavors can be identified, which includes both 
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Fig. 15. Left Panel: Deposited energies, by neutrinos interacting inside IceCube, observed in 
six years of data.”8 The grey region shows uncertainties on the sum of all backgrounds. The atmo- 
spheric muon flux (blue) and its uncertainty is computed from simulation to overcome statistical 
limitations in our background measurement and scaled to match the total measured background 
rate. The atmospheric neutrino flux is derived from previous measurements of both the 7, K, and 
charm components of the atmospheric spectrum.?? Also shown are two fits to the spectrum, as- 
suming a simple power-law (solid gray) and a broken power-law (dashed gray). Right Panel: The 
same data and models, but now showing the distribution of events with deposited energy above 
60 TeV in declination. At the South Pole, the declination angle 6 is equivalent to the distribution 
in zenith angle 8 related by the identity, 6 = 6 — 7/2. It is clearly visible that the data is flat in 
the Southern Hemisphere, as expected from the contribution of an isotropic astrophysical flux. 


muon tracks as well as secondary showers produced by charged-current interactions 
of electron and tau neutrinos and neutral current interactions of neutrinos of all 
flavors. A sample event with a light pool of roughly one hundred thousand photo- 
electrons extending over more than 500 meters is shown in the left panel of Fig. 12. 
The starting event sample revealed the first evidence for neutrinos of cosmic ori- 
gin.5®8? Events with PeV energies, and no trace of accompanying muons from an 
atmospheric shower, are highly unlikely to be of atmospheric origin. The present 
seven-year data set contains a total of 60 neutrino events with deposited energies 
ranging from 60 TeV to 10 PeV that are likely to be of cosmic origin. 

The deposited energy and zenith dependence of the high-energy starting events”* 
is compared to the atmospheric background in Fig. 15. The expected number of 
events for the best-fit astrophysical neutrino spectrum following a two-component 
power-law fit is shown as dashed lines in the two panels. The corresponding neutrino 
spectrum is also shown in Fig. 16. It is, above an energy of 200 TeV, consistent with 
a power-law flux of muon neutrinos penetrating the Earth inferred by the data shown 
in Fig. 13. A purely atmospheric explanation of the observation is excluded at 8c. 

Both measurements of the cosmic neutrino flux using cascades?! and throughgo- 
ing muons yield consistent determinations of the cosmic neutrino flux; see Fig. 17. 
The data are consistent with an astrophysical component with a spectrum close to 
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Unfolded spectrum for six years of HESE neutrino events starting inside the detector. 


The yellow and red bands show the 1ø uncertainties on the result of a two-power-law fit. Superim- 
posed is the best fit to eight years of the upgoing muon neutrino data (pink). Note the consistency 


of the red and pink bands. Figure from reference 78. 
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The flux of cosmic muon neutrinos’® inferred from the eight-year upgoing-muon track 
analysis (red solid line) with lo uncertainty range (shaded range; from fit shown in upper-right 
inset) is compared with the flux of showers initiated by electron and tau neutrinos.?! The measure- 
ments are consistent assuming that each flavor contributes the same flux to the diffuse spectrum. 
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Fig. 18. The flux of high-energy neutrinos (per flavor) above 1 TeV. At low energies, IceCube’s 
measured atmospheric muon and electron neutrinos??? are shown in light and dark grey, respec- 
tively. The modeled prompt flux®® from atmospheric charm production in cosmic ray showers is 
shown in black. The red line shows the fitted spectra for the HESE analysis covering 7.5 years 
of data and the shaded region shows the uncertainties.’ The dark red points show segmented 
flux data points in the HESE 7.5-year analysis. The results for the 10-year muon neutrino flux 
measurement is shown as the blue shaded area.’ At higher energies: the predicted flux of cosmo- 
genic neutrinos?” shown in green, as well as the upper limit on the extremely high energy flux of 


neutrinos from IceCube®® and Pierre Auger Observatory.?” 


In summary, IceCube has observed cosmic neutrinos using both methods for 
rejecting background. Based on different methods for reconstruction and energy 
measurement, their results agree, pointing at extragalactic sources whose flux has 
equilibrated in the three flavors after propagation over cosmic distances?’ with 
VeiYy iter iii. 

Figure 18 summarizes the measurements of the cosmic neutrino flux using start- 
ing events and throughgoing muons from the northern sky. An extrapolation of this 
high-energy flux to lower energy may suggest an excess of events in the 10— 100 TeV 
energy range over and above a single power-law fit. This conclusion is however not 
statistically compelling.”? 

In Fig. 19 we show the arrival directions of the most energetic events in the eight- 
year upgoing v,, +”, analysis (©) and the six-year HESE data sets. The HESE data 
are separated into tracks (@) and cascades (®). The median angular resolution of 
the cascade events is indicated by thin circles around the best-fit position. The 
most energetic muons with energy E, > 200 TeV in the upgoing vy, + v, data 
set accumulate near the horizon in the Northern Hemisphere. Elsewhere, muon 
neutrinos are increasingly absorbed in the Earth before reaching the vicinity of the 
detector because of their relatively large high-energy cross sections. This causes the 
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Fig. 19. Mollweide projection in Galactic coordinates of the arrival direction of neutrino events. 
We show the results of the eight-year upgoing track analysis” with reconstructed muon energy 


E,, 2 200 TeV (©). The events of the six-year high-energy starting event (HESE) analysis with 


deposited energy larger than 100 TeV (tracks ® and cascades @) are also shown.”*:100:10l The 


thin circles indicate the median angular resolution of the cascade events (@). The blue-shaded 
region indicates the zenith-dependent range where Earth absorption of 100 TeV neutrinos becomes 
important, reaching more than 90% close to the nadir. The dashed line indicates the horizon and 
the star (x) the Galactic Center. We highlight the four most energetic events in both analyses by 
their deposited energy (magenta numbers) and reconstructed muon energy (red number). 


apparent anisotropy of the events in the Northern Hemisphere. Also HESE events 
with deposited energy of Eaep > 100 TeV suffer from absorption in the Earth and 
are therefore mostly detected when originating in the Southern Hemisphere. After 
correcting for absorption by the Earth, the arrival directions of cosmic neutrinos 
are isotropic, suggesting extragalactic sources. In fact, no correlation of the arrival 
directions of the highest energy events, shown in Fig. 19, with potential sources or 
source classes has reached the level of 30.°° 

We should comment at this point that there is yet another method to conclu- 
sively identify cosmic neutrinos: the observation of very high energy tau neutrinos. 
Below an energy of 100 GeV, tau neutrinos are abundantly produced in the atmo- 
sphere by the oscillation of muon into tau neutrinos. Above that energy, they must 
be of cosmic origin, produced in cosmic accelerators whose neutrino flux has ap- 
proximately equilibrated between the three flavors after propagating over cosmic 
distances. Tau neutrinos produce two spatially separated showers in the detector, 
one from the interaction of the tau neutrino and the second one from its decay; 
the mean tau lepton decay length is A, = [cr ~ 50 m x (E,/PeV), where m, T, 
and E are the mass, lifetime, and energy of the tau, respectively. Two such can- 
didate events have been identified.1°? An event with a decay length of 17m and a 
probability of 98% of being produced by a tau neutrino is shown in Fig. 20. 

Yet another independent confirmation of the observation of neutrinos of cosmic 
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Fig. 20. Event view of a tau neutrino. The Cherenkov photons associated with the production 
and subsequent decay of the tau neutrino are identified by the double-peaked photon count as a 
function of time for the bright DOMs, for instance, the one shown in the top-right corner. The 
best fit (solid line) corresponding to a 17m decay length and is far superior to fits assuming a 
single electromagnetic or hadronic shower (dashed lines). 


origin appeared in the form of a Glashow resonance event shown in Fig. 21. The 
event was identified in a search for partially contained events in a 4.5 year data set: 
an antielectron neutrino interacting with an atomic electron produced an event com- 
patible with an incident neutrino energy of 6.3 PeV, characteristic of the resonant 
production of a weak intermediate W7 with a mass of My = 80.38 GeV.10° 
Given its energy and direction, the event is classified as an astrophysical neu- 
trino at the 50 level. Furthermore, data collected by the sensors closest to the 
interaction point, as well as the measured energy, are consistent with the hadronic 
decay of a W~ produced on the Glashow resonance. In the observer frame, where 
the electron mass (Mme = 0.511 MeV) is at rest, the resonance energy is given 
by Er = Mĝ,/2me = 6.32 PeV for Mw = 80.38GeV. The measured energy of 
6.05 + 0.72 PeV translates into a neutrino energy of 6.3 PeV after correcting the 
visible energy produced by the hadronic decay of the W for shower particles that do 
not radiate. Taking into account the detector’s energy resolution, the probability 
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that the event is produced off resonance by deep inelastic scattering is only 0.01 
assuming a spectrum with a spectral index of y = —2.5. Its observation extends 
the measured astrophysical flux to 6.3 PeV. Assuming the Standard Model resonant 
cross section, we expect 1.55 events in the sample assuming an antineutrino:neutrino 
ratio of 1:1 characteristic of a cosmic beam dump producing an equal number of 
pions of all three electric charges. 

The observation of a Glashow resonant event heralds the presence of electron 
antineutrinos in the cosmic neutrino flux. Its unique signature illustrates a method 
to disentangle neutrinos from antineutrinos, thus opening a path to distinguish 
astronomical accelerators that produce neutrinos via hadronuclear or photohadronic 
interactions with or without strong magnetic fields .‘°° As such, knowledge of both 
the flavor and charge of the incident neutrino will add a new tool for doing neutrino 
astronomy. 


5. Multimessenger Astronomy: General Considerations 


The most important message emerging from the IceCube measurements may not be 
apparent yet: the prominent and surprisingly important role of neutrinos relative 
to photons in the nonthermal universe. To illustrate this point, we show in Fig. 22 
the observed energy flux of neutrinos, Æ?ġ. One can see that the cosmic energy 
density of high-energy neutrinos is comparable to that of gamma-rays observed 
with the Fermi satellite” and to that of the ultra-high-energy (UHE) cosmic rays 
above 10° GeV, observed, by the Auger observatory.!04 
mon origin and, in any case, provides an excellent opportunity for multi-messenger 
studies. 

The pionic gamma rays accompanying the neutrino flux is shown as the solid 
blue. It has been derived from Eq. 10 but, for extragalactic neutrinos, this gamma- 
ray emission is not directly observed because of strong absorption of photons by 
ete” pair production in the extragalactic background light (EBL) and CMB. The 
high-energy photons initiate electromagnetic showers by repeated pair production 
and inverse-Compton scattering, mostly on CMB photons, that eventually yield the 
lower energy photons observed by Fermi in the GeV-TeV range. 

The extragalactic gamma-ray background observed by Fermi" has contributions 
from identified point-like sources on top of an isotropic gamma-ray background 
(IGRB) shown in Fig. 22. This IGRB is expected to consist mostly of emission 
from the same class of gamma-ray sources that are individually below Fermi’s point- 
source detection threshold (see, e.g., reference 105). The significant contribution 
of gamma-rays associated with IceCube’s neutrino observation has the somewhat 
surprising implication that indeed many extragalactic gamma-ray sources are also 
neutrino emitters, while none has been detected so far. 

Another intriguing observation is that the high-energy neutrinos observed at 
IceCube could originate in the sources of the highest energy cosmic rays. These 
could, for instance, be embedded in environments that act as “storage rooms” for 


This may indicate a com- 
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Fig. 21. Particle shower created at the Glashow resonance. The energy, measured via Cherenkov 
radiation in the Antarctic ice sheet, is reconstructed at the resonant energy for the production 
of a weak intermediate boson W7 in the interaction of an antielectron neutrino with an atomic 
electron in the ice. Features consistent with the production of secondary muons in the particle 
shower indicate the hadronic decay of a resonant W~ boson. 


cosmic rays with energies below the “ankle” (Ecr < 1EeV). This energy-dependent 
trapping can be achieved via cosmic ray diffusion in magnetic fields. While these 
cosmic rays are trapped, they can produce gamma-rays and neutrinos via collisions 
with gas. If the conditions are right, this mechanism can be so efficient that the 
total energy stored in low-energy cosmic rays is converted to that of gamma rays and 
neutrinos. These “calorimetric” conditions can be achieved in starburst galaxies!° 
or galaxy clusters.'°’ We will discuss these multimessenger relations in more detail 
next. 
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Fig. 22. The spectral flux (¢) of neutrinos inferred from the eight-year upgoing track analysis (red 
fit) and the six-year HESE analysis (magenta fit) compared to the flux of unresolved extragalactic 
gamma-ray sources74 (blue data) and ultra-high-energy cosmic rays!°* (green data). The neutrino 
spectra are indicated by the best-fit power-law (solid line) and 1ø uncertainty range (shaded range). 
We highlight the various multimessenger interfaces: A: The joined production of charged pions 
(x=) and neutral pions (7°) in cosmic-ray interactions leads to the emission of neutrinos (dashed 
blue) and gamma-rays (solid blue), respectively. B: Cosmic ray emission models (solid green) of 
the most energetic cosmic rays imply a maximal flux (calorimetric limit) of neutrinos from the 
same sources (green dashed). C: The same cosmic ray model predicts the emission of cosmogenic 
neutrinos from the collision with cosmic background photons (GZK mechanism). 


5.1. IceCube Neutrinos and Fermi Photons 


Recall that photons are produced in association with neutrinos when accelerated 
cosmic rays produce neutral and charged pions in interactions with photons or 
nuclei. Targets include strong radiation fields that may be associated with the 
accelerator as well as concentrations of matter or molecular clouds in their vicinity. 
Additionally, pions can be produced in the interaction of cosmic rays with the EBL 
when propagating through the interstellar or intergalactic background medium from 
their source to Earth. A high-energy flux of neutrinos is produced in the subsequent 
decay of charged pions via mt + yt + v, followed by wt > et + ve + 7, and 
the charge—-conjugate processes. High-energy gamma rays result from the decay of 
neutral pions, 7? + y +y. Pionic gamma rays and neutrinos carry, on average, 1/2 
and 1/4 of the energy of the parent pion, respectively. With these approximations, 
the neutrino production rate Q,,, (units of GeV~'s~) can be related to the one for 
charged pions as 


ye E Qu, (Ev) © 3 [ErQr+ (Er) E ~4E, $ (25) 


Similarly, the production rate of pionic gamma-rays is related to the one for neutral 
pions as 


E,Q4(Ey) ~ 2[ExQzo (Ex), ~2n, : (26) 
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Note, that the relative production rate of pionic gamma rays and neutrinos only 
depends on the ratio of charged-to-neutral pions produced in cosmic-ray interac- 
tions, denoted by K, = N,+/N,0. Pion production of cosmic rays in interactions 


with photons can proceed resonantly in the processes p + y > At —> 7° + p and 
p+7— At + nt +n. These channels produce charged and neutral pions with 
probabilities 2/3 and 1/3, respectively. However, the additional contribution of 
nonresonant pion production changes this ratio to approximately 1/2 and 1/2. 


In contrast, cosmic rays interacting with matter, for instance with hydrogen in 
the Galactic disk, produce equal numbers of pions of all three charges: p+ p —> 
Nn [10 +r +r] +X, where n, is the pion multiplicity. From above arguments we 
have K, œ 2 for cosmic ray interactions with gas (pp) and K, œ 1 for interactions 
with photons (py). 

With this approximation we can combine Eqs. (8) and (26) to derive a simple 
relation between the pionic gamma-ray and neutrino production rates: 


1 Ky 
5 BQ un (By) = GE [E3Qy(Ey)] poop, - (27) 


The prefactor 1/4 accounts for the energy ratio (E,,)/(E,) ~ 1/2 and the two 
gamma rays produced in the neutral pion decay. This powerful multimessenger 
relation connects pionic neutrinos and gamma rays without any reference to the 
cosmic ray beam; it simply reflects the fact that a n? produces two y rays for 
every charged pion producing a v,, + Y,, pair, which cannot be separated by current 
experiments. 

Before applying this relation to a cosmic accelerator, we have to be aware of 
the fact that, unlike neutrinos, gamma rays interact with photons of the cosmic 
microwave background before reaching Earth. The resulting electromagnetic shower 
subdivides the initial photon energy, resulting in multiple photons in the GeV- 
TeV energy range by the time the photons reach Earth. Calculating the cascaded 
gamma-ray flux accompanying IceCube neutrinos is straightforward.!0°:10° 

An example of the relation between pionic gamma-ray (solid) and neutrino 
(dashed) emission is illustrated in Fig. 22 by the blue lines. We assume that the 
underlying 7° / 7* production follows from cosmic-ray interactions with gas in the 


universe; K, ~ 2. In this way, the initial emission spectrum of gamma-rays and 
neutrinos from pion decay is almost identical to the spectrum of cosmic rays, after 
accounting for the different normalizations and energy scales. The flux of neutrinos 
arriving at Earth (blue dashed line) follows this initial CR emission spectrum, as- 
sumed to be a power law, E~7. However, the observable flux of gamma-rays (blue 
solid lines) is strongly attenuated above 100 GeV by interactions with extragalactic 
background photons. 

The overall normalization of the emission is chosen in a way that the model does 
not exceed the isotropic gamma-ray background observed by the Fermi satellite 
(blue data). This implies an upper limit on the neutrino flux shown as the blue 
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dashed line. Interestingly, the neutrino data shown in Fig. 22 saturates this limit 
above 100 TeV. Moreover, the HESE data that extends to lower energies is only 
marginally consistent with the upper bound implied by the model (blue dashed 
line). If the underlying assumptions hold, the neutrino spectrum cannot be harder 
than E~'°, a result somewhat challenged by the data. This example shows that 
multi-messenger studies of gamma-ray and neutrino data are powerful tools to study 
the neutrino production mechanism and to constrain neutrino source models.!!° 

Above exercise relating gamma rays and neutrinos comes with more warnings 
than a TV commercial for drugs: 


e The target for producing the neutrinos may be photons. This changes the 
value of K, and the shape of the flux. Yoshida and Murase recently worked 
through this example in detail.'!! 

e Electrons accelerated along with protons may contribute photons to the 
observed flux on top of those of pionic origin. 

e The source itself may not be transparent to high-energy photons that will 
lose energy in the source event before reaching the EBL. As a result, part 
of the pionic flux will emerge below the threshold of the Fermi satellite, at 
MeV and lower energies. This will be an important consideration when we 
discuss the first identified source of high-energy neutrinos, the black hole 
TXS 0506+056. 


While the details matter, the main message should not be lost: the matching 
energy densities of the extragalactic gamma-ray flux detected by Fermi and the 
high-energy neutrino flux measured by IceCube suggest that, rather than detecting 
some exotic sources, it is more likely that IceCube to a large extent observes the 
same universe astronomers do. Clearly, an extreme universe modeled exclusively 
on the basis of electromagnetic processes is no longer realistic. The finding implies 
that a large fraction, possibly most, of the energy in the nonthermal universe orig- 
inates in hadronic processes, indicating a larger role than previously thought. The 
high intensity of the neutrino flux below 100 TeV in comparison to the Fermi data 
might indicate that these sources are even more efficient neutrino than gamma-ray 
sources, 12113 

IceCube is developing methods that enable real-time multiwavelength observa- 
tions in cooperation with astronomical telescopes, to identify the sources and build 
on the discovery of cosmic neutrinos to launch a new era in astronomy.'!*!!° We 
will return to a coincident observation of a flaring blazar on September 22, 2017, 
further on. 


5.2. IceCube Neutrinos and Ultra-High-Energy Cosmic Rays 


The charged pion production rate Q,+ is proportional to the density of the cosmic- 
ray nucleons in the beam, Qw, by a “bolometric” proportionality factor f,;. For a 
target with nucleon density n and extension £, the efficiency factor for producing 
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pions is fr ~ 1 — exp(—Klon), where « is the inelasticity, i.e., the average relative 
$ . 
its 


interaction length with o the cross section for either py or pp interactions. The pion 


energy loss of the leading nucleon going into the production of pions, (an)7 


production efficiency f» normalizes the conversion of cosmic-ray energy into pion 
energy on the target as: 
K. 
2 5 T 2 

EzQrt (Ex) = fr EEA [ExQn(En)| En=Ex/[ty ` (28) 
We already introduced the pion ratio K, in the previous section, with Kr ~ 2 for 
pp and K, ~ 1 for py interactions. The factor x, denotes, as before, the average 
inelasticity per pion that depends on the average pion multiplicity ny. For, both, 


pp and py interactions this can be approximated as 2, = K/n, ~ 0.2. The average 
energy per pion is then (Er) = «,Ey and the average energy of the pionic leptons 
relative to the nucleon is (Ep) ~ (Eq) /4 = (a,7/4)En ~ 0.05En. 

In general, the cosmic ray nucleon emission rate, Qu, depends on the compo- 
sition of the highest energy cosmic rays and can be obtained by integrating the 
measured spectra’ of nuclei with mass number A as Qn (Ew) = >), A?Q4(AEn). 
In the following we will derive a upper limit on diffuse neutrino fluxes under the 
assumption that the highest energy cosmic rays are mostly protons.!!°!!” The local 
emission rate density, Q = poQ, is at these energies insensitive to the luminosity 
evolution of sources at high redshift and can be estimated to be at the level of 
[E2 Qp(Ep)] so10.5ey ~ (0-5 — 2.0) x 10“4erg/Mpc? /yr.?5118:119 Note, that composi- 
tion measurements indicate that the mass composition above the ankle also requires 
a contribution of heavier nuclei. However, the estimated local power density based 
on proton models is a good proxy for that of cosmic ray models that include heavy 
nuclei, as long as the spectral index is close to y ~ 2. For instance, a recent anal- 
ysis of Auger!”° provides a solution with spectral index y ~ 2.04 and a combined 
nucleon density of [E?,On(Ey)]1019-5ev ~ 2.2 x 10 erg/Mpc?/yr. 

We construct the diffuse neutrino flux from the contribution of individual 
sources. A neutrino point-source (PS) at redshift z with spectral emission rate 
Qu, contributes a neutrino flux summed over flavors (in units GeV~'cm~?s~!) 


PS ey 
OB) = Fay Z Q(U + Ei), (29) 
where dy is the luminosity distance 
dz’ 


0 
Here, the Hubble parameter H has a local value of c/Ho ~ 4.4 Gpc and scales with 
redshift as H?(z) = HÊ[(1 + 2)? Qm + Qa], with Qn ~ 0.3 and Qa ~ 0.7 assuming 
the standard ACDM cosmological model.'*! Note that the extra factor (1 + z)? 


Note that the integrated number of nucleons is linear to mass number, f dEQn(£) = 


wy AS deQa(£). 
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in Eq. (29) follows from the definition of the luminosity distance and accounts 
for the relation of the energy flux to the differential neutrino flux ¢. The diffuse 
neutrino flux from extragalactic sources is given by the integral over co-moving 
volume dV. = 4m(dz/(1 + z))?dz/H(z). Weighting each neutrino source by its 
density per co-moving volume p(z) gives!?? 


c [®© dz 
WEE OOE +E), (31) 
and, assuming that Q,, follows a power law E77, 
1 2 C fz 1 2 
3 >, Epua (E) = H3 32 Ez Qva (Ev), (32) 


where Q,, = poQ,,, is the neutrino emission rate density and 
tT ale) 
0 Va + (1+ 2z)2Qm Po 


A spectral index of y ~ 2.0 and no source evolution in the local (z < 2) universe, 
p(z) = po, yields €, ~ 0.5. For sources following the star-formation rate, p(z) = 
(1+ z)? for z < 1.5 and p(z) = (1+ 1.5)? for 1.5 < z < 4, with the same spectral 
index yields €, ~ 2.6. 

We can now derive the relation between the diffuse neutrino flux from the cosmic 
ray rate density by combining Eqs. (26) and (28): 


. (34 
2.6 1044 erg/Mpc?/yr cm? ssr (34) 


2 
E E? 6, (Ep) ~ 3 x 107° fr ( &z ) (= eee | GeV 

Q 

Here, we have assumed pp interactions with K, = 2. The calorimetric limit, 
previously discussed corresponds to the assumption that f+ — 1, and is also referred 
to as the Waxman-Bahcall “bound”.!!°!!" Sources with f > 1 are referred to as 
“dark” sources; their thick target efficiently converts nucleons to neutrinos and 
rendering them opaque to high-energy photons. Accelerator-based beam dumps are 
the ultimate dark sources of neutrinos. 

It is intriguing that the observed intensity of diffuse neutrinos indicates that 
fx = 1. The correspondence is illustrated by the green lines in Fig. 22 showing 
a parametrization that accounts for the most energetic cosmic rays (green data). 
Note, that the cosmic ray data below 1019 GeV are not described by this fit and 
must be supplied by additional sources, e.g., in our own Galaxy. Assuming that 
the cosmic ray energy density is transformed in neutrinos we derive the maximal 
neutrino emission (green dashed line). Interestingly, the observed neutrino flux 
saturates this calorimetric limit. It is therefore possible that the highest energy 
cosmic rays and neutrinos have a common origin. If this is the case, the neutrino 
spectrum beyond 200 TeV should reflect the energy-dependent release of cosmic 
rays from the calorimeters. Future studies of the neutrino spectrum beyond 1 PeV 
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can provide supporting evidence for this and, in particular, the transition to a thin 
environment (fr < 1), that is a necessary condition for the emission of the highest 
energy cosmic rays by the source, implies a break or cutoff in the neutrino spectrum. 

Note that the proton model in Fig. 22 also contributes to the flux of EeV neu- 
trinos shown as a dotted green line. Ultra-high-energy cosmic rays are strongly 
attenuated by resonant interactions with background photons, as first pointed out 
by Greisen, Zatsepin and Kuzmin!?®1?4 (GZK). This GZK mechanism is responsi- 
ble for the suppression of the proton flux beyond 5 x 101° GeV (“GZK cutoff”) in 
Fig. 22 (green solid line) and predicts a detectable flux of cosmogenic neutrinos”® 
(green dotted line) in a kilometer-scale detector. In case a significant fraction of 
the highest energy cosmic rays are heavy nuclei as indicated by the Auger data, 
the GZK flux is reduced accordingly because the individual nucleons fall below 
the threshold for producing neutrinos. This may explain the failure of IceCube to 
observe GZK neutrinos at the level predicted in Fig. 22. 


5.3. Pinpointing the Astrophysical Sources of Cosmic Neutrinos: 
Relating the Point Source and Diffuse Fluxes 


The diffuse flux of neutrinos measured by IceCube provides a constraint on the flux 
from the individual sources that it is composed of. In the sense of Olbers’ paradox, 
one can investigate under what circumstances IceCube can detect the neutrino 
emission from individual, presumably nearby, point sources that contribute to the 
quasi-diffuse emission. 

One can relate the average luminosity of individual neutrino sources to the 
diffuse flux 


5 D Elhu (E) > x~ Ein 3 2o EsQva (E v) ; (35) 


a 

where the €, summarizes the evolution i the US, Ho is the Hubble constant, 
and po is the effective source density. 

Equation (35) relates the average luminosity of individual neutrino sources to 
the diffuse flux on the left hand side, which is measured by the experiment to be 
at the level of E?¢, ~ 1078 GeVcm~?s~'sr~! for energies in excess of ~ 100 TeV; 
see Fig. 18. From the measurement, we can infer the average emission from a single 
source!” as a function their local density, 


1 2 me 45 ( & E po - eae 
3 DF, (By) = 1.7% 10 (š 10 Mpe” ergs. (36) 


This relation builds in the very strong assumption that all sources corresponding 


to some astronomical classification are identical. 

For a homogenous distribution of sources, we expect, within the partial field of 
view of the full sky, fsky, of the telescope, one source within a distance dı determined 
by 

fi Ard? /3 p0 = 1. (37) 
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Fig. 23. The effective local density and (maximal) neutrino luminosity of various neutrino 
source candidates from reference 125. The green solid (green dotted) line shows the local den- 
sity and luminosity of the population of sources responsible for the diffuse neutrino flux of 
E?¢ ~ 1078 GeV cm~? s~! sr! observed with IceCube, assuming source evolution following 
the star-formation rate (£; ~ 2.6) or no source evolution (£z ~ 0.5), respectively. The grey-shaded 
area indicates source populations that are excluded by the non-observation of point sources in the 
Northern Hemisphere (fsky ~ 0.5) with discovery potential E?6P8 ~ 2x 1071? TeV cm7? s~!,!26 


In other words, dı defines the volume containing one nearby source for a homoge- 
neous source density po. Defining ¢ as the flux of a source at distance d1, given by 
Eq. (29) with z ~ 0 and dz = dj, we can write the probability distribution p(¢) of 
finding the closest source of the population with a flux ¢ as 
3 
31 Qı 2 (ey 
po) = $= @ eG) (38) 


(This result takes into account fluctuation in the source distribution; for details 
see appendix of reference 122). The average flux from the closest source is then 
(od) œ 2.7¢1, with median mea ~ 1.3461. Applying this to the closest source 
introduced above, we obtain its neutrino flux: 
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2 —1 pans 
1 = sk 5 z 0 : E: ees 
3 2 Biv. ~ 2 x 1071? (=) (5) (2) TeV cm~? s7}, 

(39) 
Interestingly, this value is not far from IceCube’s average point-source discovery 

potential at the level of 2 x 10712 TeV cm7? s7! in the Northern Hemisphere.*° 
In Fig. 23, we show the local density and luminosity of theorized neutrino 
125 The grey-shaded region is excluded by the failure to observe these 
sources as individual point sources. The green lines show the combination of den- 


sources. 


sity and luminosity for which sources produce a diffuse neutrino flux at the level 
observed by IceCube. We show the result for a source densities following the star 
formation rate (solid line) or no evolution (dotted line). We conclude that IceCube 
is presently sensitive to source populations with local source densities smaller than, 
conservatively, 1078 Mpe~*. Much lower local densities, like FSRQs, are challenged 
by the nonobservation of individual sources. Some source classes, like Fanaroff-Riley 
(FR) radio galaxies, have an estimated neutrino luminosity that is likely too low for 
the observed flux. This simple estimate can be refined by considering not only the 
closest source of the population, but the combined emission of known local sources; 
see, e.g., reference 122. 

The luminosity-density parameter space discussed above is often exploited to 
study sources identified in astronomical catalogs as potential contributors to the 
cosmic neutrino flux. However, it is critical to note that astronomical classifications 
are meant to explain the collective behavior of sources with regard to common spec- 
tral features and power. Thus, they may lack the necessary condition for a neutrino 
source: provide a beam dump with sufficient target density to produce neutrinos. 
Spectral features, deviating from a simplified assumption of the single power law, 
together with the variations of the power and target densities, modifies the allowed 
parameter space. In addition, note that the source luminosity distribution is ne- 
glected by assuming an average luminosity for all sources. Finally, the evolution of 
the sources plays a key role in the allowed parameter space. As such, the built-in 
assumption of a specific cosmological evolution that is common to all sources is 
problematic. 

In addition to the strong dependence of these bands on the spectrum of the 
sources and their evolution, we should note that the all-sky untriggered searches 
for the sources of high-energy neutrinos allocate a large number of trial factors that 
penalizes the significance for any potential signal excess. Therefore, the presence 
of outlying sources in the neutrino data could result in the violation of such bands 
when additional information provided by other messengers distinguishes a specific 
source or direction in the sky. 

Is it possible that the sources of the extragalactic cosmic rays are themselves 
neutrino sources? We previously derived the emission rate density of nucleons from 
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the measured cosmic-ray spectrum?” 118 
Lp = po EZQp(Ep) = (1 — 2) x 10“ erg Mpc? yr. (40) 
Combining this with Eq.28 we can derive the diffuse neutrino flux: 
1 €,K, z -2 — 
3 2 Popua (Er) = fe ae 274) x 10 8 GeV cm~? s7! sr. (41) 


The equation has been rewritten and some notation adjusted in order to accommo- 
date both pp and yp interactions. Counting particles we derive that 


1 
3 5 E Qua (EL) = E:Qr ; (42) 


from which the relation for energy follows by multiplying both sides with F,; an 
additional factor of 1/4 multiplies the right-hand side accounting for the ratio of 
neutrino to pion energy within the approximations routinely used throughout. 

The requirement f, < 1 may limit the neutrino production by the actual sources 
of the cosmic rays as pointed out by the seminal work by Waxman and Bahzcall.!!® 
For optically thin sources, fr < 1, neutrino production is only a small by-product of 
the acceleration process. The energy loss associated with pion production must not 
limit the sources’ ability to accelerate the cosmic rays. On the other hand, optically 
thick sources, fr ~ 1, are more efficient neutrino emitters. Realistic sources of this 
type need different zones, one zone for the acceleration process (fr < 1) and a 
second zone for the efficient conversion of cosmic rays to neutrinos (fr ~ 1). An 
example for this scenario are sources embedded in starburst galaxies, where cosmic 
rays can be stored over sufficiently long timescales to yield significant neutrino 
production. A more relevant example is the multimessenger source TXS 0506+056 
which we will discuss in detail further on. 

For €, ~ 2.4 and K, œ 1-2, the upper bound resulting from Eq. (41) and 
fxr = 1 is at the level of the neutrino flux observed by IceCube. Therefore, it is 
possible that the observed extragalactic cosmic rays and neutrinos have the same 
origin. A plausible scenario is a “calorimeter” in which only cosmic rays with energy 
below a few 10 PeV interact efficiently. An energy dependence of the calorimetric 
environment can be introduced by energy-dependent diffusion. If D(E) is the dif- 
fusion coefficient, then the timescale of escape from the calorimeter is given by the 
solution to 6D(E)t = d?, where d is the effective size of the region. Typically, 
we have D(E) x E° with 6 ~ 0.3 — 0.6. In the following, we again consider the 
case of protons. Taking opp ~ 8 x 107° cm? at Ep = 100 PeV and the diffusion 
coefficient of D(Ep) ~ Daey (Ep/1GeV)!/3, the pp thickness can be expressed as 
Tpp © ClNgasT pp OF 


d 4 Deev =i E S13 n 
~ 0.18 p —” _ ). 43 
"pp (az) (aez) (o a (wax) (43) 


Here, we have used feasible parameters of starburst galaxies.!°°!!° Therefore, de- 
pending on the calorimetric environment, it is possible that the flux below a few 
PeV is efficiently converted to neutrinos and contributes to the TeV—PeV diffuse 
emission observed by IceCube. 
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6. The First-Identified Cosmic-Ray Accelerator: the Rotating Su- 
permassive Black Hole TXS 0506+056 


The qualitative matching of the energy densities of photons and neutrinos, dis- 
cussed in the previous section, may suggest that the unidentified neutrino sources 
contributing to the diffuse flux might have already been observed as strong gamma- 
ray emitters. Theoretical modeling!?”!?* and recent data analyses!??-!?! show that 
Fermi’s extragalactic gamma-ray flux is dominated by blazars. However, a dedi- 
cated IceCube study!*? of Fermi-observed blazars shows no evidence of neutrino 
emission from these source candidates. The inferred limit on their quasi-diffuse flux 
leaves room for a significant contribution to IceCube’s diffuse neutrino flux at the 
10% level, and increasing towards PeV. The multimessenger campaign launched by 
the neutrino alert IC170922 not only identified the first source of cosmic neutrinos, 
it will also shed light on this apparent contradiction. 

Since 2016, the IceCube multimessenger program has grown from Galactic 
supernova alerts and early attempts to match neutrino observations with early 
LIGO/Virgo gravitational wave candidates to a steadily expanding set of automatic 
filters that selects in real time rare very high energy events that are potentially cos- 
mic in origin. Within less than a minute of their detection in the deep Antarctic 
ice, the astronomical coordinates are sent to the Gamma-ray Coordinate Network 
for potential follow-up by astronomical telescopes. 

On September 22, 2017, the tenth such alert, IceCube-170922A,!°° reported 
a well-reconstructed muon that deposited 180 TeV inside the detector, and with 
an estimated neutrino energy of 290 TeV. Its arrival direction aligned with the 
coordinates of a known blazar, TXS 0506+056, to within 0.06°. The source was 
“flaring” with a gamma ray flux that had increased by a factor 7 in recent months. 
A variety of estimates converge on the probability of order 1078 that the coincidence 
is accidental. The identification of the neutrino with the source reaches the level 
of evidence, but no more. What clinches the association is a series of subsequent 
observations, culminating with the optical discovery of a switch of the source from an 
off to an on state two hours after the emission of IC190922, conclusively associating 
the neutrino with TXS 0506+056: 


e The redshift of the host galaxy was measured to be z ~ 0.34.14 It is 
important to realize that nearby blazars like the Markarian sources are at 
a redshift that is ten times smaller, and therefore TXS 0506+056, with 
a similar flux despite the greater distance, is one of the most luminous 
sources in the Universe. It must belong to a special class of sources that 
accelerate proton beams in dense environments, revealed by the neutrino. 
The source is somehow special, and this eliminates any conflict between 
its observation and the lack of correlation between the arrival directions of 
IceCube neutrinos and the bulk of the blazars observed by Fermi.!°? 

e Originally detected by NASA’s Fermi!’ and Swift'*° satellite telescopes, 
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the alert was followed up by the MAGIC and other air Cherenkov tele- 
scopes.!°” MAGIC detected the emission of gamma rays with energies 
exceeding 100 GeV starting five days after the observation of the neutrino, 
indicating that the source is a relatively rare TeV blazar given the large 
redshift of the source. After correcting for the absorption on the microwave 
background, we conclude that the source is a TeV blazar. 

Given where to look, IceCube searched its archival neutrino data up to and 
including October 2017 for evidence of neutrino emission at the location 
of TXS0506+056. When searching the sky for point sources of neutrinos, 
two analyses have routinely been performed: one that searches for steady 
emission of neutrinos and one that searches for flares over a variety of 
timescales. Evidence was found for 19 high-energy neutrino events on a 
background of less than 6 in a burst lasting 110 days. This burst dominates 
the integrated flux from the source over the last 9.5 years for which we have 
data, leaving the 2017 flare as a second subdominant feature. We note that 
this analysis applied a published prescription to data; the chances that this 
observation is a fluctuation are small. 

Radio interferometry images of the source revealed a jet that loses its tight 
collimation beyond 5 milliarcseconds running into target material that cre- 
ates the opportunity for producing the neutrinos. The origin of this target 
material is still a matter of debate, and speculations include the merger with 
another galaxy affecting the jet from the dominant supermassive black hole. 
Alternatively, in a structured jet, the accelerated protons may collide with 
a slower moving and denser region of jetted photons. The jet may interact 
with the dense molecular clouds of a star-forming region, or simply with su- 
permassive stars in the central region of the host galaxy.13®139 The VLBA 
observations not only identify a striking feature in the jetted material, but 
also reveal that the neutrino burst occurs at the peak of enhanced radio 
emission at 15 GHz, which started five years ago, see Fig. 24. The radio 
flare may be a signature of a galaxy merger; correlations of radio bursts with 
the process of merging supermassive black holes have been anticipated.!*° 
The robotic optical telescope MASTER network has been monitoring the 
source since 2005 and found the strongest time variation of the source over a 
period of two hours after the emission of IC170922, with a second variation 
following the 2014-15 burst.!*! The blazar switches from the off to the 
on state two hours after the emission of the neutrino. The observation 
conclusively associates the source with the neutrino. 


Additionally, it is important to note the striking fact that all high-energy spectra, 
for both photons and neutrinos and for both the 2014 and 2017 bursts, are consistent 
with a hard E~? spectrum, which is expected for a cosmic accelerator. In fact, the 
gamma-ray spectrum hints at flattening beyond that during the 110-day period of 
the 2014 burst. 14143 
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Fig. 24. TXS 0506+056 radio light curve from Owen Valley Radio Observatory (OVRO) at 15 
GHz (red). THe dashed line illustrates the pattern of the radio flux density. The 2014/15 110-day 
neutrino flare (yellow band) and the IceCube-170922A episodes are shown. Radio data indicate 
that the neutrinos arrive during periods of enhanced radio emission. 


The problem with modeling TXS 05060+056 as a blazar is that it is not a 
blazar at the times of neutrino emission as shown by the multimessenger observation 
itemized above. It is evident that a source that is transparent to high-energy gamma 
rays is unlikely to host the target material to produce neutrinos with the opacity 
for yy interactions two orders of magnitude larger than that for py interactions. 
We will work through these arguments quantitatively further on. 

In this context, an alert recorded by IceCube on July 30, 2019, provides striking 
support for the idea that cosmic neutrinos are produced by temporarily gamma- 
suppressed blazars. IC190730 and IC170922 are the two highest energy alerts 
recorded so far. A well-reconstructed 300-TeV muon neutrino is observed in spatial 
coincidence with the blazar PKS 1502+106.'*4 OVRO radio observations!*° show 
that the neutrino is coincident with the highest flux density of a flare at 15 GHz that 
started five years ago,'*° matching the similar long-term radio outburst observed 
from TXS 0506+056. Even more intriguing is that the gamma ray flux observed by 
Fermi shows a clear minimum at the time that the neutrino is emitted; see Fig. 25. 
At the time that the target crosses the jet and produces a neutrino beam, the high- 
energy gamma rays are absorbed and their energy cascades to energies below the 
Fermi threshold, i.e., MeV or lower. For further discussion, see reference 147. 

Other IceCube alerts have triggered interesting observations. Following up on a 
July 31, 2016, neutrino alert, the AGILE collaboration, which operates an orbiting 
X-ray and gamma-ray telescope, reported a day-long blazar flare in the direction 
of the neutrino one day before the neutrino detection.'** Even before IceCube 
issued automatic alerts, in April 2016, the TANAMI collaboration argued for the 
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Fig. 25. Temporal variation of the y-ray and radio brightness of PKS 1502+106. a: Fermi- 


LAT likelihood light curve integrated between 100 MeV and 300 GeV (marked by black dots with 
error bars). b: OVRO flux density curve of PKS 1502+106 plotted with light blue dots, that is 
superimposed by the radio flux density curve binned to the Fermi-LAT light curve (marked with 
dark blue squares). The detection time of the neutrino IC-190730A is labeled by a vertical purple 


line. 


association of the highest energy IceCube event at the time, dubbed “Big Bird,” 
with the flaring blazar PKS B1424-418.'*° The event was produced at a minimum 
of the Fermi flux,!*” as was the case for PKS 1502 +106. AMANDA, IceCube’s 
predecessor, observed three neutrinos in coincidence with a rare flare of the blazar 
1ES 1959+650, detected by the Whipple telescope in 2002.1°° However, none of 
these identifications reach the significance of the observations triggered by IceCube- 


170922A. 


Note that a subset of blazars, around 1 ~ 10% of all blazars, bursting once in 10 
years at the levels of TXS 0506+056, can accommodate the diffuse cosmic neutrino 
flux observed by IceCube. The energy of the neutrino flux generated by the flaring 
blazars is also at the same level as the flux in extragalactic cosmic rays. In order to 
calculate the flux of the high-energy neutrinos from a population of sources, we use 
the formalism developed in the previous section® to relate the diffuse neutrino flux 
to the energy injection rate of the cosmic rays and their efficiency to transfer this 
energy from protons to neutrinos. For a population of sources with source density 


p and neutrino luminosity Lp, the diffuse neutrino flux is 


ped — xfer 
dE 4r 


(44) 


page 49 
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which can be simplified into 
2dN ž € 
dE 4r 
where € it the result of the integration over the redshift history of the sources; see 
Eq. 68. We revise the equation to account for the duration of the flares At, the 
total time of observation Tobs, and the fraction of sources F: 


24AN c 


At 
E pn 128 


corresponding to 


7 oa F Ry £ Ly 
1 11 T 2 1 1 — 
3 x 10 eVcm îs “sr dg \ 1.3 Gpe 0.7 1.2 x 10% erg/s 


p At 10yr 
47 
(ar) Ta) ee 


which results in F = 0.05. In summary, a special class of BL Lac blazars with a 
neutrino luminosity similar to TXS 05060+056, can accommodate the diffuse flux 
of high-energy cosmic neutrinos observed by IceCube. These high-energy neutrino 
flaring sources constitute 5% of the sources. 

The energetics in neutrino production from these sources also matches the energy 


flux of the highest energy cosmic rays, i.e., 


d d 
P ayi (50 one), (48) 


From the observed cosmic ray injection rate of (1 — 2) x 1044 erg Mpc~#yr~! we can 


determine the average pion efficiency of the neutrino sources that makes the match 
of Eq. 47 possible: 


Ly p At 10yr F 
1.2 x 1047 erg/s 1078 Mpc-3 110d Tops 0.05 


dE /dt 
~ =(1—e 7) “a : 
(1 — 2) x 1044 erg Mpe~3yr-! 


(49) 


We find that fr > 0.8. 

In the section on blazar jets, we discussed how this can be achieved with a proton 
beam with low boost factor interacting with the blue photons in the active galaxy. 
The jet producing the neutrinos is not transparent to TeV photons, only to photons 
with tens of GeV; this is indeed what is observed by Fermi at the time of the 2014 
flare. 

That TXS 0506+056 belongs to a special class of sources is reinforced by the 
fact that conventional blazar modeling of the multiwavelength spectrum has been 
unsuccessful despite the opportunity of tuning 14 free parameters. Blazar models at 
best meet the requirements for the 2017 burst by invoking Eddington biast! for the 
single neutrino produced. The answer to this puzzle could be a new subclass of high- 
energy gamma-ray sources, labeled as blazars by astronomers, that is responsible 
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for producing the high-energy cosmic neutrinos and cosmic rays observed.'°? In 
fact, if every source labelled as a blazar produced neutrinos at the level of TXS 
0506+056, they would overproduce the total diffuse flux observed by IceCube by 
close to two orders of magnitude. 

The observation that the energy content in neutrinos and very high energy 
cosmic rays are similar, underscores the fact that the cosmic rays must be highly 
efficient at producing neutrinos, requiring a large target density that renders them 
opaque to high-energy gamma rays. A consistent picture emerges when the source 
opacity exceeds a value of 0.4,'°? resulting in a gamma-ray cascade where photons 
lose energy in the source before cascading to even lower energies in the extragalactic 
background light. Most of their energy falls below the Fermi threshold by the time 
they reach Earth. The 2014-15 burst cannot be, and is not, accompanied by a large 
gamma-ray flare. This is consistent with the premise that a special class of efficient 
sources is responsible for producing the high-energy cosmic neutrino flux seen by 
IceCube. 

The nature of the special class of sources has not been settled. One straight- 
forward explanation could be a subclass of blazars selected by redshift evolution: 
powerful proton accelerators producing neutrinos may have been active in the past 
but are no longer today. This accommodates the large redshift of TKS 0506-+056,!°° 
which would be the closest among a set of sources that only accelerated cosmic rays 
at early redshifts. It is a requirement that the sources must have a sufficient target 
density in photons or protons in order to produce neutrinos at the high flux level 
observed from TXS 0506+056.!°? If neutrino sources only represent a few percent 
of those labelled by astronomers as blazars, attempts by IceCube and others to find 
correlations between the directions of high-energy neutrinos and all Fermi blazars 
must inevitably be unsuccessful. 

Merger activity in blazars in not uncommon. In merging galaxies there is plenty 
of material for accelerated cosmic rays to interact with the jet. This fresh material 
provides optically thick environments, and allows for rapid variation of the Lorentz 
factors. A cursory review of the literature on the production of neutrinos in galaxy 
mergers is sufficient to conclude that it can indeed accommodate the observations of 
both the individual sources discussed above and the total flux of cosmic neutrinos; 
see, e.g., references 154-156. 

There is another, possibly related, development associated with the search for 
steady neutrino sources that recently delivered a new time-integrated neutrino sky 
map covering ten years of IceCube data.!°’ Evidence at the 3c level suggests that 
the neutrino sky map might no longer be isotropic. The anisotropy results from 
four sources-TXS 0506+056 among them-that show evidence for clustering at the 
4c level. The strongest of these sources is the nearby active Seyfert galaxy NGC 
1068 (Messier 77) with starburst activity. There is evidence for shocks near the core 
and for molecular clouds with densities of more than 10° cm~’. Similarly to TXS 
0506+056, a merger onto the black hole is observed - either with a satellite galaxy 
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or, more likely, with a star-forming region.'!°> This major accretion event may be 
the origin of the increased neutrino emission. 

The assumption that the evolution of most radio galaxies is governed by merg- 
ers14>159 supports the proposal of a redshift evolution selection of a special class 
of sources to account for the cosmic neutrino flares observed by the IceCube Col- 
laboration. Alternative scenarios proposed to describe the evolution of the radio 
galaxies based on their morphology, such as twin active galactic nuclei!®° and back- 
flow, 161 
flares. In this context, previous studies of the neutrino emission from radio galaxies 
are also relevant as they require a high density close to the source.*!:!°? Further- 
more, far-infrared observations support the argument that the central activity of 


the galaxies, such as starburst activity, is fueled by the concentration of matter that 
163 


also provide a potential production mechanism for the cosmic neutrino 


may be preceded by a merger. 

In a major merger scenario, when two galaxies with supermassive black holes 
at their centres with different spins merge, a spin-flip is expected to occur due to 
the spin-orbit precession and energy dissipation from gravitational radiation. !4®164 
This is accompanied by precession of the jet, and the activity may therefore not 
always look like a blazar because the jet might be only briefly pointing into our 
direction. It has been suggested that the precessional phase of the merger of two 
black holes, occurring prior to the spin-flip, is visible as a superdisk in radio galaxies 
where the precessing jet emerges as a superwind separating the radio lobes in the 
final stages of the merger.!°° 

As discussed in the previous sections, the absence of a strong anisotropy of neu- 
trino arrival directions raises the possibility that the cosmic neutrinos originate from 
a number of relatively weak extragalactic sources. It is indeed important to keep in 
mind that the interaction rate of a neutrino is so low that it travels unattenuated 
over cosmic distances through the tenuous matter and radiation backgrounds of 
the Universe. This makes the identification of individual point sources contributing 
to the IceCube flux challenging. !??166-168 
that IceCube is capable of localizing the sources by observing multiple neutrinos 
originating in the same location. Not having observed strong neutrino clusters in 


Even so, it is also important to realize 


the present data raises the question of how many events are required to make such 
a model-independent identification possible. The answer to this question suggests 
the construction of a next-generation detector that instruments a ten times larger 


volume of ice.!°? 
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7. Modeling Cosmic Beam Dumps 


In this section we will introduce the generic theoretical framework used to calcu- 
late the neutrino flux associated with an astronomical source where protons are 
accelerated to energies exceeding the threshold for producing pions that decay into 
neutrinos. We will work through several examples. We start with the generic beam 
dump where protons accelerated near the black hole hole at the center of an active 
galaxy interact with the gas (protons) in their vicinity, i.e the equivalent of the 
1 ~ proton/cm? in our own galaxy. Next, we will discuss two typical blazar models 
where accelerated protons interact with photons associated with the jet, or, alter- 
natively, with strong external radiation fields near the black hole. Finally, we turn 
from jets to fireballs and discuss the production of neutrinos when protons interact 
with photons in the relativistically expanding fireball following a stellar collapse, 
such as in a gamma-ray burst. 

Neutrinos are produced when pions, and, at higher energies, kaons and charm 
particles decay. We start by calculating the number of pions produced at the source 
when an accelerated proton beam interacts with gas of density n in the rest of the 
galaxy. We introduce the source function qr (Er) defined such that qr (Er) dEr 
represents the rate at which pions are produced within the energy range E, and 
Er + dE, per unit time: 


dN. T dN. 
=" =] dE, | dr’ q (Ep, 7) —"(E 
te = pcg = | Ee f i aol pet) SE (Ex) (50) 
where: 
qp(Ep, T") = exp(—7’) dN, /dEp. (51) 


Here dN,/dE, is the wnattenuated proton production rate at the source, in units 
GeV~'s~!; it is often labeled as jp(Ep) in the literature. The proton rate is absorbed 
on a target with total optical depth 7 producing secondary pions with an energy 
distribution dN,,/dE,, normalized to unity; it can be parametrized in terms of the 
multiplicity of ny and the average energy < Er > of the pions produced: 
ae = nyt 0 (Er — (Er)) , (52) 
The transmittance T of the target is 


l 
T= exp(- | dr'a (r")), (53) 


where / is the path length of the beam in a target and a is the attenuation coefficient 
(the inverse of the mean-free path A) which is determined by the cross section and 
the density n 


a=no, (54) 
and, assuming isotropy, 


T=nlo =No, (55) 
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where N is the column density. 
Typically, one also makes the approximation that the proton-proton cross section 
is independent of energy, o ~ 3- 107° cm? and therefore the integrals over T and 


Ep separate: 


tgs = (1 - exp(=7)) f aB, T nrs (Er Siis (56) 


In most astrophysical situations, where the optical depth is small, one can replace 
the production efficiency (1 — exp(—7)) > T, with 7T = lno. 

The integral can be performed by rewriting the delta function as a function of 
Ep. Each time a proton interacts, it deposits KE, energy into (nr) pions of average 
energy (Er); here, « is the total proton inelasticity previously introduced. Energy 
conservation implies that 


kEp Se Sate (Er) . (57) 
As before, x, is the average relative energy of a pion relative to the initial proton 
energy 
_ (Ex) O K 
C= a a (58) 
We obtain the result that 
dN, 
Int = Te i aa! (Er — £:Ep) , (59) 
or 
1 dN, (E 
„t = Non, — (| =). 60 
PU oe AEn (=) (60) 


The result is transparent: the pion source function is proportional to the intensity 
of the proton beam, the optical depth of the target and the cross section, and the 
multiplicity of the pions produced. 

We are now able to compute the rate at which neutrinos are produced, with the 
total rate given by summing the emissivities of the muon neutrino from the pion 
decay, and those of the second muon neutrino and the electron neutrino from the 
muon decay: 


dv,tot = af) + qa?) T ave + (61) 


Assuming that the total energy of the pions is distributed equally among the four 
decay leptons, 


dE, 
qu; (En) = Qn (4E) Tr. 4qr (4EL,) (62) 
dE; 
for each neutrino, vi = Ve OY Ve, Vu, Vu- Because IceCube does not distinguish 


neutrinos and antineutrinos, we will not separate them. In this approximation the 
number of neutrinos produced per flavor in the energy bin dE, originate from the 
original pion in the energy bin dE, = 4dE,; therefore, q (Er) dE, = qr (4 Ep )dEz. 
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We next have to translate the results for neutrino production at the source to 
the point source flux Eyy, (Ep) observed at Earth, where a labels the neutrino 
flavor. For a single source the emissivity at Earth is q,(E,)/4mr?, where r is the 
distance to the source. So far IceCube has not pinpointed such a flux, instead it 
discovered a diffuse flux from a, yet unidentified, source population with uniform 
density p(r) in the Universe: 


gv, (Ev) =z ferot)? Ey) (63) 


a ? 
where the first factor 1/47 is introduced in order to define the diffuse flux with the 
conventional units GeV~!em-?s~'sr—!. Therefore, 
1 V E, 
$v, (Ev) = — J irda). (64) 
j An Arr? 


or 


va (Ev) == f ero) )qu (E). (65) 


For illustration we assumed a a Universe, an approximation that is approx- 
imately valid for nearby sources. Integrating over the cosmology of the Universe is 
done by changing the integration from dr > cdt > cdz(dt/dz), with dz/dt = H(z), 
the Hubble scaling factor. We thus can rewrite the integral in covariant form 

c dz 


bv, (Ev) = r] TE 


For the standard ACDM cosmological model, the Hubble parameter scales as 
H? (z) = HÊN + 2)8Qm + Qal, with On œ 0.3, Qa œ 0.7, and the Hubble dis- 
tance c/Ho ~ 4.4 Gpc.!?1 

In the following, we will typically assume that the neutrino emission rate qv, 


p (2) q (1+ 2) By) . (66) 


follows a power law E~7 and the flavor-averaged neutrino energy density can then 
be written as 


p D Pin (E) = Ge pp Pg D Pio B (67) 


where we introduce the redshift factor 
E= E (+2)? p(z) 
Va + (F zm (0) ” 

A spectral index of y ~ 2.0 and no source evolution, p(z) = po, yields €, ~ 0.6, 
whereas the same spectral index and source evolution following the star formation 
rate yields p ~ 2.4. An earlier discussion on how to introduce the cosmological 
evolution of the sources can be found in section 5.2. 

The identical procedure can be followed for the production of neutral pions 
that subsequently decay into two gamma rays. This will lead to Eq. 8, previously 
introduced: 


(68) 


Krp Ny 


52E Bt aoe lama G 
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This is a powerful relation because it only depends on the ratio of neutral to charged 
pions produced which is determined by isospin. 

We complete this section with a discussion of the multimessenger relations be- 
tween neutrinos, photons and protons. Eq. 60 can be rewritten as 


Ne Ex 
lrt = T dp ( ) , (70) 


Tr Tr 
or, 
2 Ky 2 
En n+ (Ex) = fr ee [Es ap(Ep)| Ep=Ex/tx ’ (71) 


with, as before, Kr = 2 for pp and K, = 1 for py interactions. The charged pion 
production rate q,+ is proportional to the rate of the protons accelerated in the 
cosmic accelerator, qp, by the “bolometric” proportionality factor xf, previously 
introduced, and « is the inelasticity factor of 0.5(0.2) for pp(py) interactions. In 
general, the “proton” emission rate, qp, has to be generalized to the composition 


of the cosmic rays that may not all be protons. Using superposition one can relate 
the spectra of nuclei with mass number A as gy(En) = } 4 A*qa(AEn). 

As previously pointed out, the relation between the photons and neutrinos do 
not depend on the initial proton beam. With the usual approximations 


E, 1 1 E, 1 1 


T= = n= = = Ta ? 2 
bep gn ay ae an D 
the neutrino production rate q,, can be related to the one for charged pions: 

1 

3 Sees (E,) > [En Qn+ (Ex )lg wap, £ (73) 


Using Eqs. ws and (73), we arrive at the final relation for neutrino production: 


Q 2 Etv (E h i+ K, [Ez (Ep)] Ep=4E,/Er ` (74) 


The production rate of gamma rays from the decay of neutral pions can be obtained 


in exactly the same way. 

From the two equations for the productions of neutrinos and gamma rays, one 
can eliminate qp to obtain a model-independent relation that only depends on the 
relative contribution of charged-to-neutral pions, 


Er 
7 ə Edva (E xI [E+a(E,)] Ey=2E, ` (75) 


Here, the prefactor 1/4 accounts for the energy ratio E,/E, ~ 1/2 and the two 
gamma rays produced in the neutral pion decay. The relation simply reflects the 
fact that a 7° produces two y rays for every charged pion producing a v, +, pair, 
which cannot be separated by current experiments. This is the simple counting that 
was used to derive the energy density in gamma rays in the universe accompanying 
the flux of cosmic neutrinos observed by IceCube, shown in Fig. 22. 
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7.1. Active Galaxies: A Worked Example 


An active galaxy presents multiple opportunities for the acceleration of cosmic rays 
in the inflows and outflows or jets associated with the supermassive black hole. The 
high-energy cosmic rays will subsequently produce neutrinos in interactions with a 
variety of targets such as the hydrogen and molecular clouds in the galactic disk. 
It is therefore useful to start by considering a generic beam dump where a beam 
of protons with an initial rate dN,/dE, interacts with a target of density n over a 
distance l using the general formalism introduced above. For analytic calculations 
one can use a simple parameterization of the average pion multiplicity that allows 
for its observed increase with proton energy above the threshold pp => pp+ata7:1 


E, — Ep \ 1^ 
Nat =2/ a , (76) 


Also the increase of their average energy can be adequately parametrized as: 


1 
(Ex) = 5 (Ep - mpe? )?/4 GeV. (77) 
Given the number of protons produced by the accelerator per energy and time 
interval 
dN, E,- Wipe y ” 
— = An) 2 78 
dE, ° ( GeV ) it) 
we obtain a pion rate at the source using Eq. 56: 
—3(7-3) 
6E, 3 2 
art (Er) y 26nl Apo (=) . (79) 


The final result for the total neutrino emission rate at the source is given by 


24E,\ 781+ 
GeV ) ? 
which provides us with an estimate of the total neutrino flux at the source in terms 
of three key quantities: Ap and y, the normalization and spectral slope of the 
luminosity of the accelerator, and the column density of the target N ~ In. The 
spectral index is routinely taken to be y = 2, a value suggested by diffusive shock 
acceleration and, in any case, typical for the spectra of gamma rays observed for 


qv,tot z3x 102 nl Ap Oo ( (80) 


nonthermal sources. 

The column density N = In is taken from astronomical information, with l being 
the distance that the cosmic rays travel through a target of density n. For instance 
l could be the diffusion length of the cosmic rays propagating in the magnetic field 
of the galaxy. Diffusion allows the protons to interact with the typical density of 
hydrogen of n ~ 1cm~® over an extended path length determined by the diffusion 
time. The diffusion distance before escaping galaxy is given by 


dai f =2 y D(Ep) tesc , (81) 
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where D is the diffusion coefficient and tese is the escape time. In this case, the 
size of the target is identified with l = ctese. For instance l ~ 10 kpc for our own 
Galaxy. 

An interesting exercise can be performed by fixing the proton luminosity A, from 
the requirement that the aggregate diffuse proton flux from all AGN reproduces the 
cosmic ray flux observed at Earth; see, for instance, reference 41. Assuming all other 
parameters entering in the calculation to be similar to those observed in our own 
Galaxy, one obtains a diffuse neutrino flux that is consistent with the one observed 
by IceCube. The accompanying flux of photons also accommodates the diffuse high- 
energy photon flux observed by Fermi. This attractive solution to the cosmic ray 
problem is unfortunately not supported by other aspects of the IceCube data, most 
prominently the lack of correlation of the arrival directions of high-energy neutrinos 
with Fermi sources.!?? Also, this attractive solution to the cosmic ray problem 
implies that radio galaxies are responsible for the majority of Fermi photons, which 
may be an issue with an extragalactic diffuse flux apparently dominated by blazars. 

An alternative exercise is to determine A, from the radio luminosity of AGNs, L, 
that results from the synchrotron radiation of electrons accelerated along with the 
protons, with Le = KL and « > 1. A, is obtained from the assumptions that protons 
and electrons are connected by a constant fraction fe, with Le = feLp, which is on 
the order of 0.1.*° With these assumptions we obtain a relation between the total 
cosmic ray and radio luminosities: 


dN, xL 
Ly = / — dE, & =, 82 
7 dE, p fe (82) 
and 
— _ xL —1 —2 
Ap = A,(L, 2) = fe (In (Emax/ Emin )| GeV (83) 


for y = 2; generalization for y Æ 2 is straightforward. Observed values of x can be 
found in reference.*° 

Having related the proton flux to the radio luminosity, we obtain the neutrino 
for a single AGN from Eq. 80. The diffuse flux that can be confronted with the Ice- 


Cube data is obtained using the formalism for summing over the sources previously 


introduced: 
dv,tot dn aan dV 
p, = 2 dzdL. 4 
Ji ndi) dV dL dz” (84) 


Here, dz is the luminosity distance, dnagn/dV dL is the radio luminosity function 
of the AGN, and dV/dz is the comoving volume at a fixed redshift z. The radio 
luminosity function is usually separated into the product of a luminosity-dependent 
and a redshift-dependent function, dnagn/dV dL = g(L)f(z). The result matches 
the IceCube observations; more details can be found in reference 40 for column 
densities typical for the relatively dense targets near the AGN cores, as had been 
suggested by reference 171. 
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The two calculations illustrate how two different mechanisms manage to accom- 
modate the IceCube result. Hooper*! normalizes the proton flux to the cosmic ray 
flux and generates the neutrino flux by diffusing the protons confined in the Galaxy, 
while Tjus et al.*° relate the proton flux to the radio emission of the galaxies and 
produce the neutrinos in the dense matter near the supermassive black hole. 

High-energy neutrino emission from photohadronic (py) interactions, where the 
accelerated particles interact with dense radiation fields, can be modeled in a similar 
way. Following the scheme previously introduced, the optical depth of py interaction 
is given by 

Clesc 


(85) 


Ta = 
PY 
Apy 


where Ctesc is the distance traveled by the accelerated protons in the photon 
target, and Apy is the mean free path in that target 
dn. 
-1 
dota f TE Corn) (86) 
Here, dn/de, is the spectrum of the target photons and a;,, is the cross section 
for py interaction. The square of the center of mass energy s in the interaction is 
given by 
s= m? +2mpe, (87) 
where ¢’ is the photon energy in proton’s rest frame. The target photon spec- 
trum can be approximated as a power law, or it may be a monochromatic “line”, 
approximated by a delta function at a specific energy. By considering the photo- 
meson interaction kinematics, one can find the neutrino flux for a photohadronic 
beam dump by taking in to account the initial cosmic ray energy, the optical depth, 
and the kinematics of the photoproduction of pions: 


dN, dN, 
iE, = | dE TEE Er Belton (88) 
where £(Ep, Ep) encapsulates the energy distribution of the final state pions: 
do dN, 
E, E) = | dE," (Ex; Ep) — (4 Er). 89 
(Ep E) = | dE, E (Er Ep) E (Er: En) (89) 


dN,,/dE, is the energy distribution of neutrinos produced in the decay of charged 
pions with energy Er. 

As was done for the pp case, Eq. 50, the initial proton spectrum can be modeled 
as a power law. The integral on s in Eq. 89 can be simplified when particles energies 
meet the required threshold for the A-resonance. We will use this approximation 
later when exploring the neutrino emission from blazars. For more details on the 
kinematics of the py interaction, see references 172 and 173. 

In reality, AGNs are episodic sources producing gamma rays in bursts, with 
episodes where the flux increases by orders of magnitude for periods from seconds 
to months, sometimes years. A correlation of the arrival of IceCube neutrinos in 
coincidence with such bursts can provide a smoking gun for their origin as illustrated 
by the observation of TXS 0506+056. We will this possibility next. 
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7.2. Neutrino Flux Associated with the Plasma Blobs in the Jets of 
Active Galaxies 


Particle accelerators provide a textbook example of special relativity. So do cosmic 
accelerators, but they can do much better. Instead of accelerating protons, they 
can boost really large masses to relativistic velocities. The radio emission reveals 
that the plasma in the jets of active galaxies flows with velocities of 0.99c. A 
fraction of a solar mass per year can be accelerated to relativistic Lorentz factors 
of order 10 leading to luminosities close to the Eddington limit. In the collapse 
of very massive stars 10°! ~ 105? erg/s is released in a fireball that expands with 
velocities of 0.999 c. It is therefore more convenient to introduce the Lorentz factor 
T = [1 — (4)?]-1/?. The plasma flow in the jets can be efficiently extracted from 


Fig. 26. A blob of plasma moves with a velocity v during a time t at an angle 0 relative to the 
line of sight of the observer. 


the black hole in terms of “blobs” of plasma whose radiation is clearly visible in 
radio and X-ray maps. When observing relativistic bulk motion moving close to 
the speed of light, the different paths of the photons reaching us must be taken into 
account. Fig. 26 shows a blob B moving with a velocity v at an angle @ relative 
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to our line of sight to the black hole A positioned at the base of the jet. After 
a time At the relative separation of the the two photons on the sky is vsin and 
the photon from A is delayed by a time (c — vcos@)/c. This leads to an apparent 
velocity of the photon from point A: 


v sind 


USTe ap (90) 


For small viewing angles the apparent velocity exceeds the speed of light when 
v—c. This is a purely geometric effect arising from the Doppler contraction of the 
photons. The maximal superluminal velocity Vapp = I'v is reached for 6 = Tv/c. 
Some important implications: 


e Ifthe blob A radiates isotropically, the radiation will be strongly anisotropic 
in the frame of the observer with one half the photons observed in a semi- 
aperture sinf = 1/T. 

e The difference At between the photons will be contracted in the frame of 
the observer. After a travel time t, of the blob in Fig. 26 the photons are 
separated by a time 


At = ta (1 — v cosb) = Tt, (1 — v cosb) > Tt, (91) 


with the last equality valid for the maximal superluminal effect. The factor 
[It (1 — vcos@)]~* is referred to as the boost factor 6. 

e The energies of the photons are boosted because frequency v behaves as 
the inverse of time. 

e The photon luminosity is actually boosted by a factor r? at maximum. 
This follows from the fact that the number of photons has to be conserved 
between the two frames. Origin of the T? dependence of the luminosity fol- 
lows from the frame invariance of the number of particles in a relativistically 
invariant volume d°ad°p; see the textbooks by Rybicky and Lightman!”* 
and Dermer,!”° for instance. 


Overall what is happening in the blazar is not quite as spectacular as what we see, 
observing boosted energies and luminosities over contracted times relative what is 
actually the reality in the frame of the blazar. 

Blazars are the brightest sources of high-energy gamma rays in the Universe 
making up most of the diffuse extragalactic gamma ray flux observed by the Fermi 
satellite. Their energy density E24, is very similar to the one measured by Ice- 
Cube in cosmic neutrinos, suggesting a possible common origin. Their engines 
must not only be powerful, but also extremely compact because their luminosities 
are observed to flare by over an order of magnitude over time periods that are occa- 
sionally as short as minutes. Blazars are AGN with the jet directly pointing at our 
telescopes thus further boosting the energy of the gamma rays and contracting the 
duration of the episodes of emission. The drawings of a blazar, shown in Figs. 7 and 
27, display its most prominent features: an accretion disk of stars and gas falling 
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onto the spinning supermassive black hole as well as a pair of jets aligned with the 
rotation axis. Large magnetic fields originating from the inflow of particles on the 
black hole are wound up along the rotation axis and launch a pair of jets that are 
the site of the acceleration of electron and proton beams as conceptualized in Fig. 7. 


oe 


black accretion disk 
hole 


Fig. 27. Blueprint for the production of high-energy photons and neutrinos near the super-massive 
black hole powering an AGN. Particles, electrons and protons, accelerated in sheets or blobs moving 
along the jet, interact with photons radiated by the accretion disk or produced by the interaction 
of the accelerated particles with magnetic fields. The jet is pointing at the Earth in the subclass 
of AGN dubbed blazars. 


The energy spectrum of gamma rays radiated by a blazar jet exhibits the classic 
double-hump structure with lower energy gamma rays originating from synchrotron 
radiation by a beam of accelerated electrons, and a higher energy component result- 
ing from inverse Compton scattering of the same photons by the electrons. Other 
sources of photons may also be upscattered to high energy. In Fig. 28 we show 
the energy spectrum of TXS 0506+056 measured in the multimessenger campaign 
following the neutrino IC170922. What is special about this spectrum is that the 
production of neutrinos reveals the co-acceleration of protons along with electrons. 
While the relative merits of the electron and proton blazar have been hotly de- 
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bated, the opportunities for studying the first identified cosmic ray accelerator are 
wonderfully obvious. 
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Fig. 28. Multiwavelength spectrum of TXS 05060+056 from the multimessenger observations 
following the IceCube neutrino alert IC170922!7° 


Because of their reduced synchrotron radiation, protons, unlike electrons, effi- 
ciently transfer energy in the presence of the magnetic field in the jet. They provide 
a mechanism for the energy transfer from the central engine over distances of par- 
secs as well as for the observed heating of the dusty disk over distances of several 
hundred parsecs. When protons are accelerated along with the electrons in blobs 
of plasma opportunities exist for the production of neutrinos by pion photoproduc- 
tion, e.g., in interactions of the plasma with photons radiated off the accretion disk 
and the dusty torus in the AGN; see Fig. 27. Alternatively, the neutrinos may be 
produced inside the plasma in interactions of protons with photons radiated by the 
co-accelerated electrons; see Fig. 29. 

Confronted with the challenge of explaining a ~ E~? high-energy photon emis- 
sion spectrum reaching TeV energies, often radiated in bursts of a duration of less 
than one day, models have converged on the blazar blueprint shown in Fig. 29. 
Particles are accelerated by shocks in blobs of matter traveling along the jet with a 
bulk Lorentz factor of IT ~ 10 and higher. The duration of the burst is associated 
with the size of the blob. The T factor combines the effects of special relativity and 
the geometry of the moving source; we will routinely assume that maximal Doppler 
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Fig. 29. A “blob” of accelerated protons and electrons moves towards the observer with a boost 
factor I. Neutrinos are produced in the interactions of protons with synchrotron photons radiated 
by the electrons. [Note that this is not a Feynman diagram, the photon is real.] 


factor because the value of I is often treated as a free parameter in any case. In 
the following, primes will refer to a reference frame attached to the blob, which is 
moving with a Doppler factor I relative to the observer. In general, the transfor- 
mation between blob and observer frame is R’ = TR and F’ = E/T for distances 
and energies, respectively. Blazar bursts appear more spectacular to the observer 
than they actually are because, in the frame of the blob, distances are larger and 
times longer. 

Episodic high-energy emission is associated with the periodic formation of these 
blobs. The blobs are clearly identified in X-ray images of AGN jets, e.g., in M87. In 
order to accommodate bursts lasting a day, or less, in the observer’s frame, the size 
of the blob must be of order T c At ~ 107? parsecs or less. The blobs are actually 
more like sheets, thinner than the jet’s size of order parsec. The observed radiation 
at all wavelengths is produced by the interaction of the accelerated electrons and 
protons in the blob with the multiple radiation fields in the complex AGN structure, 
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for instance synchrotron photons produced by electrons or the ambient radiation in 
the AGN which often has a significant prominent component concentrated in the 
so-called “UV-bump” of ~ 10 eV photons. 

In order for a proton accelerated in the jet to produce pions on target photons 
of energy Ey, a process dominated by the photoproduction of the A resonance, it 
must exceed the threshold energy: 


m2 om m? 
Poia R (92) 
p 4E! 
or, 
2 2 
MA —m 
Ep > P, (93) 
4 4E, 


in the observer’s frame. Furthermore, the blob must be transparent to gamma rays 
that are actually detected at Earth from the burst. Therefore the center-of-mass 
energy s must be below the threshold for y + y interactions in the blob, or 


s < (Me), (94) 
or, 
E- Ez obs < MZ. (95) 


The kinematic constraints imply that T must exceed a few (~ 30) for target 
photon energies 10eV (1 KeV) typical for UV (synchrotron) photons. These values 
of T render the jet transparent to 50 GeV (500 GeV) photons, for example. 

From the observed luminosity L, we deduce the energy density of target photons 
in the shocked region of size R’: 


p te Re a E 
Ppa ÍTR” SOP fn(TrcAt)3 ~ Ane? T4At2 ’ 


u (96) 
where At is the duration of the flare. Here we just boosted energies and contracted 
times by a factor I from the primed blob frame to the observer frame. 

The fraction of energy fr lost by protons to pion production when traveling a 
distance R’ through a photon field of density n}, = u,/E‘ is given by the number 
of energy loss lengths of the proton 


E 
i oat PAY 97 
Oey “dE y/dt ve 
that fit inside the photon target of size R’, or 
R' R' 
h= sk (98) 
Cipy Apy 


where X,,. is the proton interaction length, CpyA—nn+ ~ 10778 cm? and the inelas- 
ticity x ~ 0.2 the energy transferred from the proton to the pion in each interaction. 
To be clear, we are not asking how many times the proton interacts inside the target 
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region, but after how many interactions it has lost its energy and no longer produces 
pions. The final result is given by 

T R' L 1  3koA 
fr = Acer Eph (eA) 4re ` 
If f, approaches unity, pions will be absorbed before decaying into neutrinos, re- 
quiring the substitution of fy by 1 — e7/. 

For a total injection rate in high-energy protons E%,Qwn, the total energy in 
neutrinos is 1/2f;t47 ER} Qn, where ty~10!° Gyr is the Hubble time. The secondary 
v, have energy E, = x2, Ey, with x, œ 0.05, the fraction of energy transferred, 
on average, from the proton to the neutrino via the A-resonance. The factor 1/2 
accounts for the fact that 1/2 of the energy in charged pions is transferred to v,,+0,. 
The neutrino flux is 


(99) 


dN c 1 1 


Z T A N (20) 2 
P, = Bu TE, An E, la fm ty EnQwn], (100) 


We presented arguments, centered on the 2014-15 neutrino flare of TXS 
0506+056 with a duration At ~ 110 days, that the observed cosmic neutrino flux 
required an efficiency for producing pions, f,, close to unity. The requirement can 
be accommodated when protons in neutrino-producing jets interact with the 10eV 
photons in the galaxy via the A-resonance; py > A — aN. The pion efficiency 
of the jet depends on the Lorentz factor of the jet (T), the target photon energy 
(E), luminosity of the target photons (L,), and the duration of the flare (At), for 
instance, 

1 3KOA 
E,,T?(cAt) 4re ’ 


ee ee Ly 10eV 1 110d 
© =**\4=6) x 10 erg/s)\ E, /\T2)\ At 
3KOA 
—> 102 
x ( Arc? ) cu?) 
This suggests that sources with small Lorentz factors and a UV luminosity ex- 


ceeding O(104°) are required for the production of high-energy cosmic neutrinos. 
177 


(101) 


and 


This level of UV luminosity has been reported in reference. 

The above estimate was performed assuming a monochromatic target photon 
spectrum, i.e., a line spectrum. A more realistic scenario would take into account 
the energy distribution of the target photon, dn./de,, introduced in Eq. 86. In this 


case,!"® 


co co 

tz (Ep) = aa | de’ Tlel f dee *dn.,/dey, (103) 
pve, e'/2Tp 

where ¢' is the photon energy in the rest frame of proton, I’, is the proton Lorentz 


factor in the comoving frame, Kp is the proton inelasticity, and eip = 145 MeV is 
the threshold photon energy for photomeson production. 
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The target photon can be modeled as a power-law, x €~°7 with a normalization 
that is fixed bolometrically by the photon luminosity of the source. Equation 103 
can be simplified when introducing the A-resonance approximation: 


2koa AEA 3L, Ep ae 
E; 


tesc 
fpr (Ep) © > 


~ , 104 
tny lta, ča 4ncRI?E, (104 


where ča ~ 0.34 GeV is the energy value where the cross section peaks, Aéy ~ 
0.2 GeV is the width of the resonance peak, and E? ~ 0.5 T? mpc? Ea /Ey. One can 
verify that the pion efficiency is equal to the one obtained by Eq. 101 when a, = 1. 


7.3. An Alternative Approach for the Modeling of Neutrino Pro- 
duction by Proton Beams on External Photons 


There is an alternative approach for modeling the production of neutrinos when 
a jet of protons accelerated near the black hole interacts with dense photon fields 
surrounding it. In this model the duration of the emission reflects the size of the 
radiation field rather than that of the blob; see Fig. 30. There are multiple oppor- 
tunities for the protons to find dense photon targets, from photons associated with 
pair plasma to the broad-line photons. The framework is illustrated in Fig. 30. The 
key feature of this model is that in the frame of the protons the target photons 
acquire boosted energies to produce neutrinos. 


Broad Line 
Region 


accretion disk 


Fig. 30. A jet of accelerated protons and electrons move towards the observer with a boost factor 
T. Neutrinos are produced in the interactions of protons with dense photons fields in the galaxy, 
e.g., photons in the broad line region close to the central black hole. 


The calculation is similar, but in this case a continuous jet of protons interacts 
with a shell of photons 7R’7AR’. As before, the opacity of the photon target is 
given by: 

AR’ 


PY 


l-e =k =K opn AR. (105) 


February 3, 2022 1:35 ws-rv961x669 Neutrino Physics and Astrophysics main page 68 


68 IceCube Neutrinos 


Here Tey is the density of photons in the jet frame: 


1 dN L' At’ 
rE) = ag | E er VE (106) 
or, with V’ = 4r R? AR’ and R' = cAf', 
z 1 L; 
Iaei RO ang E aog 
In the observer frame: 
1 1 1 L 
e "ny Tre I? cAt E (108) 


In the last transformation the luminosity transforms as T?; this was discussed earlier 
in the section when we introduced superluminal motion. The observed luminosity 
of the Galaxy L, should be reduced to only include the fraction fgrr associated 
with the broad line region. Also cAt is identified with Rgrr, the distance over 
which the jet interacts with broad line photons. Therefore: 


_ 1 1 1 ferry 
~ “rv Inc I? Reir E ` 

Guided by the 2014-15 burst of TXS 0506+056, one may ask the question at 
what point the source becomes opaque to photons but, at the same time, provides 
a sufficiently dense photon target to produce neutrinos. Identifying Rgrr with the 
110 day duration of the burst yeilds the reasonable requirement that fearrl, > 
104’erg/s for 10eV photon energy. Also, the value of Rgrpr is consistent with the 
common parametrization!”? that 


Regier = 10" cmy Las, (110) 


where L45 is that L in units of 10*°erg/s. 

We note that the same phenomenology can be reproduced with the blob model 
(Eq. 101) with a luminosity that is reduced by a factor fgrr/3 relative to the 
luminosity in the broad line model. 


l-e 


(109) 


However, for the BLR parameters above f < 1, and, in any case the neutrino 
energies too low to accommodate the neutrino energies observed in the 2014-15 
flare. A phenomenology that yields appropriate values of f requires [ > 5.6 and 
LBLR = 10% erg/s. 

The so-called structured jets model represent a variant of the previous model 
where gamma rays and neutrinos are produced when protons with boost factor T 
catch us with slower moving plasma with boost factor I’. In the BLR formalism 
one must substitute I? by TT”. 


7.4. Fireballs and GRB 


While we have confidence that some of the gamma rays produced in our Galaxy 
by the decay of neutral pions, there is no straightforward gamma-ray path to the 
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neutrino flux expected from extragalactic cosmic-ray accelerators as already illus- 
trated by our discussion of active galaxies. We presented model calculations to 
show that AGNs can plausibly accommodate the observations of TXS 0506+056, 
the multimessenger connection between photons and neutrinos is indirect with the 
flare indicating the presence of a powerful jet seen in gamma rays, but with the neu- 
trinos produced in periods where gamma rays are absent. Their energy is spread 
over the electromagnetic spectrum by absorption in the source when the neutrino- 
producing target is present. Although we showed how different blueprints for the 
beam dump can plausibly fit the TXS 05060+056 data as well as the diffuse neu- 
trino flux observed, it is certainly premature to conclude that AGN are the sources 
of the extragalactic cosmic rays. 

As discussed in the introduction, there is attractive alternative. Massive stars 
collapsing to black holes and observed by astronomers as GRBs have the potential 
to accelerate protons up to 100EeV energy. Neutrinos of 100 TeV — PeV energy 
should be produced by pion photoproduction when protons and photons coexist in 
the GRB fireball.'°° The model is promising because the observed cosmic-ray flux 
can be accommodated with the assumption that roughly equal energy is shared by 
electrons, observed as synchrotron photons, protons and neutrinos. This is indeed 
what IceCube has observed. 

The phenomenology that successfully accommodates the astronomical observa- 
tions assumes the creation of a hot fireball of electrons, photons, and protons that 
is initially opaque to radiation. The hot plasma expands by radiation pressure, 
and particles are accelerated to a Lorentz factor I that grows with the expansion 
until the plasma becomes optically thin and produces the GRB display. From this 
point on, the fireball coasts with a Lorentz factor that is constant and depends on 
its baryonic load. The baryonic component carries the bulk of the fireball’s kinetic 
energy. The rapid variations with time observed in the burst spectrum can be suc- 
cessfully associated with successive shocks (shells), of width AR, that develop in 
the expanding fireball. The rapid temporal variation of the gamma-ray burst, ty, is 
on the order of milliseconds and can be interpreted as the collision of shells with a 
varying baryonic load leading to differences in the bulk Lorentz factor. Electrons, 
accelerated by first-order Fermi acceleration, radiate synchrotron gamma rays in 
the strong internal magnetic field and thus produce the spikes observed in the burst 
spectra. 

The interpretation!®! of the IceCube observations of GRB is based on the fireball 
model. The basic assumption is that the energy is dissipated in fireball protons and 
electrons, with the latter observed via their synchrotron radiation. The neutrino 


flux is given by 

dN, € 1 dN. 

2 =(") 5r, |B (E 111 
oft -(2)}. [ee]. am 


syn 


where €p and ce are the energy fractions in the fireball in protons and electrons,'*! 
respectively. The neutrino observations determine the baryon loading in the GRB 
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fireball, €p / €e; no prediction is made. The abundance of protons in the fireball is 
determined by, or, at present, limited by the neutrino observations.’ A prediction 
for the diffuse flux for GRB emerges when assuming that they produce the observed 
cosmic rays. 

Although simulations of GRB fireballs have reached a level of sophistication,!*? 
a simple energy estimate is sufficient to predict the neutrino flux associated with 
GRB fireballs assuming that they are the sources of the extragalactic cosmic rays. 


Ro ~100km cos = 


v 
Cc 


y 


At = 1 msec 


Fig. 31. Kinematics of a relativistically expanding fireball on the left, and the resulting shell(s) 
expanding under radiation pressure on the right. 


The relativistic kinematics summarized in Fig. 31 relates the radius R’ and width 
AR’ of the expanding fireball to the observed duration of the photon burst cAt: 


R! =T°(cAt) (112) 
AR =TcAt (113) 


From the observed GRB luminosity L4, we compute the photon energy density in 
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u! = (L,At) as = Ly (114) 
Y 4&TRP?AR!  4rAt?c37® 


The pion production by shocked protons in this photon field is, as before, calculated 
from the interaction length: 
er E H= E 115 
Na WET NA (z =3 a (115) 
Also as before, 7, is the cross section for py > A — nat. As in the calculation of 
AGN, the fraction of energy going into 7-production is 
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The characteristic photon energy in the problem is 1 MeV, the energy where the 
typical GRB spectrum exhibits a break. The contribution of higher energy photons 
is suppressed by the falling spectrum, and lower energy photons are less efficient 
at producing pions. Given the large uncertainties associated with the astrophysics, 
it is an adequate approximation to neglect the explicit integration over the GRB 
photon spectrum. The proton energy required for the production of pions via the 
A-resonance is 


2 2 
MA = Mp 
E! = Lm (119) 
Therefore, 
r \* (1 MeV 
Ep = 1.4 x 101° = 12 
p= 14x 10 ev (35) ( = ) (120) 
1 
Ey = nT 10" eV. (121) 
We are now ready to calculate the neutrino flux: 
dN, cw. C Up c 1 f1 ; 
= es = „tHE 122 
dE, 4E, 4E, 4E, B p f ma 


where the factor 1/2 accounts for the fact that only 1/2 of the energy in charged 
pions is transferred to v,+v,. As before, E is the injection rate in cosmic rays be- 
yond the ankle (~4 x 1044 erg Mpc™? yr~!) and ty is the Hubble time of ~ 10!° Gyr. 


Numerically, 
dN, 7x 10ev] [ f es 
V = 2 1 —14 —2 —1 —1 T 
7 Ten S E, 0.14| |10 Gyr 


x | L | (123) 


1044 erg Mpc? yr7! 
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The only subtlety here is the r? dependence of the shell radius R’; for a simple 
derivation, see reference 6. The result for the diffuse neutrino flux is insensitive 
to weather the GRBs are beamed or not. Beaming yields more energy per burst, 
but fewer bursts are actually observed. The predicted rate is also insensitive to the 
neutrino energy E, because higher average energy yields fewer vs, but more are 
detected. Both effects are approximately linear. Neutrino telescopes are essentially 
background free for such high-energy events and should be able to identify neutrinos 
at all zenith angles. 

For typical choices of the parameters, r ~ 300 and t, ~ 1077s, about 100 events 
per year are predicted in IceCube, a flux that was already challenged'** by the 
limit on a diffuse flux of cosmic neutrinos obtained with one-half of the instru- 
ment completed in one year.'*4 The energy density of extragalactic cosmic rays of 
~ 1044 TeV Mpc~? yr7! depends on the unknown transition energy between the 
Galactic and extragalactic components of the spectrum. In the end, the predictions 
can be stretched by reducing that number, but having failed by now to observe 
high-energy neutrinos in spatial and temporal coincidence with over 1500 GRB ob- 
servations, IceCube has set a limit that is less than 1% of the PeV cosmic neutrino 
flux that is actually observed. 

One lesson learned from TXS 05060+056 is that sources that emit high-energy 
gamma rays may be cosmic accelerators, however, at the time of neutrino flares the 
source is opaque to gamma rays. The failure of IceCube to observe neutrinos from 
GRB may be the consequence of using GRB alerts based on gamma ray observations. 
Only a few GRB have been identified that are underluminous in gamma rays and 
where the jet “choked” in the dense remnant producing high-energy neutrinos while 
the photons lose energy and are shifted to low energy. Also, GRB may produce 
neutrinos even if they are not the sources of the highest energy cosmic rays. For 
instance, they may produce neutrinos of lower energy in events where the boost 
factor is limited, T < 10. Candidate events have been observed and are referred to 
as “low luminosity GRBs.” 155 

There is also the possibility that high-energy gamma rays and neutrinos are 
produced when the shock expands further into the interstellar medium. This mech- 
anism has been invoked as the origin of the delayed high-energy gamma rays. 
Adapting the previous calculation to the external shock is straightforward.'®° The 
timescale is seconds rather than milliseconds, and the break in the spectrum shifts 
from 1 to 0.1MeV. Although jf, is reduced by two orders of magnitude, in the 
external shocks higher energies can be reached, and this increases the detection 
efficiency. In the end, the observed rates are an order of magnitude smaller than 
in internal shocks, but the inherent ambiguities of the estimates are such that it is 
difficult to establish with confidence their relative neutrino yields. 
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7.5. The “Guaranteed” Cosmogenic Neutrinos 


The production of neutrinos in the sources that accelerate high-energy cosmic rays 
depends on the source and on its environment. In order to efficiently accelerate cos- 
mic rays, any loss mechanism, including pion production in py and pp interactions 
that produce neutrinos, must be suppressed as it reduces the acceleration time. 
Efficient accelerators are likely to be inefficient beam dumps for producing neutri- 
nos. High-efficiency neutrino production can however be achieved by separating the 
sites of acceleration and neutrino production, as in TXS 0506+056. Similarly, af- 
ter acceleration, extragalactic cosmic rays propagate over cosmological distances of 
more than 10 Mpc and can efficiently produce neutrinos on the EBL, predominantly 
microwave photons. 

In this section, we will briefly discuss the production of neutrinos in the interac- 
tions of ultra-high-energy cosmic rays with EBL photons. Soon after the discovery of 
the cosmic microwave background, Greisen, Zatsepin and Kuzmin!?°:!*+ (GZK) re- 
alized that extragalactic cosmic rays are attenuated by interactions with background 
photons. Protons interact resonantly via py ~ At — mtn with background pho- 
tons with mean energy € ~ 0.33 meV at energies Ep ~ (mA — m2)/4e ~ 500 EeV. 
The density and width of the Planck spectrum leads to a significant attenuation 
of proton fluxes after propagation over distances on the order of 150 Mpc above 
an energy EqgzxK ~ 50 EeV. This is known as the GZK suppression of the high 
energy cosmic ray spectrum. Also heavier nuclei are attenuated at a similar en- 
ergy by photodisintegration of the nucleus by CMB photons via the giant dipole 
resonance.!°* 

The pions produced in GZK interactions decay, resulting in a detectable flux of 
cosmogenic neutrinos first estimated by Berezinsky and Zatsepin®® in 1969. This 
guaranteed flux of neutrinos became one of the benchmarks for high-energy neutrino 
astronomy leading early on to the concept of kilometer-scale detectors. The flux of 
cosmogenic neutrinos peaks at EeV energy. It can be calculated from the observed 
cosmic ray spectrum and composition, modulo parameters: 


e the slope of the (power-law) spectrum of the cosmic rays sources that is 
injected in the EBL prior to absorption, 

e the minimum and maximum energy of the spectrum, 

e the cosmological evolution of the sources, i.e., the value of €. 


Ahlers etal.°°:!°° calculated the range of predictions for the GZK neutrino flux by 
varying all parameters within reasonable limits. They also determined the most 
likely flux. The exercise was performed assuming that the highest energy cosmic 
rays are protons which yields the largest neutrino fluxes. 15188189 An observable flux 
is predicted under these assumptions, particularly if the extragalactic component 
of the proton spectrum extends below the ankle. 

An example of models satisfying these assumptions are referred to as “dip mod- 
els,” the ankle results from the absorption of protons by Bethe—Heitler pair produc- 
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Fig. 32. All-flavor upper limit (90% C.L.) from IceCube nine-year analysis (solid black line) 
and IceCube seven-year analysis (dashed line). The limits are derived using a log-likelihood ratio 
method. Differential upper limits for per energy-decade ÆT 1 cosmogenic neutirno flux are also 
shown for Auger,?”2°! and ANITA??? after proper normalization considering the energy bin width 
and neutrino flavor. The predicted cosmogenic neutrino flux by Kotera et al.??? (cyan) and by 
Ahlers et al.t°9 (dark blue) and a cosmic neutrino model by Murase et al.?°* (green) are shown 
for comparison. All model expectations are shown assuming primary protons for cosmic ray 
composition. 


tion on CMB photons. A fit to the observed cosmic-ray spectrum requires relative 
strong source evolution with redshift!°°!%° that further enhances pion production. 
However, the corresponding electromagnetic emission via neutral pions as well as 
e~ pairs is constrained by the isotropic gamma-ray background (IGRB) observed 


by the Fermi-LAT satellite’*+!°* and limits the neutrino intensity of these proton- 
109,195-199 


dominated scenarios. 

Recent upper limits on the cosmogenic neutrino flux have resulted from the 
failure of IceCube to observe EeV neutrinos;®° see Fig. 32. The upper limit on 
the GZK flux can be accommodated by introducing a heavy nuclear composition 
into the spectrum of the highest energy cosmic rays. For nuclei, resonant neu- 
trino production still proceeds via the interaction of individual nucleons in the nu- 
cleus with background photons, but the threshold of the production is increased to 
Ecr 2 AEgzx for nuclei with mass number A. Therefore, efficient cosmogenic neu- 
trino production now requires an injected cosmic-ray flux that extends well above 
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FEoezx. Especially for a heavier nuclear composition of the primary flux, the produc- 
tion of neutrinos on photons of lower energies in the EBL becomes relatively more 
important.9°:!97-203,205-211_ The interaction with optical photons produces neutrino 
fluxes in the PeV energy range, the overall level is however reduced because of the 
lower density of the EBL photons. The PeV neutrino flux observed by IceCube 
cannot be accommodated by the production of neutrinos in the EBL.?!>?13 How- 
ever, it has been recently pointed out®! that the cascading of tau neutrinos of GZK 
origin through the Earth creates the opportunity to contribute events to the PeV 
flux at a subdominant level. 


8. Galactic Sources 


The rationale for building kilometer-scale neutrino detectors is that their sensitivity 
is sufficient to reveal generic cosmic-ray sources with an energy density in neutrinos 
comparable to their energy density in cosmic rays“ and gamma rays. The expecta- 
tion was that this relation would also be satisfied by the sources of Galactic cosmic 
rays," 14 

Galactic sources constitute one of the main components of the high-energy emis- 
sion observed at Earth. The most manifest component of the high-energy radiation 
from the Galaxy is the measured gamma-ray emission from sources in the Galactic 
plane. Today, the observed high-energy gamma-rays from Galactic sources reaches 
100 TeV, indicating that Galactic sources are capable of accelerating particles to 
very high energies. While the cosmic rays spectrum extend to energies far beyond 
these energies, Galactic sources has been understood as the principal contributors 
to the cosmic ray spectrum up to the “knee” in the cosmic ray spectrum. 

The rationale for searching for Galactic sources of high-energy cosmic neutri- 
nos rely upon the search for PeVatrons. Galactic cosmic rays are thought to reach 
energies of at least several PeV, the knee in cosmic ray spectrum. The interaction 
of cosmic rays with energies greater than PeV with the dense environments in the 
Galaxy leads to production of charged and neutral meson that would eventually de- 
cay to high-energy neutrinos, as well as gamma-rays, and should be seen in neutrino 
detectors. 

In the context of multimessenger connection, identification of the source of high- 
energy neutrinos and very high energy cosmic rays can benefit from the information 
brought by other channels of observation. For Galactic source, since they are in our 
neighborhood, the very high energy gamma rays that should accompany the high- 
energy neutrinos would not be affected by the extragalactic background light (EBL) 
and therefore, they can be seen in similar energy ranges that neutrino telescopes 
operate effectively. 

Possible neutrino emission from the Galactic plane appears as a subdominant 
component to the observed cosmic neutrino flux. However, recent studies imply 
that identification of the Galactic source of high-energy neutrinos is likely in near 
future. 
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As discussed earlier, the necessary condition for a source to be considered as 
a particle accelerator up to a certain energy is to satisfy the Hillas criterion.?!° 
That is, the gyroradius of the particle shall not excess the size of the astrophysical 
object. Otherwise, it cannot contain the particle. The other relevant parameters 
determining the maximum energy an accelerator can achieve are the acceleration 
mechanism and energy losses at the source. Galactic sources such as supernova 
remnants can meet Hillas’ condition and reach energies beyond the cosmic ray knee 
if the magnetic filed is amplified compared to the average interstellar medium’s 
average value. Possible Galactic wind termination shocks can also accelerate cosmic 
rays beyond the knee in Milky Way; see reference 216 for details. 


8.1. Galactic Neutrino-Producing Beam Dumps: Supernova Rem- 
nants 


Based on energy considerations, the suggestion that supernova remnants are the 
sources of the cosmic rays in our own Galaxy originated in the 1930s.?!" The 
observed energy density of the cosmic rays in our Galaxy is pp ~ 107-7 erg cm™?. 
Galactic cosmic rays do not exist forever; they diffuse within the microgauss fields of 
our Galaxy and remain trapped for an average containment time of ~ 3 x 10° years, 
depending on their energy. The power needed to maintain a steady energy density 
therefore requires accelerators delivering 10+! erg/s. This happens to be 10% of 
the power produced by supernovae releasing 10°! erg every 30 years (10°! erg cor- 
respond to 1% of the binding energy of a neutron star after 99% is initially lost 
to neutrinos). This coincidence is the basis for the idea that shocks produced by 
supernovae exploding into the interstellar medium are the accelerators of Galactic 
cosmic rays. These sources should be visible to the present generation of gamma 
ray instruments. A generic supernova remnant releasing an energy of W ~ 105° erg 
into the acceleration of cosmic rays will inevitably generate TeV gamma rays by 
interacting with the hydrogen in the Galactic disk. 

The emissivity in these pionic gamma rays can be calculated from neutral pion 
production rate,’? with the substitution qro = tqnt in Eq. 56. We first define the 
production rate Qro, per unit volume and time, by converting dl into cdt in the 


7’ integration in Eq. 56. Within the A resonance approximation with ny = 1 and 


k = 0.8: 
Qro = CN Opp for, Np (Ep) ô (Ero — KEp) , (124) 
or 
1 Ezo 
Qro (Ero) = C Opp PRAG ( = ) : (125) 


The emissivity Qo is proportional to the number density of hydrogen targets and 
the number of cosmic rays which is determined by the spectrum 


Ar dN, 
earls E-— 12 
"= ae fa dE’ (125) 
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with np ~ 4 x 107! (Ep/TeV) ~" cm~. The emissivity of photons is obtained 
from the fact that every neutral pion produces two photons with half its energy, 


therefore 
Q; (Ey) = 2 x 2Qro (2Ey) = Qro, (127) 
where the latter equality holds for an E~? spectrum. Our final result is 
E. 
Qy (Ey) ~ ~ Opp Mer (z > =) (128) 


or, assuming an E~? spectrum, 
Q(> 1 TeV) ~ 107? (2) cm~? s7}, (129) 
Mm 


Notice the transparency of this result. The proportionality factor in Eq. 128 is de- 
termined by particle physics, where App = (Napp) ~" is the proton interaction length 
(Opp ~ 40 mb) in a density n of hydrogen atoms. The corresponding luminosity is 


Lal >1 TeV) 2 Qy a (130) 
PE 
where W/pz is the volume occupied by the supernova remnant; given the ambient 
density pg ~ 10~” erg cm~ of Galactic cosmic rays,* a supernova with energy W ~ 
10°° erg in cosmic rays occupies the volume W/ppg. We here made the approximation 
that the density of particles in the remnant is not very different from the ambient 
energy density. 
We thus predict?!*:?!9 a rate of TeV photons from a supernova remnant at a 
nominal distance d of 
J ANy gp — B ITV) 
E>1Tey dE, And? 


TeV W n d \? 
~ 1071? — 107" . 131 
a j cm?s/ \ 105 erg Genes) 1kpc ga 


As discussed in the introduction, the position of the knee in the cosmic ray 
spectrum indicates that some Galactic sources accelerate cosmic rays to energies of 


several PeV. These PeVatrons therefore produce pionic gamma rays whose spectrum 
can extend to several hundred TeV without cutoff. For such sources the gamma-ray 
flux in the TeV energy range can be parametrized in terms of a spectral slope ay, 
an energy Ecut,y where the accelerator cuts off, and a normalization k4: 


aN, (E) ies E 
ZON S k, ( = Lee a 132 
dE 1 TeV SE Ecut,y ( ) 


The estimate in Eq. 131 indicates that fluxes as large as dN,/dE, ~ 107!?-10714 
(TeV! cm~? s71) can be expected at energies of O (10 TeV). 

Such sources are well within the sensitivity of existing high-energy gamma ray 
detectors such air Cherenkov telescopes and large acceptance ground-based detec- 
tors like Milagro and HAWC. They may have been first revealed by the highest 
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energy all-sky survey in ~ 20 TeV gamma rays using the Milagro detector.??? The 
survey revealed a subset of sources, located within nearby star-forming regions in 
Cygnus and in the vicinity of Galactic latitude | = 40°; some of them cannot be 
associated with known supernova remnants or with nonthermal sources observed at 
other wavelengths. Subsequently, directional air Cherenkov telescopes were pointed 
at three of the sources, revealing them as PeVatron candidates with an approx- 
imate E~? energy spectrum that extends to tens of TeV without evidence for a 
cutoff,??! in contrast with the best studied supernova remnants RX J1713-3946 and 
RX J0852.0-4622 (Vela Junior). 

Speculations emerged that some Milagro sources may be molecular clouds il- 
luminated by the cosmic-ray beam accelerated in young remnants located within 
~ 100pc. Indeed, one expects that multi-PeV cosmic rays are accelerated only 
for the short time period when the shock velocity is high, i.e., towards the end of 
the free expansion of the remnant. The high-energy particles can produce photons 
and neutrinos over much longer periods when they diffuse through the interstellar 
medium to interact with nearby molecular clouds.??? An association of molecular 
clouds and supernova remnants is definitely expected in star-forming regions. 

Ground-based and satellite-borne instruments with improved sensitivity are able 
to conclusively pinpoint supernova remnants as the sources of cosmic-ray accelera- 
tion by identifying accompanying gamma rays of pion origin. The Fermi Large Area 
Telescope has detected the pion-decay feature in the gamma-ray spectra of two su- 
pernova remnants, IC 443 and W44.??° In contrast, GeV gamma-ray data from 
Fermi LAT have challenged the hadronic interpretation of the GeV-TeV radiation 
from one of the best-studied candidates, RX J1713-3946.224 On the basis of new 
observations?* the HESS Collaboration has recently argued that the center of the 
Galaxy hosts a promising PeVatron candidate. Detecting the accompanying neutri- 
nos from supernova remnants or the Galactic center would provide incontrovertible 
evidence for cosmic-ray acceleration. 

Particle physics is sufficient to calculate the neutrino flux from potential Pe- 
Vatrons with the production of a v, + , pair for every two gamma rays seen by 
Milagro. The neutrino flux can also be calculated using the formalism discussed at 
the beginning of this section, with qualitatively the same outcome. Operating the 
complete IceCube detector for several years should yield confirmation that some of 
the Milagro sources produce pionic gamma rays; see Fig. 33 and the next section. 

The quantitative statistics can be summarized as follows, for average values of 
the parameters describing the flux, we find that the IceCube detector has the po- 
tential to confirm sources in the Milagro sky map as sites of cosmic-ray acceleration 
at the 3c level in less than one year and at the 5c level in three years.?!8 We here 
optimistically assumed a E~? spectrum, and that the source extends to 300 TeV, or 
10% of the energy of the cosmic rays near the knee in the spectrum. These results 
agree with previous estimates.??5:??6 There were intrinsic ambiguities in these early 
estimate of an astrophysical nature that may reduce or extend the time required 
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Fig. 33. Simulated sky map of IceCube in Galactic coordinates after five years of operation of 
the completed detector. Two Milagro sources are visible with four neutrino events for MGRO 
J1852+01 and three for MGRO J1908+06 with energy in excess of 40 TeV. These, as well as the 
background events, have been randomly distributed according to the resolution of the detector 
and the size of the sources. 


for a 5ø observation.?!* In particular, the poorly known extended nature of some 
of the Milagro sources and the value of the cutoff represent a challenge for IceCube 
analyses that are optimized for point sources. The absence of any observation of an 
accumulation of high-energy neutrinos in the direction of these sources will seriously 
challenge the concept that gamma-ray telescopes are seeing actual sources of cosmic 
rays. 

These predictions have been stable over the years and have recently been up- 
dated in the context of new gamma-ray information, in particular from the HAWC 
experiment; see reference 221. The predicted fluxes for Galactic sources in the 
Northern Hemisphere are close to IceCube’s current upper limits on a point source 
flux of 10-!* TeV~!cm~?s~!;?7 see also reference 228. 


8.2. Identifying Galactic Sources of High-Energy Neutrinos 


The search for Galactic neutrino sources is centered on the search for PeVatrons, 
sources with the required energetics to produce cosmic rays, at least up to the knee 
in cosmic ray spectrum. Supernova remnants have been the primary candidates. 
Accelerating cosmic rays beyond PeV energies, generic PeVatrons should emit very 
high energy gamma rays whose spectrum extends to several hundred TeV without 
a cut off. 

Cosmic rays interacting with matter or radiation at the source, or while prop- 
agating in the Galactic environment, would result in secondary particles that can 
be used to pinpoint the origin of the cosmic rays. The threshold for the production 
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of pions in hadronuclear interaction leads to a feature in gamma-ray flux known as 
the pion bump, appearing at ~ 70 MeV. Such a feature provides direct evidence 
for hadronic interactions. Interestingly, the Fermi Large Area Telescope (Fermi 
LAT) identified this feature in gamma-ray spectra of supernova remnants IC 443, 
W 44, and W 51.??29 Although identification of the pion bump is an intriguing 
evidence for the presence of hadrons in supernova remnants, the particular sources 
are unfortunately not Pevatron candidates. 

As discussed in the previous section, sources of very high energy can be used 
to identify Pevatron candidates with a neutrino telescope; observation of neutrino 
emission would provide the incontrovertible evidence of a cosmic ray source. The 
first very high energy survey of the Galaxy was performed by the Milagro collabora- 
tion.??° This survey revealed the brightest gamma-ray emitters in the Galactic plane 
and pointed towards possible sites of cosmic ray acceleration. Milagro was a large 
volume water Cherenkov telescope with a wide field of view. The sources identi- 
fied by the Milagro survey were subsequently scrutinized by imaging air Cherenkov 
telescopes (IACTs) such as HESS*, VERITAS®, and MAGIC‘ which provided a 
wealth of data on the gamma-ray emission from the Milky Way. Today, with the 
commissioning of the HAWC! gamma-ray observatory, a next-generation instru- 
ment following Milagro, we have reached an unprecedented sensitivity to very high 
energy gamma ray emitters in the Milky Way. 

The initial Milagro survey of the Galaxy presented the first view of the Galaxy 
at 10 TeV and identified a handful of TeV emitters. Today, there are more than 
100 sources identified with TeV emission in the Galactic plane, among them pulsar 
wind nebula (PWN), supernova remnants (SNR), binaries, molecular clouds, and 
shell SNR. However, the majority of the detected sources are classified as uniden- 
tified. The TeVCat** database provides a listing of all TeV sources, including a 
bibliography of the observations.?°° 


The Milagro survey revealed six sources??? 


as potential Pevatrons because of 
their hard spectrum that extended to very high energies without evidence for an 
energy cut-off. Besides the Crab, three of these sources—MGRO J1908+06, MGRO 
J2019+37, and MGRO J2031+4—stood out. The other two, MGRO J2043+36 and 
MGRO J2032+37, were candidate sources with relatively high level of gamma-ray 
emission. Finally, despite its large flux, MGRO J1852+01 did not reach the 5c 
detection threshold. 

Four of these sources were located in the Cygnus region, a star forming zone 
which contains giant molecular cloud complexes and populous associates of massive 
young stars. It provides the substantial requirements for particle acceleration and 
interaction and, as a result, the production of neutrinos. The other two sources were 


tHigh Energy Stereoscopic System 

8Very Energetic Radiation Imaging Telescope Array System 
Major Atmospheric Gamma-Ray Imaging Cherenkov 
High Altitude Water Cherenkov 

**tevcat.uchicago.edu 
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close to the inner galaxy and were characterized by a large flux at TeV energies. 
They were located at low declination, which is the area in the sky IceCube is most 
sensitive to. Therefore, provided that the emission was hadronic, they had the 
highest likelihood of observation in IceCube; see Fig. 33. 

The early predictions based on the measured fluxes and extensions reported by 
the Milagro collaboration determined that IceCube should be able to see neutrinos 


218,225,231 Follow-up observations 


from these sources after five years of operation. 
by Milagro reported a low-energy cut-off in the spectrum of these sources??? which 
reduced their chance of observation in IceCube.”°? The prospect for observing 
these sources in IceCube data is indeed entangled with the uncertainties in the 
spectrum and the extension of the sources. Sources with a hard spectrum, close to 
E~-?, are easier to separate from the soft spectrum of the atmospheric background. 
Also, the number of atmospheric background events increase for extended sources. 
Understanding the precise extension of the source is necessary to obtain the optimal 
sensitivity.?°4 

When IACT and ARGO-YBJ experiments performed follow up observations of 
these sources discrepancies emerged in the spectral measurements and tensions in 
the determination of their extension. A clear picture is hard to assess because each 
experiment faces systematic limitations. The wide-field water Cherenkov telescopes 
like Milagro are more suited to measure gamma-ray flux from extended sources, 
while IACTs are pointing observatories and are therefore limited in exploring ex- 
tended regions. Their small field of view prevents them from separating signal from 
background in off-source regions. On the other hand, IACTs benefit from good 
energy resolution and are likely to provide more reliable spectral measurements. 
However, their sensitivity declines for energies beyond 10 TeV. Here, we revisit the 
prospects for observation of MGRO J1908+06 as one of the preeminent candidate 
sources; for a detailed discussion on the rest of sources, see reference 221. 


8.2.1. MGRO J1908+06 


Figure 34 compares the high-energy gamma ray spectra from MGRO J1908+06 
reported by the air-shower detectors Milagro??? °°" and ARGO-YBJ,?°* as well as 
by the IACTs. The HESS collaboration measures a flux with a hard spectrum and 
no evidence of a cut-off for energies < 20 TeV. It is systematically below the Milagro 
and ARGO-YBJ data.?°° With superior angular resolution, it may be that HESS 
detects the flux from a point source that cannot be resolved by the Milagro and 
ARGO-YBJ observation. MGRO J1908+06 has been observed by VERITAS?*° 
with a flux that is consistent with the one reported by HESS. Finally, the flux 
recently reported by HAWC points towards a similar normalization.?*" 

The experiments disagree on the extension of the source. ARGO-YBJ finds an 
extension of 0.5°,74° while HESS measures 0.34°.7°° VERITAS reports 0.44°.74° 

MGRO J1908+06 is currently classified as an unidentified source. Fermi-LAT 
observes the pulsar PSR J1907+0602 within the extension reported by the Milagro 
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Fig. 34. Gamma-ray flux of MGRO J1908+06 from different experiments: the purple data points 
are measured by HESS,?*° in red the flux from VERITAS,?4? and in cyan the flux from HAWC.?41 
The blue data points show the original measurements by Milagro,?*°»?°° the orange line and 
the shaded orange region show the best fit and the lo band??’ . The dotted area identifies 
the ARGO-YBJ lo band.?3° The green lines show the spectra obtained for a ay = 2 with a 
normalization fixed to the best fit reported by HESS. The result is shown for 3 cut-off enrgies: 
Ecut,y = 30, 300, and 800 TeV (short-dashed, solid, and long-dashed lines, in green). Figure from 
reference 221. 


collaboration.?*+ On the other hand, the large size and the hard spectrum of TeV 
photons that persist away from the pulsar location, are not characteristic of a PWN 
scenario;7“° MGRO J1908+06 is perhaps a SNR. 

We use the multimessenger relation between gamma-ray and neutrino produc- 
tion rates in conjunction with the parametrization of Eq. 132 to estimate the flux of 
high-energy neutrinos from MGRO J1908+06. One is tempted to associate the cut- 
off with the knee in the cosmic spectrum. There is however some question whether 
supernova remnants can reach several PeV particle, and ~ 300 TeV photon energies. 
On the other hand, it is also believed that the the Galactic cosmic ray spectrum 
extend well beyond the knee. The green lines in Fig. 34 show the result of the 
calculation assuming a, = 2 and fixing the normalization to the best fit reported 
by HESS. We allow the cut-off energy to vary: Eeut,y = 30, 300,and 800 TeV. 
In Fig. 35 we show the corresponding number of events as a function of the en- 
ergy cut-off which we varied between 30 and 800 TeV. Also shown is the number 
of background events for both a point source and an extended source with a size 
of 0.34° extension as reported by HESS. The resulting p-value for observation of 


February 3, 2022 1:35 ws-rv961x669 Neutrino Physics and Astrophysics main page 83 


F. Halzen & A. Kheirandish 83 
MGRO J1908+06 
10! a 
F extended 
z iL point—like 
> E 
= 
= 
2 107! 
10°? 
10° 104 10° 10° 


E, [GeV] 


Fig. 35. Event distribution from MGRO 1908+06 for the spectra parametrized by HESS mea- 
surement.The grey band encodes the uncertainty on the cut-off energy. The black (gold dashed) 
line shows the background from atmospheric neutrinos for extended (point-like) sources. Figure 
from reference 221. 


neutrino emission in the direction of MGRO J1908+06 is shown in Fig. 36. Dif- 
ferent assumptions on the high-energy cut-off of the spectrum and the extension of 
the source result in different likelihoods for the identification of MGRO J1908+06. 
Overall, provided that the gamma-ray flux from MGRO J1908+06 extends to en- 
ergies beyond 100 TeV, we anticipate a 30 observation in 10 to 15 years of IceCube 
data assuming the photon flux reported by HESS. 

For all the reasons mentioned above the MGRO 1908+06 source has been in- 
cluded from the beginning in IceCube’s triggered search for point-like and extended 
sources.°0:227,254,245 A recent search for steady sources in the muon neutrino flux 
with 8 years of IceCube of data?“° found this source as the second “warmest” source 
in the catalog, with a pre-trial p-value of 0.0088. The fitted flux upper limit is con- 
sistent with the expectations discussed above. Figure 37 shows the p-value of the 
neutrino excess in the vicinity of MGRO J1908+06. This search optimistically 
assumed no extension of the source. 

A recent search for PeVatron candidates in the HESS Galactic plane survey 
data?*° establishes a lower bound on the cut-off energy of MGRO J1908-+06 that 
has further supported the idea that the source may be a PeVatron.?*” 

Besides monitoring also MGRO J1908+06 and MGRO 2019+37, IceCube per- 
formed a stacking analysis that searched for the collective neutrino emission from 
the six sources initially identified. The stacking likelihood?“ matches the neutrino 
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Fig. 36. The calculated p-value for MGRO J1908+06 as a function of time for different spectral 
assumptions compatible with the observations reported by HESS;?°° see Fig. 35. The grey band 
encodes the uncertainty resulting from different values of Ecut,y and the source morphology. Green 
lines assume an extended source. Figure from reference 221. 


arrival directions with a collection of sources or a catalog.?*9 
Progress in the search for Galactic sources has been made possible by the com- 
missioning of HAWC as we will discuss further on. 


8.2.2. Pulsar Wind Nebulae 


Progress in instrumentation has revealed more than 150 TeV sources in the Galac- 
tic plane, the majority are classified as pulsar wind nebula. Given that they are 
collectively responsible for a significant fraction of the high-energy emission from 
the Galaxy, their contribution to the neutrino flux may be important. 

Confined inside supernova remnants, pulsar wind nebulae are bright diffuse neb- 
ulae whose emission is powered by pulsar winds emerging from the rapidly spinning 
magnetized pulsars in their center. Relativistic electron-positron pairs are consid- 
ered to be the primary nonthermal component of the pulsar winds. Synchrotron 
emission by relativistic pairs is the dominant process for the low-frequency emis- 
sion in radio, optical, and X-ray, while inverse Compton scattering of synchrotron 
photons is responsible for the TeV emission.? 

The possible role of hadrons in PWN was first suggested as an alternative mech- 
anism for generating the very high energy gamma ray emission from the Crab 
Nebula. Early proposals included the possibility that protons, accelerated in the 


outer gap of the pulsar, interact with the nebula?! or, alternatively, that heavy 
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Fig. 37. P-value map of neutrino excesses near MGRO 1908+06 in point source search with eight 
years of IceCube upgoing muons. Figure from reference 245. 


nuclei accelerated in the pulsar magnetosphere interact with soft photons.?”? More 
recent modeling focuses on particle acceleration in the termination shock followed 
by neutrino production in cosmic ray interactions with the dense environment in 
the nebula, see e.g., references 253—258. Minor contamination of ions at the ter- 
mination shock would result in a significant release of energy in hadrons.?°? This 
is particularly important in the context of what is known as the o problem. The 
parameter o is defined as the ratio of the wind’s Poynting to kinetic energy flux. 
Theoretical modeling of the pulsar magnetosphere and wind point at large ø values. 
However, the magnetohydrodynamic simulations cannot match the shock size and 
expansion speed at same time. This conflict can be resolved if the majority of the 
pulsar wind energy is carried by hadrons and, in addition, results in the efficient 
acceleration of particles in termination shocks.?°° 

Searches for neutrino emission from pulsar wind nebulae have targeted individual 
73,227,254 as well as their cumulative emission in a stacking search.749:?°° So 
far, neither have yielded evidence for a correlation of the arrival directions of cosmic 
neutrinos to pulsar wind nebulae and have thus constrained the hadronic component 
of the termination shock. The results are relevant for evaluating the prospects for 
observing PWN with KM3NeT; many of these sources are located in the Southern 


sources 
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Fig. 38. Differential upper limits on the hadronic component of the gamma-ray emission from 
TeV pulsar wind nebulae. The light grey lines show the observed gamma-ray spectra of the 
sources while the dark grey line presents the sum. Red, orange, magenta and blue steps show the 
differential upper limit on the hadronic gamma-ray emission. The upper limits are obtained by 
converting 90% CL differential upper limit on the neutrino flux, and each color corresponds to a 
given weighting method. Figure from reference 260,261. 


Hemisphere; see reference 258. 

The IceCube stacking analysis presented in?*? targeted a set of 9 supernova 
remnants with associated pulsar wind nebulae that are gamma ray emitters. In a 
recent study, IceCube has examined the correlation of high-energy muon neutrinos 
with a larger set of 35 TeV pulsar wind nebulae with gamma-ray emission above 1 
TeV.76°:26! The selection of the sources were based on HAWC, HESS, MAGIC, and 
VERITAS observations and the associated pulsars listed in the Australia Telescope 
National Facility (ATNF) catalog.7°? In the absence of any significant correlation, 
upper limits were derived on the total diffuse neutrino flux from these sources. 
Figure 8.2.2 shows the results for four distinct hypothesis tests of the multimessenger 
relation between the neutrino and gamma ray fluxes. At the highest energies, 
the stacking analyses find upper limits at the level of the high-energy gamma-ray 
emission and thus start to constrain the role of hadrons in pulsar wind nebulae. We 
will return to this search in Section 8.4. 
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8.2.3. The HAWC Very High Energy Gamma Ray Survey of the Galaxy 


The characteristic energy scale of gamma rays from PeVatrons reaching cosmic 
ray energies beyond the knee is 100 TeV. The ground-based, wide-field telescope 
HAWC gamma-ray observatory has provided an unprecedented view of the Galaxy 
by reaching this energy. HAWC’s survey of the Galactic plane has revealed more 
than 20 new TeV sources, most of them classified as unidentified. 


10° — v, flux derived from 2HWC sources 


= combined y, flux from 2HWC sources 
— 90% C.L. upper limit 
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Fig. 39. Upper limit (90% C.L.) of the stacking search for a flux of muon neutrinos (black) 
associated with sources in the 2HWC catalog (PWN excluded). The thin orange line is the 
projected muon neutrino flux assuming all HAWC photons are produced in hadronic interactions. 
The combined flux (red) shows sum of the individual fluxes. Figure from reference 263. 


In order to find common sources of very high energy gamma rays and neutri- 
nos, the IceCube and HAWC collaborations designed and performed a joint anal- 
ysis.2°° Specifically, the analysis examined the correlation of the arrival directions 
of IceCube muon neutrinos passing through the Earth and sources in the 2HWC 
catalogue,?°* as well as the correlation of the arrival direction of neutrinos with the 
morphology of the very high energy gamma ray emission. By excluding the sources 
identified as PWN, the stacking analysis of 2HWC sources focused on unidentified 
sources. The latter analysis incorporated the full morphology of the gamma-ray 
emission from the plane visible in the Northern Sky. Finally, special regions were 
singled out for study: the Cygnus region, MGRO 1908+06, and the inner Galaxy, 
including MGRO J1852+01 which has been confirmed by HAWC. 

Using eight years of IceCube data, constraints are obtained on the contribution 
of pionic photons to the HAWC gamma-ray flux in the Cygnus region and in the 
section of the Galactic plane visible to both experiments. The upper limits on the 
neutrino flux from the other searches are above the flux derived on the basis of the 
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assumption that all HAWC photons are of hadronic origin. The upper limits from 
the stacking analysis of the 2HWC sources is shown in Fig. 39. More data and 
improved selections will be required. This is also the case for MGRO J1908+06. 

An under-fluctuation of the data for in the Cygnus region led to a stringent 
constraint on the hadronic component of the high-energy emission. 


8.2.4. The Cygnus Complex 


The Cygnus X complex is a nearby star-forming region in our Galaxy. It has been 
extensively studied as a promising nearby source of very high energy cosmic rays 
and neutrinos.??3-76° 2” The complex contains massive molecular clouds, populous 
associates of massive young stars, and luminous HII regions.?"!-?"° The observation 
of a hard gamma-ray spectrum, in combination with dense molecular clouds?"*:?”° 
and a large number of young OB stars?’° supports the idea that Cygnus X is a 
plausible source of cosmic rays and high-energy astrophysical neutrinos. 

High-energy gamma ray emission from this region was tentatively detected by 
EGRET*.?” Later, a hard gamma-ray spectrum was established by HEGRA??°° 
and confirmed by the detection of TeV emission by the Milagro collaboration.?”* 
More recent observations with the MAGIC,?” Fermi,?°:?8° ARGO-YBJ,?*! and 
VERITAS?” instruments have led to a better understanding of the point-like and 
extended emission in this region. 

Because of its extension, the identification of individual sources has been chal- 
lenging for the IACT experiments. In contrast, 4 out of 6 sources initially detected 
by Milagro belonged to this region. MGRO J2031+41 is particularly interesting 
because its location coincides with the Cygnus cocoon. Because of uncertainties 
associated with the origin of the flux from the Cygnus cocoon and y-Cygni, a com- 
plete picture of high-energy emission from MGRO J2031+41 is missing. A combined 
analysis of ARGO-YBJ and Fermi data finds a hard spectrum,”*? which, assuming 
a point source, implies observation by IceCube after ~ 15 years.?71 

An important feature of the Cygnus region is the observation of an extended 
excess of high-energy gamma rays with a hard spectrum referred to as the Cygnus 
Cocoon. It was first detected by the Fermi collaboration?®° and confirmed by the 
observation at TeV energies by Milagro, ARGO-YBJ, VERITAS, and HAWC in a 
region that is partially coincident with the Cygnus Cocoon. The true nature of the 
Cocoon is still a mystery. 

We have described the observation in the Cygnux X region with a model of cos- 
mic ray interactions developed for the central molecular zones of starburst galax- 
ies;?°° see Fig. 40. While at lower energies (Æ ~ 10 MeV) bremsstrahlung by 
electrons dominates the gamma-ray spectrum, the contribution of pionic photons 
grows with energy and dominates at ~ 10-100 GeV. High-energy gamma ray emis- 
sion from inverse Compton scattering is negligible at these energies because of the 


* Energetic Gamma Ray Experiment Telescope 
tHigh-Energy-Gamma-Ray Astronomy 
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Fig. 40. Observed and modeled gamma-ray fluxes from the Cygnus X region. Different compo- 
nents of the spectrum include the gamma-ray emission calculated with the YEGZ model, pulsars 
(PSRs), active galactic nuclei (AGN), supernova remnants (SNRs), unidentified sources (UnID), 
the isotropic gamma-ray background (ISO), and the Cygnus Cocoon. Observational points show 
data from Fermi (black stars),2°° ARGO-YBJ (black circles) ,?°° HAWC (black square),?64 and 
Milagro (black triangle).?7S The vertical dotted black line identifies 300 GeV beyond which the 
spectrum for Fermi 3FGL sources are extrapolated. Figure from reference 284. 


steep spectrum of the cosmic ray electrons. The comparison of the the modeled 
gamma-ray spectrum with the observations from Fermi, Milagro, ARGO-YBJ, and 
HAWC are shown in Fig. 40. The spectral fits provided in the 3FGL Fermi cata- 
logue only cover the region between 100 MeV and 300 GeV. For the Cocoon, we 
have extrapolated the spectral fit to higher energies to compare with TeV-energy 
gamma-ray observations. Finally, we only included the off-pulse emission for the 
three brightest pulsars; see reference 284 for details. 

Assuming that the Cocoon is a single source and that the spectrum is domi- 
nated by gamma-rays from neutral pion decay, the model?** predicts a neutrino 
flux that is slightly larger than the IceCube differential discovery potential at 1 
PeV.??’ Therefore, detecting neutrinos from the Cocoon may be possible provided 
the spectrum extends to PeV energies without steepening. 

We next turn to a region within the Cocoon that surrounds a large molecular 
gas cloud complex located left of Cyg OB2 and centered on (l = 81.5°, b = 0.5°),: 
CygX-North. Its emission is expected to be dominated by hadronic interactions, 
most likely from an unresolved source, either a supernova remnant or a pulsar wind 
nebula. Here the model predicts a neutrino flux below IceCube’s discovery potential. 
A summary of the expected neutrino flux from the cocoon and CygX-North is shown 
in Fig. 41. Clearly, the prospect for identification of neutrino emission from the 
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Fig. 41. The neutrino energy flux predicted for the diffuse gamma ray flux from a YEGZ model 
with a soft cosmic ray spectrum with p = 2.6 — 2.8, the Cygnus Cocoon, and the CygX-North 
molecular cloud complex. Also shown is the IceCube differential discovery potential for a point 
source.??” Note that the IceCube sensitivity to extended sources is lower, and thus this plot 
represents the most optimistic case for detection. Figure from reference 284. 


Cygnus X region critically depends on the slope of the cosmic ray spectrum. 

Arguing that diffusion losses play a significant role in modeling Cygnus X, more 
pessimistic neutrino fluxes are anticipated.?°? This transport model describes the 
broad multiwavelength spectrum but predicts a neutrino flux which approaches the 
sensitivity of IceCube only at very high energies, nominally > 50 TeV. 

Also of interest was the Milagro source MGRO J1852+01 because of its high 
gamma ray flux. However, it did not pass the statistical threshold for detection 
and inclusion in the Milagro survey. HAWC observes a high intensity flux from the 
source region that is mostly centered on the source 2HWC 1857+027. The large 
flux reported by Milagro in a 3° x 3° region most likely originated from this source. 
Given the large flux and the vicinity to the horizon where high-energy neutrinos 
reach IceCube without absorption by the Earth, this region is a sweet spot in 
IceCube sensitivity. Interestingly, the largest excess in the joint IceCube-HAWC 
search is associated with 2HWC 1857+027. 

Overall, the joint IceCube-HAWC analysis did not find any significant correla- 
tion. However, it may have provided clues about the potential sources of high-energy 
neutrino in the Galaxy. There are two main outcomes: first, the hadronic compo- 
nent of the emission from Cygnus region has been constrained. Second, in the inner 
Galaxy the results are consistent with expectations. That is, while IceCube has not 
reached the sensitivity for identifying emission from those Galactic sources, their 
observation is within reach as the data accumulates. 

In a more recent development, the HAWC collaboration reported the measure- 
ment of the gamma-ray flux from the Galaxy above 56 TeV.7°° HAWC’s very high 
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energy survey identifies 9 sources, four of which coincide with Milagro PeVatron 
candidates. Additionally, this survey resolves two more sources above —5° declina- 
tion, where IceCube is most sensitive. 


8.2.5. Binaries 


Microquasars are X-ray binaries of a collapsed object with prominent jets rotat- 
ing around a massive star. They have been identified as potential sites of cosmic 
ray acceleration and neutrino emission.?°” °° High-energy photons are produced 
in the interaction of the jet with the stellar wind and radiation fields associated 
with the high-mass companion star. As usual, an alternative possibility is that the 
high-energy photons result from the inverse Compton scattering of stellar photons 
by relativistic electrons. In hadronic models, protons in the jet interact with cold 
protons in the stellar wind.??!?°? Not only are neutrinos produced, cooling of the 
other secondary particles produced in decay of the charged pions, i.e., muons and 
electrons, alter the broadband emission spectrum of gamma-ray binaries. Studies of 
Cygnus X-1 suggest that protons in the jet could carry a significant fraction of their 
kinetic energy.?°*:?°4 Microquasars were actually suggested as the dominant con- 
tributors to the Galactic component of the high-energy neutrino flux.?°° This was 
partially motivated by the spatial clustering of the first sample of cosmic neutrino 
events near LS 5039, a known microquasar which has been historically suggested as 
a potential cosmic ray accelerator.?°° 

X-ray binaries have periodic outbursts which makes them ideal targets for time- 
dependent searches. These benefit from lower background rates and can identify 
the source with fewer signal events. The periodic emission and X-ray flares have 
been exploited in time-dependent searches for neutrino emission.7°° ?°° 

The discussion presented so far has been centered on potential Galactic sources 
of high-energy neutrinos in the Northern Sky. In the future, with the commissioning 
of KM3NeT the sensitivity to neutrino sources in the southern sky will be signifi- 
cantly improved. The HESS collaboration has detected very high energy gamma ray 
emission near the Galactic center??? with a spectrum that is difficult to account for 
by leptonic modeling. Observation of neutrino emission could establish this source 
as a Galactic PeVatron. The neutrino emission may originate from the Sgr A*, the 
low-luminosity AGN located at the center of the Milky Way.°°°3°! Observations 
of X-ray echos have revealed a high level of activity in its past.°°7-39° Accelerated 
cosmic rays will interact with the very dense environment near the center of the 
Galaxy and generate a flux of high-energy neutrinos. 

RX J1713.7-3946, a known supernova remnant, and Vela X°°* are promising 
sources in the field of view of KM3NeT. It is worth mentioning that the IceCube 
search for point sources using cascades identified RX J1713.7-3946 as the most 
significant source.°?° With cascades IceCube has a reduced sensitivity, by about 
a factor of 5, but covers the Southern Hemisphere. HAWC observes high-energy 
gamma sources in the fraction of the Galactic plane in the Southern Hemisphere 
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that is within its field of view. If hadronic in origin, these sources are within the 
sensitivity of KM3NeT after 10 years of operation.°°° 


8.3. Diffuse Galactic Neutrino Emission 


Over their diffusion time in the Milky Way, very high energy cosmic rays interact 
with the dense environment of the Galactic plane that provides the target material 
for producing neutrinos. Therefore, diffuse gamma-ray and neutrino emission from 
the Galactic plane is considered a guaranteed non-isotropic component in the high- 
energy sky.42307-310 

The intensity of the Galactic diffuse neutrino emission is expected to vary along 
the plane tracing the gas density and the distribution of the sources. The intensity 
and morphology of the diffuse Galactic neutrino flux can be further constrained by 
modeling the observed diffuse Galactic gamma-ray emission profile using the local 
measurements of the cosmic ray spectrum.°®!! 3!% 

The diffuse component of gamma-ray emission from Galactic plane has been 
measured by Fermi-LAT in the range 100 MeV to tens of GeV.314315 A tool rou- 
tinely used to model these observations is the GALPROP®!° code. The diffuse 
Galactic gamma-ray component is derived by assuming a distribution of potential 
sources of cosmic rays and by modeling their propagation and diffusion in the in- 
terstellar medium. While good agreement can be achieved with the gamma-ray 
measurements below 10 GeV, the modeling fails to account for the observed flux 
from the inner Galaxy above that energy.*'* The so-called KRA} models?!73!9 
propose to resolve this by modifying the diffusion coefficient as a function of en- 
ergy. Energy cut-offs at 5 and 50 PeV are chosen to accommodate the cosmic ray 
measurements by KASCADE*”? and KASCADE-Grande.**! Figure 42 shows the 
anticipated morphology of neutrinos predicted in the KRAŽ model integrated over 
energy. 

IceCube has searched for the Galactic component of the diffuse cosmic neutrino 
flux?49:311,322,923 based on a variety of assumptions for the energy dependence of 
the Galactic spectrum, including the one predicted by the KRA, models. 

In a recent study, using ten years of ANTARES muon neutrino and cascade 
data and seven years of IceCube muon neutrino data, a joint search search has been 
performed for a diffuse Galactic component of the cosmic neutrino spectrum.°?+ The 
results of the maximum likelihood analysis finds a non-zero component, however, 
the excess falls short of statistical significance because the upper limits obtained 
from the test fall above the actual Galactic flux; see Fig. 43. 

It has been suggested’? that the extended region of the Fermi bubbles may be 
a site of cosmic ray diffusion and high-energy neutrino emission.*!!°?° Gamma ray 
emission with a hard spectrum up to ~ 100 GeV has been observed.*?" IceCube 
and ANTARES have searched for a diffuse Galactic flux from this region but only 
established upper limits.°°°:°?5 
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Fig. 42. Diffuse gamma-ray emission from the Galactic plane. The color map shows the energy 
integrated neutrino flux, assuming the KRA? model,?!” as a function of their arrival direction in 
equatorial coordinates. Figure from reference 324. 


8.4. Closing In on Galactic Sources of High-Energy Neutrinos 


The arrival directions of high-energy neutrinos indicate an extragalactic origin. 
However, a subdominant contribution from a Galactic component is expected, as 
discussed above. 

Galactic sources are in fact a guaranteed contributor to the flux of cosmic neu- 
trinos. Cosmic rays in the Galaxy are believed reach energies of few PeV. Their 
interaction with dense stellar and interstellar gas and radiation results in the pro- 
duction of neutrinos and gamma rays. 

In this section, we have discussed the prospects for identifying neutrino emission 
from Galactic sources and reviewed the status of past searches. The discussion was 
centered around major gamma-ray emitters in the Galaxy because they provide, if 
not a measurement, upper limits on the production of pionic gamma rays. 

The subdominant Galactic component of the high-energy neutrino flux is esti- 
mated to contribute at most 10% of the total measured cosmic neutrino flux. How- 
ever, identification of Galactic sources with strong neutrino emission is expected 
with upgraded and next-generation detectors. While a significant contribution to 
the flux at very high energies is not expected,*?? features in the cosmic neutrino 
spectrum and the excess flux observed in IceCube below ~ 10TeV may already 
indicate a Galactic contribution from the southern sky with a harder spectrum. 

In addition, according to our best estimates, the searches for point-like and ex- 
tended sources in the Galaxy will soon reach the required sensitivity to be successful. 
In this context, continued improvement of the gamma ray data will be helpful. 


February 3, 2022 1:35 ws-rv961x669 Neutrino Physics and Astrophysics main page 94 


94 IceCube Neutrinos 


KRAy model J 
Combined UL KRAy° : 
- Combined UL KRAy7® ] 
ANTARES UL KRA® 
- IceCube UL KRAy? 


IceCube starting events 


E?d®/dEdQ [GeV cm~?s7! sr7}] 


IceCube up-going vyu 


E [TeV] 


Fig. 43. Combined upper limits at 90% confidence level (blue lines) on the three-flavor neutrino 
flux for the KRA} with the 5 and 50 PeV cutoffs (black lines). The boxes represent the diffuse 
astrophysical neutrino fluxes measured by IceCube using an isotropic flux template with starting 
events (yellow) and upgoing tracks (green). Figure from reference 324. 


HAWC observations of the gamma ray flux at energies exceeding 50 TeV con- 
tinue to support the idea that sources like MGRO J1908+06 may be sources of 
high-energy neutrinos. The projected Southern Wide-field Gamma-ray Observa- 
tory (SWGO)*”® will provide a deeper survey of the Galactic plane in the Southern 
Hemisphere. Soon the Cherenkov Telescope Array (CTA)?%331 and LHAASO??? 
will provide enhanced sensitivity reaching into the PeV regime, essential to identify 
the PeVatrons in the Galaxy. Based on current estimates, CTA is expected to re- 
solve ~ 100 supernova remnants.’ It is also worth mentioning that proposed MeV 
telescopes like eAstrogram and AMEGO have the capability of identifying the pion 
bump in the broadband spectrum of Galactic sources. This will contribute to the 
superior modeling of the diffuse emission because unresolved sources with harder 
intrinsic cosmic ray spectra can have a significant contribution to diffuse gamma-ray 
and neutrino fluxes.°!!:334 

Anticipated improvements in angular resolution and, specifically, in identifying 
more starting muon tracks, will improve IceCube’s sensitivity to Galactic sources in 
the southern sky.°°°:?8° Addition of cascade data sets, despite their reduced angular 
resolution, has proved to add substantial capabilities in this search. 

Future neutrino detectors will considerably improve the likelihood of finding 
Galactic sources of high-energy neutrinos. Development of neutrino detectors in 
the Northern Hemisphere such as KM3NeT and Baikal??? will boost the sensitivity 
towards Galactic center and other regions of interest in the southern sky. The next- 
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generation IceCube detector is expected to improve the sensitivity to sources by at 
least a factor of 5.122:169,338 Increased sensitivity to point-like and extended sources 
will be achieved improved angular resolution, better energy resolution, higher statis- 
tics, and reduction of the background. 

A potential more exotic Galactic component contributing to the high-energy 
neutrino flux is dark matter. The decay and annihilation of dark matter parti- 
cles create distinctive signatures in the observed neutrino spectrum; see e.g., refer- 
ences 339-344. An enhancement in the direction of the Galactic center could result 
from the higher local concentration of dark matter. In contrast, the possible inter- 
action of cosmic neutrinos with dark matter in the Galaxy would create a deficit in 
the neutrino flux in the same direction.**° 

More important, IceCube has a unique capability of identifying the annihilation 
of dark matter particle trapped in the sun. We will discuss this in the next section. 


9. Beyond Astronomy 


We now turn to the study of the neutrinos themselves. Over more than a decade of 
operations IceCube has focused on three traditional methods to study neutrino 
physics: the study of neutrino oscillations exploiting the atmospheric neutrino 
beam, the search for MeV neutrinos from supernovae explosions in our Galaxy 
with their extraordinary potential to reveal neutrino physics, and the search for 
dark matter particles that have accumulated in the sun to annihilate in Standard 
Model particles that decay into neutrinos. 


9.0.1. Atmospheric Neutrino Oscillations 


Deepcore®**® is the IceCube tool to study the oscillations of neutrinos produced in the 
atmosphere. It instrumented a volume with improved photocathode coverage at the 
bottom center of the detector with eight additional strings that were deployed inside 
the regular IceCube grid. The PMTs on these strings have higher quantum efficiency 
and are located in the clearest ice below 2.1km. The DeepCore instrumented volume 
is roughly 10’ mê with a module density about five times that of the surrounding 
IceCube array that is used as a veto removing with high efficiency the cosmic ray 
muon background. 

347-349 use samples of neutrinos of more than 100,000 
events per year from the decay of pions and kaons produced by cosmic ray air 
showers in the Earth’s atmosphere. The initial atmospheric flux is dominated by 
Ve and v, and contains negligible numbers of v,.°°° The atmospheric neutrinos 
interacting in the DeepCore subarray of IceCube travel distances ranging from L = 
20km (vertically downward-going) to L ~ 1.3 x 104km (vertically upward-going 
or the full diameter of the Earth). For vertically upward-going neutrinos, the first 
peak of the v, > vr oscillation probability occurs at 25 GeV. This is comfortably 
above the E, ~ 5 GeV threshold of DeepCore.*°' The energy corresponding to the 
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oscillation maximum is also well above the kinematic threshold of E,, = 3.5 GeV 
for charged current v;-nucleon interactions. 

Results from two separate analyses with different strategies for event selection 
and background estimation are shown in Fig. 44. The main analysis A targets 
a high statistics sample of neutrino events of all flavors resulting in 13462 events 
in three years. The background estimation is simulation-driven. In contrast, the 
confirmatory analysis B is optimized for higher rejection of background events and 
its atmospheric muon background estimation is data-driven. 
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Fig. 44. The 90% allowed region using the data sample from Analysis A in blue compared to 


other experiments.?49:352-355 The best fit point from Analysis A is shown as the blue cross mark. 
The IceCube 2017 result,?49 represented in black, uses the data sample from Analysis B. The top 
and right plots are the 1-d Ax? profiles of the measured oscillation parameters. 


Fig. 44 shows the 90% allowed region of atmospheric neutrino oscillation pa- 
rameters for the analyses based on Analysis A, and the allowed regions reported by 
other experiments. Overall, the 90% allowed regions for the two analyses are sta- 
tistically consistent, and both results compare favorably with the latest published 
90% contours from other neutrino experiments.°°? °°° 

The shift between the two IceCube contours in both Am}, and sin? 623 can 
be explained by statistical fluctuations alone, as 65% of events in Analysis A are 
unique with respect to Analysis B. Additionally, detailed investigation show that 
differences in the analyses, specifically differences in the parametrization of detector 
effects, in the event selection and in reconstruction, can lead to small systematic 
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shifts of < 0.50 as well. 

While the measured oscillation parameters agree with previous experiments, it 
is important to keep in mind that they have been measured at a roughly ten times 
higher energy. The measurement is therefore also sensitive to any new neutrino 
physics, an important consideration when the precision of the IceCube measure- 
ments will be further improved. 

An upgrade project is underway that will result in the deployment of 7 new 
strings that will further increase the light collection in DeepCore .*°° By decreasing 
the separation between IceCube strings from 125 to 75 meters and the distance be- 
tween light sensors along a string from 17 to 7 meters, DeepCore reduced the energy 
threshold from the TeV-and-above scale suitable for neutrino astrophysics to the 
GeV scale required for measurements of neutrino properties at the highest energies. 
The seven new strings will be deployed with vertical and horizontal spacings that 
are three times smaller than DeepCore and will also include advanced calibration 
devices. 

This new instrumentation will dramatically boost IceCube’s performance at the 
lowest energies, increasing the samples of atmospheric neutrinos by a factor of ten. 
As a result of this upgrade, IceCube will yield the world’s best measurements in 
neutrino oscillations as well as critical measurements that could provide evidence 
for new physics in the neutrino sector. 


9.0.2. Supernova Neutrinos 


A core-collapse supernova radiates the vast majority of the binding energy of the 
resulting compact remnant in the form of neutrinos of all flavors. Information about 
the astrophysics of the collapse and subsequent explosion, and about the physics of 
neutrinos, is encoded in the time, energy, and flavor structure of the neutrino burst. 
When supernova 1987A exploded just a couple of tens of events provided sufficient 
information not only to probe the physics of the explosion, but also to shed new light 
on neutrino properties such as their mass, magnetic moment, and lifetime. Neu- 
trino oscillations produce significant modifications of the neutrino spectra. Matter 
transitions in the expanding remnant have significant time- and energy-dependent 
effects on the spectra of either neutrinos or antineutrinos, depending on the neu- 
trino mixing parameters and the neutrino mass hierarchy. Furthermore, more exotic 
“collective” effects resulting from neutrino-neutrino interactions can result in sig- 
nificant modulation of the spectra, such as “spectral swaps,” which exchange flavor 
spectra above or below a particular energy threshold. 

Although designed as a high-energy Cherenkov detector with a nominal thresh- 
old of 100 GeV (~ 5 GeV in DeepCore), IceCube is still capable of detecting MeV- 
energy neutrinos from core-collapse supernovae by counting them, albeit without 
recording the direction of individual events. A large burst of below-threshold MeV 
neutrinos produced by a supernova will nevertheless be detected as a collective in- 
crease in photomultiplier counting rates on top of their low dark noise in sterile ice. 


February 3, 2022 1:35 ws-rv961x669 Neutrino Physics and Astrophysics main page 98 


98 IceCube Neutrinos 


The origin of this increased rate is Cherenkov light from shower particles produced 
by supernova electron antineutrinos interacting in the ice, predominantly by the 
inverse beta decay reaction. Additionally, a new supernova data acquisition that is 
capable of recording multiple photons from a single neutrino event provides a proxy 
for measuring their energy. 

In ice, positron tracks of about 0.6cm x (£,/MeV) length radiate 178 x 
(Ee+/MeV) Cherenkov photons in the 300-600 nm wavelength range. From the 
approximate Æ? dependence of the neutrino cross section and the linear energy 
dependence of the track length, the light yield per neutrino scales with E?. With 
absorption lengths exceeding 100 m, photons travel long distances in the ice so that 
each DOM effectively monitors several hundred cubic meters of ice. Typically, only 
a single photon from each interaction reaches one of the photomultipliers, which 
are vertically separated by roughly 17 m and horizontally by 125 m. The Deep- 
Core subdetector, equipped with a denser array of high efficiency photomultipliers, 
provides higher detection and coincidence probabilities. Although the rate increase 
in individual light sensors is not statistically significant, the effect will be clearly 
seen once the rise is considered collectively over many independent sensors. For 
a supernova at 10kpc, a signal of one million photoelectros will trace a detailed 
movie of the time evolution of the neutrino signal. The additional measurement 
of the rate of two-photon correlations is sensitive to the energy of the supernova 
neutrinos. With a two megaton effective volume for supernova neutrinos IceCube 
is the most precise detector for analyzing their neutrino light curve." 

While IceCube does not resolve individual neutrino events, it provides a high- 
statistics and detailed movie of the supernova as a function of time. Rather than 
just monitoring the noise in the DOMs, IceCube buffers the information from all 
photomultiplier hits.” Every photon will be recorded to an accuracy of 2 ns in the 
case of a supernova. This low-level data provides several advantages: the complete 
detector information is available and the data-buffered at an early stage of the 
data acquisition system on the so-called string hubs-will be available in the unlikely 
case that the data acquisition fails, for instance in the case of an extremely close 
supernova, which could exhaust the system. The automatized hitspooling has been 
working reliably for several years, including the automatic transfer of data in case 
of serious alarms to the North. The data from these triggered events have been 
carefully studied and compared to the results from the standard supernova data 
acquisition. 

The excellent sensitivity to neutrino properties such as the neutrino hierarchy 
as well as the possibility of detecting the neutronization burst, a short outbreak 
of anti- Ve released by electron capture on protons soon after collapse are discussed 
in.°°” Also, the formation of a black hole or a quark star as well as the characteris- 
tics of shock waves are investigated to illustrate IceCube’s capability for supernova 
detection. In Fig. 45 we illustrate the sensitivity of IceCube to the hierarchy for a 
supernova at the most likely distance of 10 kpc. 


February 3, 2022 1:35 ws-rv961x669 Neutrino Physics and Astrophysics main page 99 


F. Halzen & A. Kheirandish 99 


-+-+ No Oscillations 
— Scenario A (NH) 


Scenario B (IH) 


a m 
ahal 
td Saat 

n a 


DOM Hits (20ms binning) 
(e>) [0>] Ww Q3 Q3 D P 
O S] p= fo?) œ O S] 
oO O oO oO oO oO oO 
O oO O oO oO oO oO 
(=) oO oO oO oO oO oO 


0 0.5 1 1.5 2 2.5 3 
Time Post-Bounce [s] 


Fig. 45. Expected supernova neutrino events as a function of their arrival times for a supernova 
at 10 kpc distance. Scenarios A (normal hierarchy) and B (inverted hierarchy) are contrasted 
while the case for no oscillations is shown as a reference. 


9.0.3. Search for Dark Matter 


A neutrino “telescope” is a powerful tool to search for the particle nature of dark 
matter. Although IceCube, was primarily designed as a discovery instrument, 
searching for dark matter was from inception one of its priority missions. Its Deep- 
Core infill that lowers IceCube’s threshold over a significant part of the volume was 
initially proposed by Per Olof Hulth as a way to enhance IceCube’s capabilities for 
detecting dark matter particles with lower masses. In this context it is worthy of 
note that the AMANDA detector, the forerunner and proof of concept for IceCube, 
received a significant fraction of its initial funding from the Berkeley Center for 
Particle Astrophysics to search for dark matter. 

IceCube searches indirectly for dark matter by looking for neutrinos from con- 
centrations of weakly interacting massive particles (WIMPs) in the Sun,*°? °°! in 
the Milky Way*°?-°° and in nearby galaxies.°°’ The neutrinos are secondary prod- 
ucts of the annihilation of pairs of WIMPs into standard Model particles which 
includes decays into neutrinos, either directly or via their leptonic decays. 

If “weakly” interacting massive particles (WIMPs) make up the dark matter, 
they have been swept up by the Sun for billions of years as the Solar System moves 
about the galactic halo. Though interacting weakly, they will occasionally scatter 
elastically with nuclei in the Sun and lose enough momentum to become gravita- 


February 3, 2022 1:35 ws-rv961x669 Neutrino Physics and Astrophysics main page 100 


100 IceCube Neutrinos 


tionally bound. Over the lifetime of the Sun, a sufficient density of WIMPs may 
accumulate in its center so that an equilibrium is established between their capture 
and annihilation. The annihilation products of these WIMPs represent an indirect 
signature of halo dark matter, their presence revealed by neutrinos which escape the 
Sun with minimal absorption. The neutrinos are, for instance, the decay products 
of heavy quarks and weak bosons resulting from the annihilation of WIMPs into 
XX > TTT , bb or WtW-. Neutrino “telescopes” are sensitive to such neutrinos 
because of their relatively large neutrino energy reflecting the mass of the parent 
WIMP. 

The beauty of the indirect detection technique using neutrinos originating in 
the sun is that the astrophysics of the problem is understood. The source in the 
sun has built up over solar time sampling the dark matter throughout the galaxy; 
therefore, any unanticipated structure in the halo has been averaged out over time. 
Other astrophysical uncertainties, such as the dark matter velocity distribution, 


368,369 Given 


have little or no impact on the capture rate of dark matter in the sun. 
a WIMP mass and cross section (and the assumption that the dark matter is not 
exclusively matter or antimatter), one can unambiguously predict the signal in a 
neutrino telescope. If not observed, the model will be ruled out. This is in contrast 
to indirect searches for photons from WIMP annihilation, whose sensitivity depends 
critically on the structure of halo dark matter; observation requires cuspy structure 
near the galactic center or clustering on appropriate scales elsewhere; observation 
not only necessitates appropriate WIMP properties, but also favorable astrophysical 
circumstances. 

Extensions of the Standard Model of quarks and leptons, required to solve the 
hierarchy problem, naturally yield dark matter candidates. For instance, the neu- 
tralino, the lightest stable particle in most supersymmetric models, has been in- 
tensively studied as a possible dark matter candidate. Detecting it has become 
the benchmark by which experiments are evaluated and mutually compared. We 
will follow this tradition here. Supersymmetric models allow for a large number of 
free parameters and, unfortunately, for a variety of parameter sets that are able to 
generate the observed dark matter density in the context of standard big bang cos- 
mology. The neutralino interacts with ordinary matter by spin-independent (e.g., 
Higgs exchange) and by spin-dependent (e.g., Z-boson exchange) interactions. The 
first mechanism favors direct detection experiments?” because the WIMP inter- 
acts coherently, resulting in an increase in sensitivity proportional to the square of 
the atomic number of the detector material. Given the rapid improvement of the 
sensitivity of direct experiments, indirect detection experiments are therefore not 
competitive for models where the dark matter has a large spin-independent cou- 
pling. This is not the case for spin-dependent models for which IceCube achieves 
world-best limits. Within the context of supersymmetry, direct and indirect exper- 
iments are complementary. 

Although IceCube detects neutrinos of all flavors, sensitivity to neutrinos pro- 
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duced by WIMPs in the sun is primarily achieved by exploiting the degree accuracy 
with which muon neutrinos can be pointed back to the Sun. Searches optimized 
to probe WIMP masses below about 50 GeV benefit from extending the detection 
channel to all neutrino flavors. 

IceCube is most sensitive to WIMPs with significant spin-dependent interactions 
with protons because they result in large concentrations of WIMPS in the Sun, a 
nearby and readily identifiable source of protons. The signal would simply be an 
excess of neutrinos of GeV energy and above, depending on the WIMP mass, from 
the direction of the Sun over the atmospheric neutrino background in the same 
angular window. There would be no alternative astrophysical explanation of such 
a signal that represents a smoking gun for dark matter particles. In most WIMP 
scenarios, the cross section for WIMP capture (oy p) is large enough so that an 
equilibrium between capture and annihilation would have been achieved within the 
age of the solar system.°’! In this case, limits on neutrinos from the Sun can be 
expressed in terms of the capture cross section, Cy, p. If equilibrium is not reached, 
weaker limits can still be derived. 

Quantitative interpretation of the IceCube limits depends on detailed analysis 
of capture rates and annihilation channels. The analysis traditionally uses Dark- 
SUSY,°”! which builds on the classic calculations of Gould??? for capture and anni- 
hilation rates. The neutrino spectrum from annihilation of pairs of WIMPs depends 
on the dominant channel for coupling to standard model particles. The neutrino 
spectrum at production is also modified by subsequent interactions and oscillations 
of neutrinos as they propagate out from the solar core. Three flavor oscillations 
with matter effects are included as in reference 373. The range of possibilities is 
bracketed by calculating two extremes, 


Xx — bb (soft) and xx > r*77 (hard). (133) 


In recent works T+T~ serves typically as hard channel, in older works often a com- 
bination of W= and T= are used, where t~ is replaced with WF for my > my. 
The world-best limits on the spin-dependent interaction cross section of WIMPs 
with normal matter are shown in Fig. 46. Even stronger limits, or a discovery, are 


forthcoming using an 8-year data set and improved reconstructions. 


9.1. IceCube, the Facility 


During its construction phase, IceCube demonstrated a significant potential for 
facilitating a range of other research. For example, a dust logger provided mea- 
surements with millimeter precision that are valuable for event reconstruction in 
IceCube but that also provide a record of surface winds over more than 100,000 
years.374 

Already during construction of AMANDA, antennas forming the RICE detector 
were deployed in deep holes to expand the target volume in the search for cosmogenic 


375 


neutrinos.’‘? An acoustic test setup of receivers in the upper portion of four IceCube 
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Fig. 46. Limits on the spin-dependent WIMP-proton cross section (o ) from IceCube 3-year 


361 


holes was deployed in 2007 to explore the acoustic technique for detecting ultra- 
high-energy neutrinos. A retrievable transmitter (pinger) was submerged briefly 
in several newly prepared holes at various depths and distances from the receivers 
to measure the attenuation of sound in ice. The attenuation length of 300 m is 
significantly less than had been expected.?”° In contrast, the Askaryan Radio Array 
(ARA)? is taking advantage of the kilometer-scale attenuation for radio signals 
in ice to operate a detector that performs R&D for a much larger array that will 
be colocated with the next-generation IceCube detector. The ARA detectors send 
data to computers housed in the IceCube Lab (ICL) for staging and transmission 
to the north. 

The DM-Ice experiment?’* proposes to repeat the DAMA experiment in the 
quiet and sterile environment of the Antarctic ice. An interesting feature of the 
measurement is the fact that the seasonal modulation of the muon rate has the 
opposite phase relative to the motion of the Earth through the halo of dark matter 
compared to a detector in the Northern Hemisphere. A test detector to explore the 
noise environment for DM-Ice was deployed at the bottom of an IceCube string in 
December 2010. Its computers and data transmission are also hosted in the ICL. 
The experiment will be deployed in the 2022-23 Antarctic summer as part of the 
Upgrade project, which we will discuss next. 
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The enhancement of the low-energy capabilities of IceCube provided by the 
DeepCore subarray had led to the PINGU proposal?” to deploy an additional 40 
strings within the existing DeepCore detector. This aimed to lower the threshold 
below 5 GeV (< 25 m muon track length in ice). In this energy range, matter effects 
in the Earth lead to resonant oscillations of V, + Ve (D, + De) for normal (inverted) 
hierarchy**° that depend on zenith angle. By taking advantage of the fact that the 
neutrino cross section is larger than that for antineutrinos, coupled with the excess 
of v, compared to y,, a measurement sensitive to the neutrino mass hierarchy is 
possible on a relatively short timescale. PINGU would also have had sensitivity to 
v,, disappearance, v, appearance, and maximal mixing. The lower energy threshold 
would enhance the indirect searches for dark matter with IceCube as well as the 
sensitivity to neutrinos from supernova explosions. In addition, PINGU presented 
the potential for neutrino tomography of the Earth. The science goals of PINGU 
are now pursued with the IceCube Upgrade project. Although the construction of 
a PINGU-style detector can be launched in the future, the next-generation IceCube 
detector will target high-energy neutrino astronomy instead. The ORCA detector, 
with the same science goals and performance, is now under construction in the 
Mediterranean Sea.°°! 

The so-called IceCube Upgrade is also under construction. It will deploy seven 
new strings at the bottom of the detector array that have been designed as an 
incremental extension of the DeepCore detector and as a test bed for the technology 
of the next-generation detector. 

The new strings will increase the light collection in what is already a denser 
region of the gigaton neutrino detector. By decreasing the separation between 
IceCube strings from 125 to 75 meters and the distance between light sensors along 
a string from 17 to 7 meters, DeepCore reduced the energy threshold from the TeV- 
and-above scale suitable for neutrino astrophysics to the GeV scale required for 
measurements of neutrino properties at the highest energies. The seven new strings 
will be deployed with vertical and horizontal spacings that are three times smaller 
than DeepCore and will also include advanced calibration devices. 

This new instrumentation will dramatically boost IceCube’s performance at the 
lowest energies, increasing the samples of atmospheric neutrinos by a factor of ten, 
and will enhance the pointing resolution of astrophysical neutrinos. As a result of 
this upgrade, IceCube will yield the world’s best measurements in neutrino oscilla- 
tions as well as critical measurements that could provide evidence for new physics 
in the neutrino sector. 

The new calibration devices will advance our understanding of the response vari- 
ability of the light sensors in both current and new strings, resulting in enhanced 
reconstruction of cascade events and better identification of tau neutrinos. Cascade 
signatures are the light patterns of more than 75% of astrophysical neutrinos and 
currently have a poorer pointing resolution than the so-called tracks—cascades allow 
reconstructing the direction of the incoming neutrino with a typical angular resolu- 
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tion of 5-10 degrees, while tracks point to the neutrino’s origin to within less than 
1 degree. The refined calibration of the existing sensors will also enable a reanaly- 
sis of more than ten years of archival data and increase the discovery potential of 
neutrino sources. 


9.2. From Discovery to Astronomy, and More... 


Accelerators of cosmic rays produce neutrino fluxes limited in energy to roughly 
5% of the maximal energy of the protons or nuclei. For Galactic neutrino sources, 
we expect neutrino spectra with a cutoff of a few hundred TeV. Detection of these 
neutrinos requires optimized sensitivities in the TeV range. At these energies, the 
atmospheric muon background limits the field of view of neutrino telescopes to the 
downward hemisphere. With IceCube focusing on high energies, a second kilometer- 
scale neutrino telescope in the Northern Hemisphere would ideally be optimized to 
observe the Galactic center and the largest part of the Galactic plane. 

Following the pioneering work of DUMAND,® several neutrino telescope projects 
were initiated in the Mediterranean in the 1990s.1°-!? In 2008, the construction of 
the ANTARES detector off the coast of France was completed. With an instru- 
mented volume at a few percent of a cubic kilometer, ANTARES reaches roughly 
the same sensitivity as AMANDA and is currently the most sensitive observatory 
for high-energy neutrinos in the Northern Hemisphere. It has demonstrated the fea- 
sibility of neutrino detection in the deep sea and has provided a wealth of technical 
experience and design solutions for deep-sea components. 

While less sensitive than IceCube to a diffuse extragalactic neutrino flux, 
ANTARES has demonstrated its competitive sensitivity to neutrino emission from 
the Galactic center?*?3°% 
sphere.**+ The important synergy between Mediterranean and Antarctic neutrino 
telescopes has been demonstrated by the first joint study of continuous neutrino 
sources?°? 396 

An international collaboration has started construction of a multi-cubic- 
kilometer neutrino telescope in the Mediterranean Sea, KM3NeT.**! Major progress 
has been made in establishing the reliability and the cost-effectiveness of the de- 
sign. Most importantly, this includes the development of a digital optical module 
that incorporates 31 3-inch photomultipliers instead of one large photomultiplier 
tube; see Fig. 47. The advantages are a tripling of the photocathode area per opti- 
cal module, a segmentation of the photocathode allowing for a clean identification 
of coincident Cherenkov photons, some directional sensitivity, and a reduction of 
the overall number of penetrators and connectors, which are expensive and failure- 
prone. For all photomultiplier signals exceeding the noise level, time-over-threshold 
information is digitized and time-stamped by electronic modules housed inside the 
optical modules. This information is sent via optical fibers to shore, where the data 
stream will be filtered online for event candidates. 

KM3NerT in its second phase?! will consist of two units for astrophysical neu- 


and extragalactic neutrino sources in the Southern Hemi- 


as well as neutrino follow-up campaigns of gravitational waves. 
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Fig. 47. The KM3NeT optical module (from reference 381). The optical module consists of a 
glass sphere with a diameter of 42 cm, housing 31 photosensors (yellowish disks). The glass sphere 
can withstand the pressure of the water and is transparent to the faint light that must be detected 
to see neutrinos. 


trino observations, each consisting of 115 strings (detection units) carrying more 
than 2,000 optical modules, and one ORCA detector studying fundamental neu- 
trino physics with atmospheric neutrinos. The detection units are anchored to the 
seabed with deadweights and kept vertical by submerged buoys. The vertical dis- 
tances between optical modules will be 36 meters, with horizontal distances between 
detection units at about 90 meters. Construction is now ongoing near Capo Passero 
(east of Sicily). 

A parallel effort is underway in Lake Baikal with the construction of the deep 
underwater neutrino telescope Baikal-GVD (Gigaton Volume Detector).°°” The 
first GVD cluster was upgraded in spring 2016 to its final size (288 optical modules, 
a geometry of 120 meters in diameter, 525 meters high, and instrumented volume 
of 6 Mton). Each of the eight strings consists of three sections with 12 optical 
modules. At present 7 of the 14 clusters have been deployed reaching a sensitivity 
close to the diffuse cosmic neutrinos flux observed by IceCube. 

Further progress requires larger instruments. IceCube therefore proposes as a 
next step to capitalize on its large absorption length to instrument 10 km? of glacial 
ice at the South Pole and thereby improving on IceCube’s sensitive volume by 
an order of magnitude.®38:388:389 This large gain is made possible by the unique 
optical properties of the Antarctic glacier revealed by the construction of IceCube. 
As a consequence of the extremely long photon absorption lengths in the deep 
Antarctic ice, the spacing between strings of light sensors can be increased from 125 
to close to 250 meters without significant loss of performance of the instrument. The 
instrumented volume can therefore grow by one order of magnitude while keeping 
the construction budget of a next-generation instrument at the level of the cost 
of the current IceCube detector. The new facility will increase the event rates of 
cosmic events from hundreds to thousands over several years. The superior angular 
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resolution of the longer muon tracks will allow for the discovery of the sources seen at 
the 30-level in the 10-year sky map. Construction of a next-generation instrument 
with at least five times higher sensitivity would likely result in the observation of 
cosmogenic neutrinos.'°? Obviously, higher sensitivity would also benefit the wide 
range of measurements performed with the present detector, from the search for 
dark matter to the precision limits on any violation of Lorentz invariance. 
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